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Virtually  all  the  literature  concerning  queueing  and  storage  theory 
assumes  that  the  input  to  a system  has  an  arrival  rate  which  is  constant 
over  time.  This  research  is  concerned  with  extending  results  to  the 
more  realistic  case  of  nonhomogeneous  input.  It  is  shown  that  under  a 
variety  of  circumstances,  the  presence  of  a periodic  pattern  of  input 
is  enough  to  insure  that  a quasi-limiting  distribution  of  store  content 
exists.  For  stores  with  nonhomogeneous  Compound  Poisson  input  and 
general  release  rule  a probabilistic  approach  is  developed  which  allows 
equations  to  be  written  down  easily  for  such  quantities  as  the  mean 
store  content  and  the  probability  of  an  empty  store.  To  do  this  it  is 
assumed  that  the  intensity  can  be  written  in  the  form  X(t)=Xo(l+e$(t)). 
The  formulas  are  series  expansions  in  e.  A more  in-depth  look  at  the 


first  two  terms  of  the  e-expansion  shows  that  for  periodic  <t>(t),  a 
general  expression  can  be  found  for  an  approximation  to  the  probability 
of  an  empty  store.  Applications  of  the  formulas  for  stores  with  release 
rules  r(x)=c,  r(x)=a+bx  and  r(x)=ax  are  examined.  For  at  least  some 
instances  of  the  M( t ) /M/1  queue  it  is  shown  that  the  approximations 
developed  are  very  good.  The  problem  of  a stochastic  intensity  function 

is  also  briefly  considered.  
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Equations  will  be  numbered  consecutively  throughout  each  chapter. 
References  to  equations  within  a chapter  will  be  by  number  only. 
References  to  equations  in  earlier  chapters  will  include  both  chapter 
and  equation  number.  Thus  a reference  to  equation  25  in  Chapter  3 would 
be  to  equation  3.25. 

Notation 

The  Laplace-Stieltjes  Transform  of  a function  F will  be  denoted  by 
F*(s) ; i.e.  F*(s)  = £°VsxdF(x) . 

The  (ordinary)  Laplace  Transform  of  a function  F will  be  denoted 
F°(x);  i.e.  F°(s)  = £°°e"sxF(x)dx. 

The  Stieltjes  Convolution  of  functions  F and  G will  be  denoted 
F*G(x)i  ’•••  F-G{x)  = ^ F(x-u)dG(u)  - ^G(x-u)dF(u). 

We  note  that  the  above  functions  F and  G are  assumed  to  be  zero  on 
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CHAPTER  1 

INTRODUCTION  AND  REVIEW  OF  THE  LITERATURE 

The  problem  this  research  considers  is  the  modeling  of  and  effect  on 
queueing  and  storage  systems  which  have  input  that  is  not  homogeneous  in 
time.  Most  work  concerning  queueing  and  storage  processes  has  been  for 
the  case  where  the  times  between  inputs  to  the  system  are  i.i.d.  random 
variables.  This  implies  that  the  probability  of  an  input  in  one  time 
interval  is  equal  to  the  probability  of  an  input  in  any  other  time  inter- 
val of  the  same  length.  But  in  many  practical  situations  the  input 
process  is  of  a more  general  nature,  often  varying  with  time.  Sometimes 
the  arrival  rate  seems  to  fluctuate  randomly  over  time.  More  often, 
though,  the  dependence  on  time  is  of  a more  deterministic  nature,  with 
the  arrival  frequency  depending  on,  say,  the  time  of  day,  month  or  year. 

Automobile  traffic,  airline  departures  and  landings,  restaurant  and  resort 
business,  telephones— all  these  exhibit  a cyclical  or  periodic  pattern  of 
input  and  are  quite  typical  of  the  types  of  processes  one  often  finds. 

A common  method  that  has  been  used  to  deal  with  the  problem  of  time 
varying  input  (when  it  is  recognized)  is  to  treat  periods  of  widely 
differing  inputs  separately.  For  example,  in  queueing  theory  much  has 
been  written  on  the  problem  of  the  rush  hour-referred  to  as  heavy  traffic 
and  non-stationary  queues.  Such  work  focuses  on  the  time  when  the  rate 
of  arrival  to  a queue  is  almost  equal  to  or  even  greater  than  the  ser- 
vice rate--not  the  typical  situation  the  majority  of  the  time.  If 
reasonably  practical  results  could  be  obtained  for  a process  with 
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nonhomogeneous  input,  then  we  could  model  the  entire  system  simultane- 
ously, which  would  certainly  be  preferable.  To  understand  the  complex- 
ity of  the  problem,  we  should  first  look  at  what  has  been  discovered  in 
the  literature  to  date. 

The  literature  of  queueing  theory  has  developed  quite  separately  from 
that  of  storage  and  dam  theory,  although  the  two  areas  are  often  differ- 
ent aspects  of  the  same  problem.  This  can  readily  be  seen  if  we  con- 
sider that  the  input  to  a store  (or  dam)  in  (0,t]  is  equivalent  to  the 
total  service  time,  or  work  load,  of  all  customers  to  arrive  in  ( 0 , t ] ; 
the  content  of  a store  at  time  t corresponds  to  the  (virtual)  waiting 
time  of  a customer  who  arrives  at  a queue  at  time  t,  and  the  release 
rate  of  the  content  of  a store,  per  unit  time,  is  equivalent  to  the 
amount  of  work  load  processed  per  unit  time  by  the  server  of  a queue. 
Thus,  in  many  cases,  what  has  been  discovered  about  one  process  can  be 
applied  to  the  other.  Still,  the  work  done  in  queueing  theory  on  non- 
homogeneous arrivals  has  not  been  picked  up  by  those  concerned  with 
storage  theory.  There  has  been  relatively  little  published  in  the  area 
of  nonhomogeneous  processes  and  virtually  all  that  has  been  done  has 
been  in  the  area  of  queues. 

The  first  major  work  in  this  area  was  done  by  Taka'cs  (1955)  who 
investigated  the  virtual  waiting  time  process  of  a single  server  queue. 
The  virtual  waiting  time  is  defined  as  the  time  a customer  arriving 
would  have  to  wait  to  be  served.  If  the  system  is  empty  at  the  time  of 
arrival  then  the  customer  is  served  immediately  and  the  waiting  time  is 
zero. 

For  the  single  server  queue  with  nonhomogeneous  Poisson  arrivals 
(intensity  \(t))  and  general  service  time  distribution,  Taka'cs  derived 
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the  integro-differential  equation  for  the  distribution  of  waiting  time 
at  t,  F(t,x): 

i 

3P(tiX)  - 3F(t,x)  - x(t)F(t,x)+X(t)yxH(x-y)dvF(t,y). 

3t  3x  o y 

From  this  he  derived  an  expression  for  <t>(t,s),  the  Laplace-Stieltjes 
Transform  of  the  waiting  time: 

st  - [l-ii'(s) ]A( t)  f r1  -su+[l-<Ms)]A(u)  "I 
$(t,s)  = e [l-sjoe  F(u,o)duJ  . 

where  i>{s)  is  the  Laplace-Stieltjes  Transform  of  the  service  distribution 
and  A(t)  = /*A(u)du. 

If  the  probability  the  system  is  empty  c?:n  be  determined  uniquely, 
then  the  distribution  of  waiting  time  can  be  determined.  As  Takacs  ad- 
mits, the  determination  of  this  probability  is  generally  quite  difficult. 

In  the  particular  case  where  A(t)  converges  to  a positive  constant  A for 
large  values  of  t,  and  Acx<l  (where  a is  the  expected  service  time)  the 
limiting  distribution  of  waiting  time  exists,  independent  of  the  initial 
distribution.  The  probability  of  emptiness  in  this  case  is  1-Aa,  which 
agrees  with  the  results  for  the  queue  with  homogeneous  arrival  rate  A. 

The  problem  of  determining  the  probability  of  emptiness  has  been  pur- 
sued by  a few  authors  with  varying  degrees  of  success.  Reich  (1958, 

1959)  was  able  to  show  that  the  probability  a queue  is  empty  at  time  t 
is  the  unique  solution  of  a Vol terra  equation  of  the  first  kind.  With 
this  he  was  able  to  examine  the  behavior  of  the  probability  for  large  t, 
culminating  in  a result  which  gives  the  time  average  of  the  probability 


subject  to  an  error  of  magnitude  o(l). 
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Gani  (1962)  approached  the  problem  in  the  context  of  dam  processes. 
He  showed  that  the  probability  of  emptiness  is  a function  of  the  first 
emptiness  probability,  which  in  the  case  of  inputs  of  constant  unit  size 
reduces  to  a simple  recurrence  relationship. 

Hasofer  (1964,  1965)  was  able  to  use  the  results  of  Reich  to  get 
explicit  results  for  the  M(t)/G/1  queue  by  assuming  a particular  form 
for  the  Poisson  parameter  \(t) . Taking  /H(u)du-t-zb(t)  Hasofer  was 
able  to  write  the  probability  of  an  empty  queue  at  time  t as  a power 

CO 

series  in  z,  P(x,t)=  n|  ZnFn(t).  Using  Reich's  equations,  some  intri- 
cate complex  analysis  and  the  additional  assumption  that  b(t)  and 
b'(t)  are  uniformly  bounded,  Hasofer  was  able  to  find  a general  expres- 
sion for  the  Laplace  Transform  of  the  F : 

n 

l [j_  rx+i°° 

Fn°(p)=(5o(p)  y71’  1 bn(p-o)gn(o)do 

L Jx-i°° 


n 

k= 1 r=o 


I„  ?rf 


<+i» 

bk-Fp-s)Gk<s>  2W 


,x+i«> 


b°(s-o)F°  (o)dods> 
rv  c n-k  ' 


'x-i^ 


1 x-i» 


where 


G°(p)=  UkMnll";  9°(p)=i.  (n)ri-»(n)]n;  b°(p)=  /e'pt[b(t)]ndt 
nW(n)]  n!  [W(n)] 


and  n is  the  unique  solution  of  s-l+iHs)-p=0. 


Much  neater  results  were  obtained  when  the  Poisson  parameter  was 

n 

assumed  periodic.  For  A(t)=t-z  n?1arisin(nu)t+4>r|)  Hasofer  showed  that  the 
functions  Fn  have  the  asymptotic  form 

oo 

Fn(t)  - [Akn  cos(ku)t)+Bknsin(ku>t)] 

and  the  Laplace-Stieltjes  transform  of  the  waiting  time  has  a similar 
asymptotic  form.  The  coefficients  in  the  Fourier  expansion  were  not 
given  explicitly,  although  for  the  special  case  A(t)=l-zwcos(ojt)  Hasofer 
derived  FQ  and  Fj. 

A number  of  authors  have  approached  the  problem  of  nonhomogeneous 
arrivals  by  employing  the  familiar  Kolmogorov  difference-differential 
equations  for  the  queue  size  distribution.  Clark  (1956)  explored  the 
queue  size  distribution  of  a single  server  queue  with  Poisson  arrivals 
and  negative  exponential  service  times,  where  both  parameters  are  con- 
tinuous functions  of  time.  While  no  exact  represen+'ati  m for  the  dis- 
tribution of  queue  size  was  obtained,  the  problem  was  reduced  to  finding 
the  solution  of  an  integral  equation.  Much  simpler  results  were  pre- 
sented for  the  mean  and  variance  of  the  queue  size,  which  only  depend  on 
the  (unknown)  probability  of  emptiness.  In  the  special  case  where 
traffic  intensity  is  constant  the  distribution  of  queue  size  was  found 
explicitly.  It  can  readily  be  seen  that  in  this  case  the  queue  is 
equivalent  to  one  with  constant  arrival  and  service  rates. 

Luchak  (1956)  treated  a similar  problem,  but  with  a more  general 
service  distribution.  He  approximated  a general  service  distribution 
by  a Weighted-sum  Erlang  distribution.  Starting  again  from  the  Kolmogorov 
equations,  Luchak  found  that  if  the  traffic  intensity  could  be  expressed 
as  a polynomial  in  time,  then  the  distribution  of  queue  size  could  be 


found  recursively,  but  only  in  selective  cases  could  a closed  form  be 
found.  At  the  end  of  his  article,  Luchak  noted  that  periodic  arrival 
rates  could  be  dealt  with  more  easily  than  most,  since  the  queue  size 
distribution  need  only  be  found  for  the  initial  period.  The  state  at 
the  end  of  each  period  can  then  be  used  as  the  initial  condition  for  the 
next  period.  He  suggested  that  if  the  initial  state  of  the  system  were 
taken  to  be  the  steady  state  solution  of  a queue  with  constant  traffic 
intensity  equal  to  the  average  intensity  over  a period,  then  the  "quasi- 
stationary" steady  state  would  be  reached  in  only  a few  periods.  We 
will  make  use  of  this  suggestion  in  later  chapters. 

In  recent  years  the  work  of  Luchak  and  Clarke  has  been  expanded  to 
include  more  generality.  Koopman  (1972)  was  interested  in  arrivals  and 
departures  of  airplanes  at  an  air  terminal.  He  incorporated  certain  of 
the  peculiarities  of  the  airport  into  his  model.  These  included  periodi- 
city (24  hours)  and  a limit  on  the  maximum  queue  size.  He  also  made 
arrival  and  service  time  parameters  dependent  both  on  time  and  queue 
size.  The  method  of  solution  suggested  by  Koopman  was  a direct  computer 
solution  of  the  (finite)  system  of  differential  equations. 

Moore  (1975)  like  Luchak,  used  a Weighted-sum  Erlang  distribution 
for  service  times  and  assumed  customers  could  arrive  in  groups.  He 
developed  a technique,  in  the  stationary  case,  which  focused  on  the 
embedded  Markov  chain  formed  by  looking  at  departure  times.  Then  by 
approximating  the  continuously  varying  arrival  rate  by  a step  function 
he  was  able  to  extend  this,  developing  a computer  algorithm  solution. 

E.  L.  Leese  (1967)  wrote  a paper  in  which  he  evaluated  the  computa- 
tional feasability  of  a number  of  proposed  solutions  for  the  nonhomo- 
geneous  queue.  Leese  considered  direct  solution  of  the  Kolmogorov 
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differential -difference  equations,  generating  function  techniques,  a 
Taylor  series  expansion  of  the  queue  size  probabilities,  a matrix  approach 
which  uses  a step  function  approximation  to  the  traffic  intensity,  and  a 
method  developed  by  Wragg  (1963)  which  involves  solving  an  integral  equa- 
tion. Most  of  the  objections  raised  by  Leese  revolve  around  the  excessive 
number  of  computational  steps  required.  It  should  be  noted  that  over  ten 
years  of  rapid  advances  in  computer  technology  has  abrogated  some  of  his 
original  objections.  In  fact,  we  will  later  investigate  the  uses  of  the 
matrix  approach  when  the  arrival  pattern  is  periodic. 

A number  of  authors  have  investigated  a queueing  problem  where  the 
arrival  (and  possibly  service)  rate  is  not  a deterministic  function  of 
time.  Rather,  the  arrival  rate  is  a heterogeneous  Poisson  process 
governed  by  an  "extraneous  phase  process"  which  is  a continuous  time 
Markov  chain,  i.e.  customers  arrive  in  a Poisson  stream  with  the  arrival 
rate  depending  on  the  "phase"  of  the  queue,  and  the  amount  of  time  spent 
in  a phase  is  a negative  exponential  random  variable.  Service  times  may 
or  may  not  vary  with  each  phase. 

Eisen  and  Tainiter  (1963)  and  Yechiali  and  Naor  (1971)  looked  at  such 
a process  with  two  phases,  where  service  also  followed  a negative  expo- 
nential distribution.  Eisen  and  Tainiter  found  a formula  for  the  mean 
number  in  the  system  and  mean  waiting  time  in  steady  state.  Yechiali 
and  Naor  looked  at  a number  of  special  cases.  They  found  that  if  the 
traffic  intensity  in  each  phase  of  the  queue  is  the  same,  then  the  dis- 
tribution of  queue  size  in  the  steady  state  is  what  would  be  expected  in 
the  M/M/1  queue.  They  also  claimed  the  above  result  would  hold  only  in 
the  case  of  equal  traffic  intensities,  which  we  will  show  later  is  not 
true. 
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Yechiali  (1973)  later  extended  his  result  to  include  an  m-phase  pro- 
cess. To  find  explicit  results  requires  finding  the  roots  of  a difficult 
equation  which  yields  that  nemesis,  the  probability  of  an  empty  queue. 

No  closed  form  solutions  were  presented  except  in  the  special  case  where 
traffic  intensity  of  each  phase  are  all  equal,  with  results  as  before 
for  m=2. 

Neuts  (1971)  extended  this  work  even  further  by  allowing  general  ser- 
vice time  distributions  in  each  phase.  By  using  very  powerful  tools 
along  with  a Markov  renewal  branching  process  approach  Neuts  was  able 
to  find  both  transient  and  steady  state  results. 

Three  separate  papers  have  been  published  concerning  the  existence 
of  limiting  distributions  when  arrivals  to  a queue  are  from  a nonhomo- 
geneous  Poisson  process--one  by  Sakr  (1971),  another  by  Yevdokimova 
(1974)  and  a third  just  recently  published  by  Harrison  and  Lemoine  (1977). 
All  three  assume  that  the  intensity  of  the  arrival  process,  \(t),  is 
periodic,  i.e.  A(t+x)  = \(t)  for  some  period  t and  every  t. 

Sakr  proved  that  for  a queue  with  m servers,  if  the  traffic  intensity 
per  period,  A/p  is  less  than  one  (where  A = - /?\(u)du  is  the  averaqe 

T O 

intensity  per  period  and  1/p  is  the  mean  service  time)  then  the  distri- 
bution of  queue  length  converges  to  a "quasi -stationary"  limiting  dis- 
tribution, where  the  distribution  of  queue  length  is  a function  of  the 
time  within  the  period.  Yevdokimova  expanded  upon  the  work  of  Sakr  to 
give  a representation  of  the  waiting  distribution  for  the  single  server 
queue.  He  showed  that  if  A/u<l  then  the  waiting  time  Is  the  sum  of  two 
functions,  one  of  which  is  periodic  and  independent  of  the  initial  dis- 
tribution, while  the  other  goes  to  zero  as  t goes  to  infinity. 


! 
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The  methods  used  In  both  articles  are  virtually  identical.  Both 
authors  used  two  queues  with  related  arrival  processes  to  bound  the  orig- 
inal process.  They  were  then  able  to  use  classical  theory  for  the  homo- 
geneous queue  to  get  limiting  results.  To  get  his  more  explicit  results, 
Yevdokimova  also  made  use  of  the  fact  that  the  time  points  {nT}^  such 
that  the  system  is  empty  at  nt  form  an  embedded  renewal  process. 

Harrison  and  Lemoine,  with  assumptions  identical  to  those  of  the  first 
two  articles,  proved  the  existence  of  both  time  average  limits  and  limits 
in  distribution  for  the  virtual  waiting  time  and  waiting  time  process. 

They  also  derived  relationships  between  the  limiting  distributions.  The 
method  of  proof  used  is  similar  to  that  of  Yevdokimova,  focusing  on  empty 
times  as  a regenerative  event. 

As  this  brief  literature  review  demonstrates,  very  little  progress 
has  been  made  in  developing  practical  formulae  for  determining  virtual 
waiting  time  when  nonhomogeneous  input  is  considered.  In  fact,  no  work 
at  all  has  been  done  for  storage  processes  with  time  varying  input  and 
general  release  rules  beyond  the  queueing  theory  results  corresponding 
to  r(x)=c.  This  research  will  attempt  to  develop  practical  formulae 
for  the  probability  of  emptiness  and  expected  content  of  stores  with  non- 
homogeneous Compound  Poisson  input  and  general  release  rate. 

In  Chapter  2 we  are  concerned  with  the  existence  of  the  limiting  or 
quasi-limiting  distribution  of  store  content  under  rather  general  condi- 
tions. In  particular,  we  show  that  under  a variety  of  circumstances  the 
presence  of  a periodic  pattern  of  input  is  enough  to  ensure  that  a limit 
exists. 

A probabilistic  approach  is  developed  in  Chapter  3 which  allows  us 
to  quickly  and  easily  write  down  equations  for  such  quantities  as  the 
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distribution  of  store  content,  probability  of  emptiness,  etc.  To  do 
this  we  assume  that  the  intensity  has  the  form  X(t)=XQ( l+r<(>(t) ) . The 
formulae  are  series  expansions  in  e.  A more  in  depth  look  at  the  first 
two  terms  of  the  expansion  shows  that  for  periodic  <j>(t),  a general 
expression  can  be  found  for  an  approximation  to  the  probability  of  an 
empty  store. 

Chapters  4 and  5 focus  on  applications  of  the  formulae  developed  in 
Chapter  3 for  specific  release  rules.  Chapter  4 treats  the  most  common 
rule,  r(x)=l,x>0.  Specific  formulae  are  given  for  an  approximation  to 
the  mean  store  content  and  probability  of  an  empty  store  when  the  input 
size  distribution  is  negative  exponential  or  Weighted-sum  Erlang. 

Chapter  5 considers  the  release  rules  r(x)=a+bx  and  r(x)=cx,x>0.  The 
latter  is  representative  of  a class  of  storage  processes  which  never 
empty.  We  show  that  the  expect’d  content  of  such  a store  can  be  written 
as  a second  degree  polynomial  in  e. 

Chapter  6 presents  an  alternative  matrix  approach  to  the  study  of 
the  M/M/1  queue  with  nonhomogeneous  input.  Its  particular  appeal  in  the 
case  of  periodic  input  is  pursued.  We  also  use  this  method  as  a means 
to  investigate  the  adequacy  of  the  c-approximations  in  the  special  case 
of  negative  exponential  service  time. 

In  Chapter  7 we  consider  a slightly  different  problem  from  earlier 
chapters.  There  we  are  concerned  with  the  situation  where  X(t)=X(t,<i\) 
is  a stochastic  process.  We  show  that  it  is  sufficient  to  know  the  mean 
function  of  X(t,<.»)  to  be  able  to  investigate  the  behavior  of  the  store 
content.  The  problem  presented  by  Yechiali  and  Naor  (1971)  is  also 


examined,  with  the  surprising  result  that  when  the  mean  service  time 
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does  not  vary,  the  formula  for  the  probability  of  emptiness  is  quite 
simple  and  only  involves  the  first  two  terms  of  the  e-expansion.  We 
also  correct  a mistake  by  Yechiali  in  his  papers. 


CHAPTER  2 


LIMIT  THEOREMS  FOR  STORAGE  SYSTEMS 

In  this  chapter  we  will  prove  the  existence  of  certain  limiting 
probabilities  for  storage  systems  with  nonhomogeneous  input.  We  will 
be  concerned  here  with  the  distribution  of  store  content,  Z(t).  A 
limiting  distribution  for  Z(t)  exists  which  is  independent  of  the 
initial  content  distribution  if 

}}2  TO*)  £ x | Z(0)}  = P{Z*  < x} 

for  some  random  variable  Z*.  If  instead  the  limiting  probability  is  not 
completely  independent  of  time,  but  rather  exhibits  a cyclical  pattern, 
then  we  say  we  have  a quasi-limiting  distribution,  i.e.  we  have  a quasi- 
limiting  distribution  if 

P{Z(k5+T)  < x)  | Z(0)}  = P{Z*(t ) < x) 

for  some  u>  and  random  variables  Z*(t),  0 < <2. 

2.1  Some  Results  in  the  Storage  Theory  Literature 

Limiting  distributions  were  first  shown  to  exist  for  the  queue  with 
homogeneous  Poisson  input  and  negative  exponential  service  time  distri- 
bution. The  general  method  employed  focused  on  the  states  {E j } which 
denote  the  queue  size.  When  service  is  negative  exponential  these  states 
form  a Markov  process.  When  service  time  is  not  negative  exponential 
this  is  not  the  case.  But  TakScs,  in  1955,  showed  that  this  was  not 
necessary  to  prove  the  existence  of  a limit,  since  the  waiting  time  is 
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a Markov  process  if  arrivals  from  a Poisson  process.  Even  if  inter- 
arrival times  are  merely  from  a renewal  process,  the  arrival  times  are 
regeneration  points,  and  this  is  sufficient  for  the  proof.  Thus  Takacs 
was  able  to  prove  results  such  as  the  following  theorem. 

Suppose  A(t)  is  the  intensity  of  a nonhomogeneous  Poisson  process, 
the  random  variable  is  the  service  time  of  the  n*^  customer,  and 
F(x,t)  is  the  distribution  of  waiting  time. 

Theorem  3 (Takacs):  Let  the  expected  value  of  the  randan  variable  \n 
be  a = / xdH(x)=  / ( l-H(x)  )dx  and  let  A(t)=A  be  a positive  constant. 
If  Aa^l  then  the  limiting  distribution  F(t,x)=F*(x)  exists,  is  inde- 
pendent of  the  initial  distribution  Fq(x)  and  is  uniquely  determined  by 
the  equations 

F*(0)=l-Aa 

and 

» \[  F*(x)  - /°H(x-y)dF*(y)] 
dF* 

where  - denotes  the  right  hand  derivative.  If  Aa>l  then  the  limiting 
distribution  F*(x)  does  not  exist,  however  F(t,x)=0  for  each  x. 

Of  course  this  result  can  also  be  applied  to  standard  storage  and 
dam  theory.  No  new  limiting  results  appeared  concerning  these  processes 
until  interest  in  storage  processes  with  general  release  rules  was 
aroused  in  the  late  1960's  and  early  1970's. 

A release  rule  (or  rate)  r(u)  is  a function  such  that  in  any  time 
interval  (t,t+dt)  the  amount  of  content  in  a store  which  is  released  is 
r(Z(t))dt+o(t);  i.e.  the  amount  of  content  that  will  be  released  from 
the  store  per  unit  time  is  a function  of  the  content  of  the  store  at 
that  time.  In  queueing  theory  this  is  equivalent  to  saying 
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that  the  amount  of  work  load  processed  by  the  server  per  unit  time  is 
dependent  on  the  waiting  time.  The  standard  release  rule  is  r(u)=c,  u>0 
and  r(0)=0.  This  release  rule  is  essentially  independent  of  the  content, 
as  the  store  content  is  released  at  a constant  rate  as  long  as  the  store 
is  not  empty.  Some  general  release  rules  which  have  been  considered 
which  are  dependent  on  store  content  are  r(x)=a+bx  and  r(x)=cxa. 

Storage  processes  with  general  release  rules  which  depend  on  the  con- 
tent of  the  system  have  more  recently  become  a subject  of  interest.  In 
1971  Cinlar  and  Pinsky  opened  up  the  field  with  a landmark  work  in  this 
area.  They  focused  on  a storage  process  with  input  which  has  stationary 
and  independent  increments  and  finite  jump  rate,  and  a release  rule  r(x) 
which  is  a Lipschitz  continuous,  strictly  increasing  function.  They  were 
able  to  show  that  if  the  mean  input  is  less  than  the  supremum  of  r(x) 
over  all  x,  then  the  storage  content  has  a limiting  distribution  which 
is  independent  of  initial  conditions. 

Brockwell  (1977)  extended  this  work  to  show  that  under  certain  condi- 
tions, a store  with  input  which  is  a pure  jump  Levy  process  (with  perhaps 
infinite  jump  rate)  will  admit  a stationary  distribution  of  store  content. 
In  particular,  in  Lemma  5.4  he  proved  that  for  r(x)  strictly  positive, 
nondecreasing  and  continuous  on  ( 0 ,°°) , if  r(0+)>0,  expected  input  in 
(O.t)  finite,  and  s^p  r(x)>EAj  then  a stationary  distribution 
exists  and  is  equal  to  the  limiting  distribution. 

In  1976  Harrison  and  Resnick  derived  necessary  and  sufficient  condi- 
tions for  the  existence  of  a stationary  and  limiting  distribution  of 
store  content  in  the  case  where  the  Input  process  is  Compound  Poisson 
with  finite  jump  rate,  and  the  release  rule  meets  the  following  condi- 
tions: 
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i)  r(x)  is  strictly  positive  for  x > 0 

ii)  r(x)  is  left  continuous 

iii)  r has  a strictly  positive  right  limit  everywhere  in 

0 < x < 


>v)  «*>-  4 Z-  < 


With  these  assumptions  not  only  were  Harrison  and  Resnick  able  to  find 
necessary  and  sufficient  conditions  for  the  existence  of  a stationary 
(and  limiting)  distribution  of  store  content,  but  also  a representation 
of  the  distribution.  They  proved  the  following  two  theorems: 

Suppose  X is  the  input  rate  and  F(  ) the  distribution  of  input 
size.  Define 

Q(x)=  X ( l-F(x) ) 

K(x,y)=Q(x-y)/r(x) 

Kn+i(x.y)%X  K(x,z)  Kn(z,y)dz  0 < y v " 


K*(x,y)=n'r1  Kn(x,y) 


Theorem  1 (Harrison  and  Resnick):  The  contents  process  has  a station- 
ary distribution  iff 


1 


- = 1 + / K*  (x,o)  dx  < >* 


in  which  case  the  unique  stationary  distribution  x has  density 


g(x)  = Yq  K*(x,0)  on  (0,"’)  and  > ( 0 ) = >0 - 


Theorem  2 (Harrison  and  Resnick):  is  positive  recurrent  iff 

K*(x,0)dx  < <*>  iff  X^  has  a stationary  distribution  x . In  this  case 
the  limiting  distribution  II  coincides  with  the  unique  stationary  distri- 
bution in  Theorem  1. 


..  .r—  -r-r.-rr 
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While  the  input  distribution  considered  by  Harrison  and  Resnick  is 
more  restrictive  than  before,  the  release  rule  considered  here  is  cer- 
tainly more  general  than  that  of  either  Cinlar  and  Pinsky  or  Brockwell. 
Unfortunately,  as  they  pointed  out,  assumption  (iv)  rules  out  any 
release  rule  which  makes  it  impossible  for  the  store  to  empty. 

As  mentioned  in  Chapter  1,  the  problem  of  proving  the  existence  of 

a limit  when  the  input  process  is  not  homogeneous  in  time  has  only  been 

♦ 

examined  for  the  M(t)/G/1  queue.  There  we  have  two  results.  The  first, 
stated  earlier  in  this  chapter,  is  for  an  intensity  which  approaches  a 
finite  limit.  The  second,  proved  separately  by  Sakr  (1972),  Yevdokimova 
(1974)  and  Harrison  and  Lemoine  (1977),  shows  that  for  a queue  with  non- 
homogeneous  Poisson  arrivals  with  periodic  intensity,  a limiting  distri- 
bution of  waiting  time  will  exist  if  the  average  traffic  intensity  per 
period  is  less  than  one. 

2.2  Some  Basic  Tools 

We  extend  these  results  here  to  storage  processes  under  a variety 
of  conditions  on  the  input  process  and  the  release  rule.  For  the  most 
part,  the  input  processes  we  consider  will  be  much  more  general  than 
those  that  have  been  considered  in  the  past.  One  restriction  we  add  is 
the  assumption  that  each  realization  of  the  storage  process  can  be  repre- 
sented by  the  infinite  sequence  of  inputs 

uj  * (Xj,  ur  x2,  u2,  ...  1 

where  X^  is  the  time  between  the  i-lth  and  ith  input  and  U.  is  the  size 
of  the  ith  input.  We  call  such  input  discrete  input. 

*We  adapt  the  standard  queueing  theory  notation  to  let  M(t)  denote 
a nonhomogeneous  Poisson  arrival  rate. 
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A method  of  proof  that  will  be  used  often  throughout  this  chapter  is 
to  compare  two  storage  processes  with  the  same  input  w,  i.e.,  we  will 
compare  Zj(t)«Z(t,ui,f;,r,)  and  Z2(t)=Z(t,w,^2,r2)  where  Z(t ,«>,£, r)  denotes 
the  storage  content  at  time  t of  a store  with  release  rule  r,  initial 
content  f.»  and  input  w.  In  forthcoming  proofs  various  assumptions  will 
be  made  on  f , and  r,. . But  by  assuming  that  for  a given  realization  both 
processes  receive  inputs  at  precisely  the  same  moment  and  have  the  same 
size,  we  remove  any  differences  which  could  be  caused  by  the  input 
structure.  Thus  we  can  more  readily  compare  the  contents  of  the  two 
processes. 

We  will  consider  here  two  types  of  storage  processes.  A Type  A pro- 
cess is  defined  as  a storage  process  with  release  rule  r(x)  such  that 
for  some  "a",  r(a)=0  and 

(1)  R(a,x)  = / rjjy  < * 0 < x < « 

<3 

The  function  R(a,x)  represents  the  amount  of  time  it  would  take  the  con- 
tent of  the  store  to  go  from  level  X to  level  a if  there  were  no  input. 
Thus  (1)  implies  that  for  a Type  A process  there  exists  a level  a which 
can  be  reached  from  any  level  X above  a (with  positive  probability)  in 
a finite  amount  of  time. 

A Type  B process  is  a storage  process  with  release  rule  r(x)  such 
that  for  some  level  a,  r(a)=0  and  for  every  r N 0 

(2)  R(a+c,x)  < “ a + e<x<«' 

but  there  is  an  X such  that 

(3)  R(a,x)  = - 

Li 
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Thus  for  a Type  B process  there  is  a level  a which  cannot  be  reached 
from  above  in  a finite  amount  of  time,  even  if  there  is  no  input,  but 
the  level  of  the  store  content  can  get  arbitrarily  close  to  a. 

Note  that  without  loss  of  generality  we  can  assume  a=0.  When  a>0 
we  simply  have  the  situation  of  a store  with  a "false  bottom".  For  sake 
of  simplicity,  in  the  future  we  will  assume  a=0  unless  otherwise 


specified. 

At  this  point  it  is  convenient  to  derive  a few  properties  of  Type  A 
and  B storage  processes  with  general  release  rules.  We  define  a con- 
sumption curve  as  the  graph  a realization  of  store  content  Z(t)  would 
follow  if  there  were  no  input.  These  curves  are  described  by  the  differ- 
ential equation 


(4) 


dZ(t)  = -r  (Z(t)) 
dt 


We  take  Z(t)  to  be  right  continuous,  i.e.  Z(t+0)  = Z(t).  Thus  (4)  implies 
that  Z(t-)  satisfies 


(5) 


0 

du 


rTUy 

Z(t-) 


= t - t 


From  (5)  we  can  see  that  Z(t)  is  continuous  except  where  inputs  occur, 
since  if  there  is  no  input  in  (tQ-  tQ) 


lim 

n-*<x> 


rz'lo-K> 

du  _ 

r(u)  n-**> 


[‘o-'Vs’]'0 


i.e.  Z(tQ-0)  • Z(tQ). 
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Thus  we  have  that  the  consumption  curves  are  continuous  and  strictly 
decreasing  functions  of  t passing  through  the  point  (tQ,  Z(tQ))and  every 
positive  level  less  than  Z(tQ).  The  latter  is  true  since  for  a Type  A 
or  B process 

b dx 

/ FtTJ-  ' ” 0 < a < b 

d 

We  prove  now  some  Lemmas  concerning  the  storage  processes. 

Lemma  2.1.1:  Let  Z^(t)  = Z(t,o>,  £j,r)  and  Z^(t)  = Z(t,u>,  £2>r)  the 

content  of  two  storage  processes.  If  for  some  t Z0(t  ) >Z.(t  ) then 

o c o — 1 o 

Zg(t)  >_  Zj(t)  for  all  t>t 

Proof:  Let  tj  be  the  time  of  the  first  input  after  t . We  first  show 
that  Z2(t)>_  Z^(t)  for  tQ  < t <ty  In  the  case  Zj(tQ)  £ Z2(tj)  thls  1S 
clear,  since  Zj(tj)  1 Z^(tQ)  1 Z2(tj).  Thus  we  need  to  look  at  the 
case  Z2(tj)  < Zj(t  ).  Suppose  the  conclusion  we  seek  does  not  hold  and 
Z2(t)<Zj(t)  tQ  <_  t<t^.  Then 

fVV  o> 

*‘*0  = J fTST  ' ) + 

JZ2(t)  Jzl(t0) 


But  we  also  have 


cTu) 

Zj(t) 


This  is  impossible  since  r(u)X)  and  ^2( tQ) >2j ( tQ)  and  Zj( t) >Z2( t) . Thus 
we  can  conclude  that  Z2(t)>Zj(t)  for  tQ  < t<tj. 

Suppose  at  tj  there  is  an  input  of  size  U.  We  know 

z2(tr)  > Zj(tr) 

Then 

z2(ti)  3 z2(tr)  + U - Zl(tl_)  + U = Zl(tlJ 

Thus  we  know  that  Z2(t)  >Z:(t)  both  between  inputs  and  directly  after 
inputs,  which  is  at  all  times;  i.e.  we  can  conclude 

Z2(t)  >.  Zj(t)  for  all  t 


Lenwa  2.1.2;  Suppose  r(u)  is  nondecreasing.  Then 


pb+c 

du  < / du 
r(u)  rTuT 
a+c  Ja 


I. 


a ,b,c>0 


b+c 

Proof : For  the  integral  { du  we  make  the  substitution  y=u-c. 

rTuT 

a+c 


This 


gives 


J-b+c  rt 

T du  = | dy 

r(u)  I r( y+c) 

a+c  -'a 


then 


^ - [V  ( - ( & 

c ■'a  a a 

C sfeHist}  * 


■ — — — - — - ■ — - 
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Since  r(y)  is  nondecreasing  r(y)  < r(y+c).  Thus  we  have 


/-b+c 
du 

FTuT- 

a+c 


Lemma  2.1.3:  Let  Z ^ ( t ) = Z(t,w,t.j  ,rj ) and  Z^(t)  = Z( t,u),.f.,r.,)  be  the 
contents  of  two  Type  A or  B storage  processes,  with  r^(x)vr^(x).  If 
for  some  tQ,  Z ^ ( tQ ) ^Z2( ) then  Zj(t)>Z0(t)  for  all  t^t  . 

Proof:  Let  t,  be  the  time  of  the  first  input  after  t . For  t ^t^t, 
1 r o o 1 


t - t. 


rW 

du 


r’jTu) 


Zj(t) 


Z,  (t  ) 

r 1 0 

du 

r2fu i) 
'Zj(t) 


since  rj(x)<r2(x).  The  assumption  that  Z](t0)"Z2(to)  iniP^ies  that 


-W 


n 

Zj(t) 


du 

r2Tu) 


rW 

du 


r^u) 


Zj(t) 


Suppose  that  Z^(t)<Z2(t).  Then 


rW 

I du 

rCuT 

Zj(t) 


W 

du 

rTuT 


rZ2(tc 

J rj 

Jz2(t) 


= t-t 


a contradiction.  Thus  for  t^t^tj,  Z j ( t ) > Z 2 ( t ) . 

Suppose  at  t^  there  is  an  input  of  size  U.  Then 


z^tj)  * u ♦ Zj( tj-)  > z2(tr)  ♦ 1)  - z2(t.j). 
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Since  we  have  Z1(t)>Z2(t)  in  between  times  of  Inputs  and  directly  after 
an  Input,  it  Is  clear  that  Z1(t)>Z2(t)  for  all  t>tQ. 

2.3  A General  Limit  Theorem 

We  now  proceed  to  results  concerning  the  limiting  value  of  store 
content.  The  first  theorem  shows  that  under  certain  circumstances  the 
existence  of  a stationary  distribution  is  sufficient  to  insure  that  a 
limiting  distribution  will  exist  which  is  independent  of  the  initial 
distribution. 

Theorem  2.1:  Let  Z(t)  be  the  content  of  a Type  A or  B storage  process 
with  nondecreasing  release  rule  r(x)  and  discrete  input.  Suppose  further 
that  for  any  e>0  and  initial  content  t there  exists  a proper  random 
variable  T(£,e)=T  such  that  Z(T)=e.  (If £<e  then  define  T=0).  If  the 
process  has  a stationary  distribution  then  it  has  a limiting  distribu- 
tion which  is  independent  of  the  initial  distribution,  and  the  two 
distributions  are  equal. 

Proof:  Suppose  we  have  two  storage  processes  of  Type  A or  B,  Z^(t)= 
Z(t,u>,£j,rj)  and  Z2(t)=Z(t,a>,C2,r2)  where  r^(x)=r^(x)  is  nondecreasing 
and  both  processes  satisfy  the  assumptions  of  the  theorem.  Let 
By  Lemma  2.1.1 

(6)  Z1(t)>Z2(t)  for  all  t 

It  follows  that  T(  Si  ,e)  >_T(£;2,e),  which  together  with  (6)  implies  that 

(7)  0 < Zj  (T(cre))  - Z2(T(q,e))  < e 


Suppose  that  the  first  input  after  T(f, j,e)  occurs  at  tj.  If  we  can  show 

(7)  holds  for  all  T(f, ^,e)<t<tj  then  clearly  it  will  hold  for  all  t>T(Fj,e). 
Let  T(C j,c)<t<t1>  Since  there  have  been  no  inputs  in  this  interval 

we  still  have 

(8)  e = Zi  (Ttej.e))  > Z^t)  1*1,2 


0<Z J^(t)  — Z2(t)  < e for  T(£>1  ,e)<t<tj 
Now  suppose  at  tj  we  have  an  input  of  size  U.  Then 

Zi(tl)=Zi(tl")+U  i=1 »2 

and  we  have  from  (6)  and  (8)  that 

U < Z2(ti)  < Zj(tj)  < U + e 
Hence  we  can  conclude  that  for  all  t>T(£j,c) 

0 < Zj(t)  - z2(t)  < e 

Similarly,  if  f, we  ^ave  that  for  t>T  (^.c) 

0 < Z2(t)  - Zj(t)  < e 

From  which  we  can  conclude  that 

|Z1(t)-Z2(t)  |<  e for  t>max  [t^.e),  T(C2.e)]  = T* 

We  can  now  easily  show  that  Zj(t)—  Z2 ( t ) -*0  in  probability  for 
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PdZjtt)— Z2(t)|  >ek1.C2}  - 

p{ | Zj(t) — Z2(t)|  >e  and  t > + 

P{  |Z1(t)— Z2(t)  | >e  and  t < T*\ZVZZU 
0 + P{t<T*} 

But  T(C^»e)  is  proper  for  each  i,  and  thus  so  is  T*.  which  means 
lim  P{|Z1(t) — Z2(t)|  >e  | C1»C2>  = 0 

Let  Fc  (x,t)  = P{Zi(t)  < x | Z . (o)  = £.}  1-1,2. 

Since  Zj(0)  is  independent  of  Z2(0)  we  have 

(x,t)-P(Z1(t)  < x |Clt52} 

= P{Zj(t)  < x and  | Z1(t)-Z2(t)  | < eKj.Sg) 

♦ P(Zj(t)  < x and  | Z^t)— Z2(t)  | >e  ISj.Cg) 

< P{Z2(t)  < x +e  + P{ IZjCt) — Z2(t)  | > c *^2) 

i .e. 

(9)  F (x.t)  < F?  (x+e,t)  + P{|Zj(t)-Z2(t)|  >e  J 

Suppose  now  that  has  distribution  G^.  From  (9)  we  have,  taking 
expectations 

Fj  (x.t)  < F2  (x+e ,t)  + E [P{|Zj(t)-Z2(t)|  >e  | q,^}] 

where  F.  (x.t)  s E Ff  (x.t)  1=1,2. 

1 M 


We  note  that  by  dominated  convergence 


lim  E [ P{  |Zj(t) — Z2(t)  | > e | } ]=  0 

Thus 

(10)  lim  F,(x,t)  £ lim  F„(x+e,t) 

Similarly,  by  interchanging  the  roles  of  and  and  substituting  x-e 
for  x we  can  get 

ljm  F?(x-e,t)  < 1 im  F.(x,t) 
t-X”  t"X»  ^ 

By  assumption  a stationary  distribution  H exists  for  each  storage 
process.  Thus  we  may  let  G2(x)  = H(x).  Then  H(x)  = F2(x,t)  for  all 
t,  and  from  (10) 

1 im  F. (x,t)  < H(x+c) 

t-xo  1 

Letting  c go  to  zero  we  can  get 

lim  F^x.t)  £ H(x+0) 

Simi larly 

H(x-O)  £ Vim  F. (x,t) 
t«-  1 

Thus 

H(x-O)  £ lim  F.(x,t)  £ lim  F.(x,t)  £ H(x+0) 

" tx»  1 t 1 

At  each  continuity  point  of  H we  have 

lim  F.(x,t)  = H(x) 

t-+°° 

the  desired  conclusion. 


2.4  Limit  Theorems  for  Type  A Processes  with  General  Periodic  Input 
Theorem  2.1  requires  a very  strong  assumption,  the  existence  of  a 
stationary  distribution,  which  in  itself  is  often  difficult  to  verify. 

The  next  theorem,  and  those  which  follow,  will  not  require  such  a strong 
assumption.  In  exchange,  more  structure  will  be  put  on  the  input 
process. 

We  suppose  there  is  a number  co  (called  the  period)  and  in  each  inter- 
val [kw,  (k+l)£>)  there  are  a finite  number  of  times  where  there  is  a 
positive  probability  of  an  input  to  the  system.  Let  these  points  be 
denoted . 

koi+Oj  <kG+G2<  •••  <kui+Op 

where  the  are  the  same  for  each  interval.  Each  point  has  associated 
with  it  a distribution  of  input  size  G.(x),  where  we  allow  G.(0+)>0.  The 

J J 

Laplace-Stieltjes  Transform  of  these  inputs  in  (0,t]  t=kw+T,  0<t<w,  is 

r „ ik  , 


' P * 1 1 * 

M.(s,t)  = n g.(s)  n G.(s) 

U=i  J j=i  3 


where 


01  < t < 0n 


* 

n G.(s)  is  defined  to  be  1 if  t < 0. . 
j=l  3 1 


We  also  suppose  that  in  any  interval  (t,t  + dt)  there  is  a probability 
X (t)  dt  + o(dt)  that  an  input  will  occur,  where  X (t)  is  taken  to  be  a 
periodic  function,  i.e.  X(kai+T)  * X(t),  0<T<oi.  The  size  distribution  of 
these  inputs  B (x,t)  is  also  taken  to  be  periodic.  The  Laplace-Stieltjes 


27 


Transform  of  these  inputs  into  the  system  in  (0,t)  is 

M^s.t)  = exp  {-  /g  (1-B  (s,u))  A(u)du) 

The  total  input  to  the  storage  system  is  a convolution  of  these  two 
inputs,  and  has  laplace-Stieltjes  Transform 
M(s,t)  = M^s.t)  M2(s,t) 

We  will  call  any  input  process  which  is  characterized  by  this  Laplace- 
Stieltjes  Transform,  general  periodic  input  or  g.p.i.  for  short.  The 
i^  moment  of  input  size  will  be  denoted  by 

The  g.p.i.  is  really  much  more  general  than  we  have  encountered  in 
the  literature.  Some  storage  theory  assumes  discrete-time  input,  some 
have  continuous  time  input,  but  none  have  both  in  combination  as  we  have 
here. 

Another  important  feature  of  the  g.p.i.  is  that  it  is  not  homogeneous 
in  time.  In  particular,  both  the  input  size  distribution  and  the  rate  of 
occurrence  of  inputs  is  periodic.  Most  work  in  this  area  has  required 
that  inputs  be  independent  of  time,  with  stationary  increments.  This 
implies  that  the  distribution  of  input  in  one  time  interval  is  equal  to 
the  distribution  of  input  in  any  other  time  interval  of  the  same  length, 
as  mentioned  earlier.  With  periodic  input,  we  have  that  the  distribution 
of  input  in  one  time  interval  which  begins  at  some  point  r from  the 
start  of  a period,  0 < t < wis  the  same  for  all  periods.  This  fact  allows 
regenerative  events  to  be  of  use  here,  as  they  have  been  of  use  in  the 
past  for  the  homogeneous  case. 

It  is  clear  that  a change  of  time  scale  can  transform  the  g.p.i.  to 
have  period  u * 1.  Thus,  without  loss  of  generality  we  will  assume 
throughout  the  remainder  of  this  work  that  m = 1. 


! 


I 


.1 
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Theorem  2.2:  Let  Z ( t ) be  a Type  A storage  system  with  general  periodic 
Input  and  nondecreasing  release  rule  r(x).  Define 

¥<u>  ■ m 

If 

P CO  1 ■» 

u = / Y(u)G,(du)  +/' V Y(x)B(dx,t)X(t)dt<l 

y o 1 o o 

then  the  store  content  has  a quasi-limiting  distribution  which  is  inde- 
pendent of  the  initial  distribution. 

Proof:  Let  Z(t)  be  the  store  content  at  time  t and  define  the  event  t- 
0 £ t < 1 as  follows:  eT  occurs  at  t=n+l  (n  an  integer)  if  Z(n+T)  = 0. 
Clearly,  is  a regenerative  event,  as  defined  in  Smith  (1955).  Note 
that  eT  is  periodic  with  period  1.  For  a quasi -1 imiting  distribution  to 
exist  for  Z(n+£)»  0 < S < 1 we  must  show,  for  a useful  class  of  sets  A, 
that 

Al)  c is  certain,  i.e.  c will  occur  a first  time  with  probability 

T T 

one. 

A2)  n"Q  4>(n+s)  (l-F(n+s))  < » 


where  *A(t-m)  = P(Z(t)rA  |Z(0),nt>0,  Tn  = ml 

nt  = number  of  occurrences  of  in  [0,t] 

T„  a time  of  last  occurrence  of  e before  t 
nt  T 

To  prove  Al  and  A2  we  will  use  a series  of  Lemmas.  Before  looking  at 
the  Leninas,  though,  we  need  to  introduce  the  concepts  of  empty  time  and 
busy  time.  We  define  the  empty  time  of  the  store  in  the  interval  [0,t), 
denoted  e(t),  as  the  Lebesgue  measure  of  the  set  (s  < t:  Z(s)  * 0). 
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In  a corresponding  manner  we  define  the  busy  time  in  the  interval  [0,t), 
denoted  b(t),  as  the  Lebesgue  measure  of  the  set  {s  < t:  Z(s)  > 0}. 

Notice  that  e(t)  + b(t)  s t.  In  addition,  we  can  define  a busy  cycle. 

A busy  cycle  is  defined  to  begin  when  the  store  empties,  and  ends  when 
the  store  next  empties.  The  initial  busy  cycle  begins  at  t=0  and  ends 
when  the  store  first  empties.  We  now  proceed  to  the  Lemmas. 

Lemma  2.2. 1:  Suppose  Z(t)  is  the  content  of  a Type  A storage  process  as 
described  in  Theorem  2.2.  Then 

lim  E [e(n)]  > 0 

tv-**  n 

Proof:  To  prove  this  lemma  we  need  to  construct  an  alternative  storage 
model  which  we  shall  call  a warehouse  model:  The  warehouse  system  will 
receive  the  same  input  u as  the  original  system.  If  an  input  arrives 
while  the  alternative  store  is  not  empty  then  the  input  is  stored  in  a 
separate  warehouse.  When  the  alternative  store  is  about  to  empty  the 
next  input  waiting  in  the  warehouse  is  immediately  put  into  the  store  to 
prevent  emptiness  unless  an  input  arrives  at  that  moment,  in  which  case 

A 

this  input  is  added  to  the  store,  preventing  emptiness.  Let  Z (t)  be 
the  store  content  in  the  alternative  warehouse  model,  and  Zw(t)  be  the 
warehouse  content. 

We  claim  that  the  total  content  in  the  warehouse  system  (both  store 

and  warehouse)  is  at  least  as  great  as  the  content  in  the  original  model, 
A W 

i.e.  Z (t)  + Z (t)  Z(t).  We  prove  this  by  induction  on  the  number  of 

inputs.  Suppose  inputs  occur  at  times  t^,  t2*...  with  size  U(tj),  U^), 
...  For  0 < t < tj  ZW(t)=0  and  ZA(t)=Z(t). 


W. 


Define  T0  * Inf  {t:  ZA(t)  ■ 0}. 

Suppose  tx  < TQ.  At  2W( t1)  * U( tj)  and  ZA(tr)  = ZA(t1).  Also 

Zdj)  * Z( tj-)  + U( tj) . But  we  know  ZU^)  « ZA(t]-).  Thus 

Z(tj)  « Z(tj-)  ♦ U(tj) 

- ZA(tj)  ♦ ZW(t1). 

A 

Suppose,  on  the  other  hand,  that  To  < tj.  Then  Z (tj-)  = 0 =*  Z(tj-) 
and  ZW( tj ) » 0.  Also 

ZA(tj)  - U(tj)  = Z(tj)  and  ZW(t1)  = 0, 

since  the  arriving  input  is  put  directly  into  the  store,  giving 

Z(tj)  * U(tj)  - ZA(t1)  + zw(t1). 

Thus 

Z(t)  < ZA(t)  + ZW(t)  for  t < tj 

Now  suppose  we  know  that  ZA(t)  + ZW(t)  > Z(t)  for  all  t ^ t^.  We 
want  to  prove  that 

ZA(t)  + ZW(t)  > Z(t)  for  t < tk+1 

For  definiteness,  let  U( 1 1 ) , U(t^  + 1^,...  U(t^)  be  the  Inputs  stored  in 
the  warehouse  and  define 

Tk  1 ” Tk  ” inf  **  - V 
Tk,l+1  * inf  {t  > Tk  ZA(t)  - 0) 


There  are  a number  of  cases  to  consider. 


C5i<a>  tk  < t < tk+l  < Tk 


Suppose  tk  < t<  tk+J 


Recall  that 


^k  - 


r^V 

l du  = 

I rfuT 

Z(t-) 


From  Lemma  2.1.2  we  have 


zA(t. ) 


ZA(t-j 


k' 

du 

rfuT 


/'V  /(tkW"(tk) 

du  — / du 

J a ^ uT  J A rTuL 

ZA(t-)  JZA(t-)+r(t.) 


r ii  r 

J,a,/F  l 


and  we  know  by  assumption  that  Z(tk)  < Z^( tk ) + Z1^ ( tk) 


■zA(tk>  y“<V 


r(t-)  +-¥ "(tk) 


z(tk) 


du 

A rf'u)  JJ 

ZA(t-)  + zw(tk) 


which  yields 


z(tk) 


du  - 


f(u7 

Z(t-) 


du 

A r'ft‘T  u 

ZA(t-)  + 2“(tk) 


For  the  inequality  to  hold  we  must  have  Z(t-)  < ZA(t-)  + ^(tk). 
Since  t / Tk,  ZA(t-)  = ZA(t). 

Also,  for  t c (tR,  tk+1)  Z(t-)  = Z(t)  and  ^ (\)  = Z^t). 

Thus  we  have 


Z(t)  < ZA(t)  + Z^t) 


\ < 1 < lk+l 


At  t = t, 


Z^tk+P  Z^k+1^  + U^k+P 

i ZA(tk-,)  ♦ z“(tk-j)  ♦ u(tk+1) 

- zA(tktl)  ♦ z"(V]) 


The  latter  is  true,  since  if  tk+1  < Tk  then 


z (tk+1)  * zA(tk-x) 


z“(t,«)  ■ * Ului 


or  if  tk+1  = Tk 


z (tk+i>  ■ uttk+1) 

a-d  ZW(tk7l)  - Z“(tk+1) 

since  if  an  input  arrives  just  as  the  alternative  store  empties,  this 
input  goes  into  the  store.  Thus  we  have 
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z(t)  < zA(t)  * z“(t) 


for  tk  < t < tk+,  < Tk 


lk  < ‘ i Tk  < Vl 


Suppose  t < Tk>  As  we  did  for  Case  I we  can  get 


fZ(tk’  , 

! du 

J rTu T 

Z(t-) 


,ZA(tk) 


r k'  r 

du  - 

A r(u) 


'r(t-) 


zA(tk)+z“(tk) 

du  : 

A 

zA(t-)+zw(tk) 


.Z(tk) 

du 

A r^u)  u 


'z^(t-)+zri(tk) 


Thus 


Z(TJ  = Z(T. -)  < ZA(T-)  + ZW(Tr)  = ZA(TJ  + ZW(Tj 


Case_m:  < TkJ<  ...  <Tkfj  < t < TkjJ.+1  <•••  <Tkfl  <tk+l<TkJ+l 


There  can  be  one  or  more  times  before  tk+1  where  the  store  in  the  ware- 


house system  will  empty.  We  look  at  the  interval  Tk  i < t £ Tk  2 which 
is  typical  of  the  others.  Again,  for  t < Tk  ^ 


r2<Tk,l> 

rzA(Tk,i}  . 

rzA<Tk.l>+zW<Tk,i>  r 

du  = 

du  > 

du  > 

FTDT 

J A r(u) 

A r'M  u 

ZA(t-)+ZW(Tk>1) 

Z(t-) 

^ZA(t-)  J 

n) 
.A 


du 

rTul 


W# 


k,r 


Thus  we  have 


Z(t-)  < ZA ( t- ) + Zw(Tk>1), 





— - v 


K 
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But  Z , t ) = Z(t-),  Z (Tk^)  = Z (t)  and  Z (t-)=Z  (t),  since  no  Inputs  or 
emptiness  occur  at  t.  Thus 

Z(t)  < ZA(t)  + ZW(t)  for  Tk>1  < t < Tk  2 

For  t = Tk,2  we  have  zA(Tk>2)  + zW(Tk,2)  = zA<Tk72>  + ZW(Tk~2) 
which  means 

Z(Tk,2>  - Z<Tk72>  i zA'Tk72»  * Z“<T«72>  • ZX.2>  + Z“^k,2) 

So  for  Tk>1  < t < Tfc+1 

Z(t)  < 2A(t)  ♦ z“(t). 

Clearly,  it  will  also  be  the  case  that 

Z(t)<  At)+Z“(t)  for  TkJ  < t < Tk>jtl 

We  now  look  at  the  last  case. 


tm  < 'iVliVm 


Suppose  Tk  1 < i 

t < tk+j.  As  usual 

rz(Tk,i* 

Am’ 

zfl(Tkil)fz"(Tki,) 

rz(Tk,1> 

i du  = 

I du  > 

du  > 

I du 

J rTuT 
Jz(t-) 

JZA(t?P  • 

1 A 

;ZA(t-)+ZW(TkJ) 

J A r('u)  u 

JZA(t-)+ZW(T| 

which  implies 

Z(t-)  < ZA( t- ) + ZW(Tkjl). 

For  1 < lk+i  we  have  z(t_)  = Z(t^  and  zW(Tkfl)  = z(t). 
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Also,  since  t e (Tk  j , Tk  j + j)  we  have  ZA(t-)  * ZA(t).  This  gives 
Z(t)  < ZA(t)  + ZW(t). 

For  t « 

Z(W  ■ Z<W  * u<tWl)  i zA(tk¥l>  * 2“ttk*l)  * »(Vl> 

But  we  know  that  either  ZA(tk?1)  = ZA( tk+1)  and  ZW(tk+1)  = Z^t^)  + 
U(tk+1)  if  tk+1  < Tkj+1.  or  if  = Tk,l+1  then  1 ^k+l^  = 
zW( tk+P  = z^( tk+i ) , zA(fck+i^  = U^k+1^  In  either  case 

Z Uk+1)  + z (^+1)  + u<tk+i)  “ 1 ^k+l^  + zW^lk+i^ 

Thus 

Z(t)  < ZA<t)  ♦ z“(t>  for  Tkil  <t<tktl<TkU1. 

With  these  four  cases  wo  have  shown  that 

Z(t)  < ZA  (t)  + ZW( t)  for  t < tk+1 

and  by  induction  this  then  holds  for  all  t. 

Recall  that  we  denote  the  busy  time  and  empty  time  by  b(t)  and  e(t) 
respectively.  For  the  alternative  warehouse  model  we  will  use  bA(t)  and 
eA(t).  Since  Z(t)  < ZA(t)  + Z^(t)  we  have  e(t)  > oA(t)  and  b(t)  < bA(t). 

U(tj) 

We  call  Y(U(t.))  = f du  the  standard  busy  time.  We  can  think 

J o rTuT 

of  Y as  the  amount  of  busy  time  associated  with  the  input  U(t.)  in  the 

V 

warehouse  model.  We  note  that  the  busy  time  associated  with  U(t.)  in  the 
original  model  will  be  no  bigger  than  this  standard  busy  time. 
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For  the  warehouse  model  we  have  the  inequality 

(11)  E [b\n)]  < E [y  (Z(0))]  +nyy 

The  first  term  on  the  right  hand  side  of  (11)  represents  the  stand- 
ard busy  time  associated  with  the  initial  content  Z (0 ) - The  second 
term  represents  the  standard  busy  time  from  all  inputs  to  arrive  at  the 
warehouse  system  in  (0,n).  The  inequality  indicates  that  possibly  all 
inputs  arriving  in  (0,n)  are  not  used  up  during  that  time  (i.e.  possibly 
ZW(n)  > 0)  so  the  right  hand  side  of  (11)  may  be  an  overestimate. 

Since  eA(t)  = t - bA(t)  we  have 

E [e(n) ] > E [eA(n)] 

= n - E [bA( n ) ] 

> n - npy  - E [y(Z(0))] 

Then 

E [e(nj]  > 1 - u - E [y(Z(q))J 
n y n 

and  since  by  assumption  1 - py  > 0 

lim  E [e(n)]  > 1 - u > 0 
n y 

Lemma  2.2.2:  With  a storage  process  as  defined  in  Theorem  2.2,  there 
exists  t,  0 < t < 1 such  that  eT  will  recur  with  probability  1,  and  the 
mean  recurrence  time  pt  is  finite. 

Proof:  Let  C*  = J 1 if  Z(t)  = 0 

r 'So  if  Z(t)  > 0 
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Then  /n  r.  dt  = e(n) 
o 

and  I ( n)  = E [e(n)] 
n 

= / P (Z(t)  = 0}  dt 

= /"  "e1  P {Z( j+x)  = 0}  dr 
j=o 

From  Lenina  2.2.1  we  know  that  Vim  i I (n)  > 0. 

n~x» 

Then 

lim  / "e1  P {Z( j+x)  =0}  d x > 0 
n-K»  o j=o 

We  can  think  of  the  event  as  a delayed  recurrent  event  where 
stands  for  "a  success  occurs"  (system  is  empty  at  t=n+x)  in  a sequence 
of  Bernoulli  trials.  Pursuing  this  (in  the  manner  of  Feller)  we  note 
then  that 

n-1 

UT(n)  = P (Z( j+T)  = 0} 

is  the  expected  number  of  renewals  in  [o,n]„  What  can  we  say  about 

U (n)/n? 
x 

If  = oo  or  the  probability  recurs  is  less  than  1 then 

lim  U (n)  = 0. 

x 

n-  — 

If  recurs  with  probability  one  and  the  first  occurrence  of  has 
probability  p(t)  > 0 then 

lim  U (n)  = p(x) 

* 
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In  any  case,  it  is  always  true  that  Ur(n)/n  coverges  to  a limit  as  n 
goes  to  infinity.  We  will  call  this  limie  O(t).  Since  0(x)  < p*  we 
have,  by  dominated  convergence,  that 

1 

lim  I (n)  = ljjn  £ U (n)  dr 

n-*»  n n-*»  — — 

n 

fi 

= 0 0(t)  dx 

> 0 


Thus  we  must  have  that  0(x)  > 0 on  a set  of  positive  Lebesgue  measure. 
Let  0 = (x:  0(x)  > 0).  For  x e 0 we  have  that 


which  implies  that  uT  < °»  , cT  recurs  with  probability  one  and  p (x)  > 0. 


Lemma  2.2.3:  For  a storage  system  as  described  in  Theorem  2.2,  if  x e 0 

then  the  first  occurrence  of  c is  certain. 

x 

Proof,:  Suppose  we  have  two  storage  systems,  Zj(t)  = Z(t,w,Lj,r)  and 
Zg(t)  = Z(t,co,  ^2»r)  which  are  identical  except  for  initial  store  content, 
and  ^2  > = Lemma  2.1.1  we  know  that  Z^t)  <_  Z2(t)  for  all  t. 

If  there  is  a tQ  = tQ  (Cg.uj)  such  that  ^(t^  = 0 then  Zj(tQ)  = 0 and  for 
t >,  tQ  Zj(t)  = ZgU).  We  will  prove  that  there  must  exist  such  a t . 
Suppose  not.  Then  e2(t)  = 0 for  all  t.  But  this  implies  that 
^(t)  = E2  (e(t))  = 0,  a contradiction  as  long  as  u < 1 (Lenina  2.2.1). 
Hence  there  must  exist  a t . 


( 


- — - - 


— 
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Suppose  now  that  we  have  two  processes,  Z3(t)  = Zft.u^.Cj.r)  and 
Z^(t)  = Z(t,<dj,£4.r)  such  that  Z3(x)  = 0 and  Z^(t ) = £;  for  some  t e 0. 

Construct  two  new  processes  with  t as  the  new  origin:  Z^(t)  = Z3(t+x) 

★ ★ 

and  Z^(t)  = Z^(t+T ) . We  know  there  exists  a tQ(C,w1)  such  that 

* * ★ 

Z3(t)  = Z4(t)  for  t > By  Lemma  2.2.2  we  know  that  for  some 

★ 

k > tQ(£;,u>j),  Z3(k+x)  = 0,  since  for  t e 0,e  recurs  with  probability 
one.  But  Z3(k+x)  = Z3( k ) = Z4(k)  since  k > t*U,Wj),  and 

'if 

Z^(k)  = Z4(k+x),  i.e.  Z^( k+x ) = 0.  Thus  we  have  shown  that  if  we  have 
a storage  system  which  meets  the  requirements  of  the  Lemma,  and  x e 0, 

then  c is  certain. 

x 

Proof  of  Theorem  2.2  continued:  Choose  x e 0.  From  Lemma  2.2.3  we 

know  e is  certain.  To  show  that  f <j>.  (n+f;)  (1  - F(n+F,))<  °°  note 
I n-U  n 

that 

oo  00 

n=0  <}>A(n+£)  (1  - F(n+0)<  nI0  (l-F(n+0)<  Wj  < 00 
by  Lenina  2.2.2.  Thus  in  conclusion,  by  Theorem  3 of  Smith  (1955) 
lim  P {Z(n+0  e A | Z(0)}  = I f *.(n+0  (1  - F(n+0) 

n*«> 

The  next  result  follows  directly  from  Theorem  2.2,  but  now  the  release 
rule  is  not  constrained  to  be  nondecreasing.  Instead,  we  assume  that  the 
release  rule  is  always  bounded  away  from  zero. 

Theorem  2.3:  Let  Z(t)  be  the  content  of  a Type  A storage  system  with 
general  periodic  input  and  release  rule  r(x)  such  that 

p 00  - oo 

r(x)  > a > jTj  o x d + o o x * (u)  du  ^or  each 
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x >0.  Then  there  exists  a quasi-limiting  distribution  of  store  content 
which  is  independent  of  the  initial  content  distribution. 

Proof:  Define  two  Type  A storage  processes  with  g.p.i. 

Zj(t)  = Z(t,u),£,rj)  and  Z^(t)  = Z(t,io,C,r2)  where  r^(x)  is  as  described 
in  the  theorem  and  r2(x)  = a,  x > 0.  We  will  denote  functions  of  the 
1 n process  with  subscript  i,  i+1,2.  Let  the  regenerative  event  e^(x) 
be  the  event  which  occurs  at  time  n+1  if  Z^n+t)  = 0,  t c 0,  i-1,2. 

Since  r^(x)  ^ r2(x)  we  have  Z ^ ( t ) < Z2(t).  In  particular, 

Z2(t)  = 0 implies  Z^t)  = 0.  Thus  if  we  can  show  that  c2(t)  is  certain 
and  v«2  < ~ ^ have  shown  that  c ^ ( x ) is  certain  and  pj  < °°  . •'  From 
Leninas  2.2.1  - 2.2.3  we  know  that  e2(t)  is  certain  and  p2  < <*  if  we  can 
show  that 

Z { V2(u)  G-j  ( du ) + Z1  f Y2(u)  B(du,t)  X (t)  dt  < 1. 


This  follows  from  the  definition  of  Y(u).  Since 


V“>  ' l 7^) 


dX  = 1 / 


a 


dx 


a 


we  have  that 

Z l Y2(u)  G-j  (du)  + Z1  f Y2(u)  B(du.t)  X (t)  dt 

v 00  1 . 

" 1=1  o au6j(du)+  o1  o a u B(du,t)  X (t)  dt 


by  the  assumption  of  the  Lemma.  From  Theorem  3 of  Smith  (1955)  we  can  con- 
clude that  a quasi-limiting  distribution  exists. 
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2.5  A Limit  Theorem  for  Type  B Processes  with  General  Periodic  Input 

In  looking  at  a Type  A storage  process  we  have  the  advantage  of  being 
able  to  focus  on  the  times  the  store  is  empty,  a convenient  event  which 
is  regenerative  in  nature.  The  limit  theorems  of  Harrison  and  Resnick 
(1976)  require  assumption  (iv)  to  rule  out  any  process  which  cannot 
empty  for  this  very  purpose.  But  Type  B processes  are  quite  common  and 
one  would  hope  that  results  similar  to  Theorem  2.2  or  2.3  should  exist. 
And  indeed  such  is  the  case. 

The  next  theorem  is  for  a Type  B process  with  periodic  input--a 
partner  to  Theorem  2.2.  The  fact  that  a Type  B process,  while  not 
emptying,  can  have  content  arbitrarily  close  to  zero  is  used  to  set  up 
a "false  bottom"  near  zero  where  Theorem  2.2  can  apply. 


•< 


Theorem  2.4:  Suppose  Z(t)  is  the  content  of  a Type  B storage  system 
with  general  periodic  input  and  nondecreasing  release  rule  r(x).  Define 


Y(u,e) 


.u  dx 
4 iTx). 


If,  for  every  e 


V1c 


CO  1 O) 

s Y(u,e)  G-j (du)  + f f Y(u,e)  B(du,t)  X (t)  dt  < 1 


then  the  store  content  has  a quasi-limiting  distribution  which  is  inde- 
pendent of  the  initial  content  distribution. 


Proof:  Let  Z(t)  = Z(t,r,i;,r)  be  a Type  B process  as  described  in  the 
theorem.  In  addition,  define  a sequence  of  storage  processes 
Zn(t)  = Z(t,«o,f.,rn)  which  have  the  same  input  <o  and  initial  input  f,  as 
the  Z(t)  process,  but  release  rule  rn(x),  where 


1 
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\ 


* 

* . 

r(x)  - < r(x)  x > 1/n 
j 0 x < 1/n 

Note  that  each  Z^-  process  is  a Type  A process  with  "bottom"  1/n.  Since 
u i /n  < 1 the  assumptions  of  Theorem  2.2  are  satisfied  for  each  Zn»  pro- 
cess. Thus  we  can  conclude  that  each  Z - process  will  have  a quasi- 

*1 

limiting  distribution;  i.e.,  if  H(i(x,t)  • P {Z(i(t)  < x)  then 

lim  Hn(x,k+r)  * H (x,t)  where  H is  the  quasi-limiting  distribution. 
k^*> 

Since  {r^}  is  a nondecreasing  sequence,  by  lemma  2.1.3  {Zn(t)l  is 
non increasing.  Thus  lHn(x,t))  is  a nondecreasing  sequence  of  distribu- 
tion functions  for  each  t.  This  implies  that  IH^x.t)}  is  also  a non- 
decreasing sequence,  since  if  m^n  implies 

H?1(x,k+v)  < Hm(x,k+t) 

then 

**  «** 

Hn(x,t)  = lim  Hti(x,k+r)  < lim  H (x,k+t)  * H (x,t). 
k-*«  k-^v 

Thus  we  know  that  a limit  of  the  distribution  functions  exists.  Let 
H(x,t)  * lim  H (x,t) . 

n^v 

Since  {Hn(x,t)}  is  a nondecreasing  sequence  we  have  H^(x,t)  v Hn(x,t) 
for  all  n.  H^(x,t)  is  a distribution  function,  so  lim  H^(x,t)  ■ 1.  Thus 

X-*"' 

for  every  c there  is  an  A(  such  that  Hj(A  ,t)  > 1 -e.  Then  we  know  that 
Hn(A  ,t)  > 1-e  for  all  n.  This  means  that  iHn(x,T)l  is  a tight 
sequence,  since  for  all  n 

Hn(Ac,T)  - Hn(-Ac,r)  - Hn(A£,T)  > 1 - r. 
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Thus  we  can  conclude  that  H(x,t)  is  a distribution  function  and  Hn(x,x) 
converges  weakly  to  H(x,t). 

Next  we  show  that  for  each  n,  0 < Zn ( t ) - Z(t)  < for  all  t.  Let 

en  be  the  event:  en  occurs  at  t=k+l  if  Z (k+x)  = 1/n,  0 < r < 1.  Since 
t t n — 

Zn ( t ) is  the  content  of  a Type  A process  and  < 1 we  know,  from  Lemma 

2.2.3  that  for  some  t,  e is  certain.  Thus  if  we  let  T be  the  first 

t n 

time  Zn(t)  = ^ then  Tp  is  a proper  random  variable.  For  t < T^, 

Z(t)  = Zn(t).  Let  tj  be  the  time  of  the  first  input  after  T . For 
Tn  < t < t1  we  have  Zn(t)  = and  Z(t)  < ^,  so  0 < zn(t)  - Z(t)  < ^ for 

t < tj.  At  tp  there  is  an  input  of  size  U,  say.  Then 

W = U + Zn(tl_)  and  Z(tlJ  = U + Z(V)-  So 

o i Czn(tr)  + U]  - [Z(tr)  + U]  < I . 

But  this  is  equivalent  to 

0 <*„<»,)  -Zttjlii. 

Thus  for  t < tj 

(12)  0 < Zn(t)  - Z(t)  < I . 

Since  (12)  holds  both  between  inputs  and  directly  after  an  input,  we  can 
conclude  that  it  holds  for  all  t. 

From  (12)  it  must  hold  that  for  every  n 

(13)  P (Zn(t)  < x)  < P (Z(t)  <x)<  P (Zn(t)  < x + M . 

Let  H(x,t)  = P (Z(t)  < x)  . Then  (13)  can  be  rewritten 


i 
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1 

Hn(x,t)  i 1 Mn(*  + t)  for  all  t. 

This  implies 

(14)  H (x,t)  £ lim  H(x,k+T)  < liln  H(x,k+x)  < H (x  + -.x) 
n k-«°  n n 

For  any  e,  let  n be  large  enough  that  e > Then 

Hn(x  + -,x)  £Hn(x+e,x).  Letting  n go  to  infinity  in  (14)  we  get 

H(x,x)  £ lim  H(x,k+x)_<l im  H(x,k+x)  £ H(x+e,x) 
k-«o  k-»° 

Letting  e approach  zero  we  have 

H(x,x)  £ lim  H(x,k+x)  £ lim  H(x,k+T)  £ H(x+0,r) 
k-*»  k-*°° 

At  every  continuity  point  of  H we  have  lim  H(x,k+x)  = H(x,r). 

k-«» 

Thus  we  have  shown  that  H(x,k+x)  converges  weakly  to  H(x,x),  the 
desired  conclusion. 

2.6  A Limit  Theorem  for  a Homogeneous  Storage  System 

We  close  the  chapter  with  a theorem  concerning  homogeneous  storage 
processes.  The  result  was  originally  developed  as  a response  to  the 
paper  by  Harrison  and  Resnick  (1976),  whose  conditions  for  the  existence 
of  a limiting  distribution  are  not  at  all  obvious  or  intuitive.  Also, 
it  appeared  that  the  conditions  would  often  be  difficult,  if  not  impos- 
sible to  verify.  In  addition,  the  paper  unnecessarily  eliminated  a 
rich  class  of  processes--Type  B processes.  It  seemed  apparent  that 
some  other  level  than  the  zero  level  could  be  used  as  a basis  for  a 
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regenerative  event,  allowing  Type  B processes.  For  these  reasons. 
Theorem  2.5  was  developed. 

Theorem  2.5:  Let  Z(t)  be  the  content  of  a storage  system  with  homo- 
geneous  Compound  Poisson  input,  intensity  X,  and  release  rule  r(x) 
which  satisfies  the  following  criterion: 

There  exists  a level  c such  that 
r(x)  > « for  x > c 

r(x)  < M for  x < c 

— c — 

for  some  constants  a and  M . 

c 

If  the  first  moment  of  the  input  size  distribution  (p^)  is  finite  and 
a > \ then  there  exists  a limiting  distribution  of  store  content  that 
is  independent  of  initial  conditions. 

Proof:  For  the  store  content  Zj(t)  = Z(t,u,^,r^)  as  described  above, 
let  us  say  eA  occurs  at  t if  Z^ ( t)  = c and  Z^t-)  > c.  (i.e.  Z^t) 

crosses  c from  above  at  time  t).  We  do  not  consider  £=c  an  occurrence 

of  eA. 

A 

Clearly  is  a regenerative  event.  If  we  can  show  the  following 
three  criteria  are  satisfied  then  the  desired  limiting  distribution  will 
exist. 

a 

(Bl)  is  certain 

A 

(B2)  The  time  between  successive  occurrence  of  has  finite 
expectation 

(B3)  For  some  useful  class  of  sets  C,  4>c(t)(l-F(t))  is  of  bounded 

variation  in  every  finite  interval,  c e C,  where 

*c(t)(l-F(t))  = P{Z1(T^  + t)ec,tj>t  | A,T?}and  TA  is  the  time 

A 

of  the  first  occurrence  of  e.. 
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If  c=0  then  we  have  r^x)  > a > for  x > 0.  Define 

Z3(t)  = Z(t,u>,£,r3)  where  r3(x)  = a < r^x).  From  Lemma  2.1.1 
Z^(t)  £ Z3(t).  Thus  if  Z3(t)  = 0 then  Zj(t)  = 0.  From  standard  storage 
theory  (see  Appendix)  we  know  that  the  Z3  - process  will  empty  with 
probability  one  and  the  expected  time  between  emptying  is  finite.  Thus 
we  have  (Bl)  and  (B2)  satisfied  for  the  case  c=0. 

To  prove  ( B 1 ) - ( B2 ) for  c > 0 we  will  need  to  define  some  additional 

D 

events  and  random  variables.  Let  us  say  occurs  at  t if  Z1(t-)<c  and 
Zj(t)>c.  Let  Uj  be  the  time  from  an  to  the  next  occurrence  of 
Vj  be  the  time  from  an  ej  to  the  next  occurrence  of  and  Wj  be  the 
time  an  to  the  next  occurrence  of  e".  Define  Z2(t)=Z(t,w,£,r2)  where 
r2^x)=^c’  x > 0*  All  the  above  events  and  random  variables  apply  to  the 
Z2  - and  Z3  - process  (defined  above),  denoted  by  appropriate  subscripts. 
Note  that  the  Z2  and  Z3  systems  have  homogeneous  Compound  Poisson  input 
and  constant  release  rule.  Results  on  such  processes  are  well  known. 

In  particular,  we  know  that  Z2( t)  and  Z3(t)  are  type  A processes  and  can 
empty.  Thus  for  Z2(t)  and  Z3(t)  we  can  introduce  the  event  which 
occurs  at  t if  Z.(t)=0  and  Z.(t-)>0,  j=2 ,3 . The  random  variable  X.  will 

J J J 

be  the  length  of  the  busy  cycle  in  the  Z.  - process,  j=2,3. 

J 

A A 

To  begin,  we  show  that  is  certain.  Suppose  £ > c.  For  t <_  Tr 

we  have  Z3( t ) >_  Z^t),  since  r^x)  > a = r3(x).  Thus  will  be  certain 
A A 

if  e3  is,  and  e3  will  be  certain  if  e3  is.  But  A Bj  < r3(x)  = a implies, 

a 

by  standard  storage  theory,  that  e3  is  certain.  Thus  we  have  that  Cj  is 
certain  if  £ > c. 

A B 

Let  us  next  suppose  £ < c.  The  event  will  be  certain  if  is 

D 

certain  and  Vj  is  a proper  random  variable.  Let  T°  be  the  time  of  the 
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first  occurrence  of  e^.  For  t < t!|  we  have  Z^(t)  < Z^(t).  This  is  true, 
since  while  Z^t)  < c,  r^x)  < Mc  » r2(x).  Thus  e®  will  occur  before  c® 
does,  which  implies  that  c®  is  certain  if  cjj  is. 

Let 

p = PfZ^t)  jumps  above  c during  a busy  cycle  of  the  Z^  - process) 

p = P{Z^(t)  jumps  above  c during  the  initial  busy  cycle  of  the 
Z^-process ) 
q = 1 - P,  q 3 1 - P 

If  we  assume  p > 0 then 

P{T^  < ■»)  = Ptc^  first  occurs  during  the  initial  busy  cycle) 

00  n a.  i 

+ ^ Pft^  first  occurs  during  the  k busy  cycle) 

= p + q p + q qp  + q (q)  p + ... 

=p+qp(l+q+q  +...  ) 

3 P + q P/P 
= 1 

i.e.,  if  p>0  then  the  first  occurrence  of  clj  is  certain.  Now,  in  any 
interval  of  length  y we  know  that  the  output  of  the  7^  ~ process  must 
be  less  than  or  equal  to  Mcy,  since  r^(x)  = M^..  We  know  that  with  Com- 
pound Poisson  input  there  is  a positive  probability  of  input  greater 
than  c + Mc  y.  Thus  we  have  pH)  and  i is  certain,  implying  is 
certain. 
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g 

Note  that  Z^(T. ) > c.  Thus  showing  that  Vj  is  a proper  random 

A B 

variable  is  equivalent  to  showing  that  is  certain  if  ^ = Z1  (T"), 

A 

which  we  have  already  shown.  Hence  we  have  that  is  certain. 

To  show  that  (B2)  is  satisfied  we  need  the  following  lemma. 

Lenina  2.5.1:  Let  Z^(t)  be  the  content  of  a storage  system  as  described 

A 

in  Theorem  2.5.  Let  occur  at  t if  Z^(t)=c  and  Zj(t-)>c,  and  Wj  be 

A 

the  time  between  successive  occurrences  of  Cj.  If  P2  < " then 

E(Wj)  < ». 

Proof:  We  employ  the  same  processes  and  random  variables  as  in  the  proof 
of  Theorem  2.5.  Define  further  the  following  random  variables: 

n is  the  length  of  a busy  cycle  of  Z^( t ) during  which  the  store 
content  does  not  cross  level  c during  the  busy  cycle, 
t is  the  length  of  a busy  cycle  of  Z2(t)  during  which  the  store 
content  does  cross  level  c during  the  busy  cycle. 

Since  Wj  = Uj  + we  need  to  look  at  and  Vj.  Suppose  £=c.  Then 
U1  - U2’  Thls  true  since  when  ^j(t)  < c.r^x)  > r2(x),  so  Z2(t) 
will  take  at  least  as  long  as  Z^(t)  to  rise  above  c.  Now 

U2  £ Dj  + n2  + • • • + nN  + c = 0 

where  N,  a random  variable,  is  the  number  of  times  the  process  does  not 
cross  level  c during  successive  busy  cycles,  before  finally  crossing  c. 
Let  S|y  = + n2  • • • + Hn»  so  that  0 = S^  + c,  a random  sum.  Note  that 

n.t  and  N are  mutually  independent.  Using  Walds  equation 
E(O)  - E(Sn)  ♦ E(c) 


E (N)E(n)  + E(c) 
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We  also  have 

E(02)  - E(S2)  + 2E(r,SN)  + E(c2) 
Since  t.  n and  N are  independent 
E(cSn)  + E(c)E(N)E(n) 

What  is  E(S2)? 


E(sJ)  - e[e(S^  |N)] 


■ & p<N'k>E<sk> 
k 


k=l  k 


" k^o  EtnJ)  * 2 }lx  ,£jtl  Etnjn,)) 

- kI,  P(N=k)  (k£(nZ)  + k< k- U [e< n) ]2} 

= E(n  ) klj  kP(N=k)  + (Eor  k^P(N=k)-k^kP{N=k) } 
■ E(n2)E(N)  + (En)2  [e(N2)-£(N) ] 


Thus 


E(02)  = E(n2)E(N)  + (En)2  [e(N2)  -E(N) ] + 2E(c)E(N)E(n)  + E(r.2) 

L, 

We  know  that  P(N=k)  5 q p,  where  p and  q are  as  defined  in  the  proof 

of  Theorem  2.5.  Thus  E(N)  * q/p  < •"  if  p^O.  But  in  the  proof  of 

Theorem  2.5  we  showed  that  p-'O.  Both  E(n)  and  E(0  will  be  finite  if 

E(X^)  < « and  p>0,  since  E(X^)  = q E(n)  + p E(c).  By  standard  storage 

theory  results,  E( X?)  < *>  if  Bj  < «*  and  X Bj  < , both  of  which  are  true 

2 2 

by  assumption.  Similarly,  E(n  ) and  E(t  ) will  be  finite  if  B.,  < ■”  and 

X B,  < M . Thus  we  have  that  for  f.=c 
1 c 
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E(Uj)  < E(U2)  < E(0)  < - 
E(uJ)  < E(u|)  < E(02)  < « 

Suppose  &c  and  define  Z^*(t)  = . We  know  that 

* * 

E(Uj  ) < °°,  where  Uj  is  as  defined  earlier,  but  with  reference  to  the 

* . * * 

Z1  - process,  rather  than  Z^(t).  Note  that  Z^  (U^  ) > c implies  that 

Zj(T^  + ui  ) > c»  i*e*  ui  1 ui**  Thus  we  have  that  E(Uj)  £ E(Uj  ) 
for  all  £. 

We  still  need  to  consider  V^.  Suppose  £;>c.  We  have  shown  earlier 
that  Vj<Vj,  and  clearly  we  know  that  V.j<c.  Thus  we  have  that 

E(Vj)  < E(c) 

E(V2)  < E(^2) 

and  we  know  that  these  both  will  be  finite  since  Bj  < 00 ,3^  < 00  and 

X Bj  < a. 

Next  suppose  we  drop  the  assumption  that  £=c.  Define 

Z**(t)  = Zj(t+T® .w.S.rj).  Then  Zj*(0)>c.  We  know  that  E(Vj  )<®where 

**  **  ,, 
is  as  defined  earlier,  but  with  reference  to  the  Z^  - process.  Note 

that  Z^  (Vj  )<c,  which  implies  that  Z^(V^  + T®)<c,  i.e.  . Thus 

we  have  E(Vj)  £ E(Vj  )<°°for  all  £;. 

Since  Wj  = Uj  + Vj  we  have  E(Wj)  = E(l!j)  + E(Vj)  < » and  since  Uj  and 

are  independent 

E(W2)  = E(U2)  + E(V2)  + 2E(U1)E(V1)  < - 


and  the  lemma  is  proved. 
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To  finish  the  proof  of  the  theorem,  we  must  prove  (B3).  This  can  be 
done  using  the  following  adaptation  of  Lemma  2 of  Smith  (1955). 

Lenina  2 (Smith):  If  (i)  I = (to+a,  tQ+b]  is  some  finite  t- interval, 
and  if  y.  is  a nonnegative  random  function  of  half-closed  subintervals 
I.  of  I with  the  properties 

J 

(a)  Yt  + yT  * y,  , for  any  adjacent  subintervals  I,,  I9cl  I 

ll  *2  1 1 *2  1 1 

(b)  with  probability  one  E (y j | tQ , Z(0))  < » 

(ii)  With  probability  one,  conditional  upon  tQ,  Z(0),  for  any  I • c.  I , 

ndjO  ?sA)  < yj 

Then  ^(t)  = <j>.(t)(l-F(t))  is  of  bounded  variation  in  (a,b). 

For  the  above  lemma,  t is  the  time  of  the  first  occurrence  of  some 
event  t,  A is  some  useful  class  of  sets,  S^  = (t:  Z(t)iAl,  f S^  denotes 
the  frontier  or  boundary  of  S^  andTl  denotes  the  cardinality  of  a set. 
The  original  lemma  required  that  (a)  hold  for  any  disjoint  intervals 
1 j , l^cLl,  but  since  the  additivity  property  was  only  used  in  a binary 
manner  over  the  disection  of  I,  the  property  need  only  hold  for  adjacent 
subintervals.  Since  we  are  only  concerned  with  showing  that  the  distri- 
bution of  store  content  has  a limit,  we  can  restrict  ourselves  to  the 

class  of  sets^r  which  consists  of  sets  A = lx:  x • c). 

c 

In  order  to  prove  (B3)  for  the  class  we  must  find  a set  function 
y.  which  satisfies  (i)  and  (ii)  of  Lemma  2 above.  As  suggested  in  Smith 
(1955)  we  define  the  random  variable  it.  as  follows: 

For  any  interval  (ij.j*  tj]  = 


*>.  7 


if  Z(t.  ,)  eA  and  Z(x.)  eA 

J * J 

or  Z(tj_1)  *A  and  Z(x.)  {A 

if  Z(xj_j)  {A  and  Z(t.)  cA 

if  Z ( x j _ j ) eA  and  Z(Tj)  $A 


We  also  define  a second  set  function  nT  as  the  number  of  inputs  into  the 

‘j 

system  in  the  interval  I.,  and  take  y.  = 2 n.  + 6.  . 

J 

Additivity  is  clear  for  n.  . To  see  this  for  6.  it  is  somewhat 
easier  to  think  of  5.  as  a sum  of  two  functions  a.  and  6.  where  for 

'o  = (Tj-rTj] 


-1/2 

if  tj-i  1 

1/2 

if  Vi 1 

1/2 

if  Xj  eA 

-1/2 

,f 

Clearly  6.  = a,  + $T  . For  additivity  we  examine 

l3  J 


So  we  have  6T  + 6.  = 6r  llT 

Li  l3+l  J+l 


That  (b)  is  true  is  obvious,  since  Ey.  < 2 E n.  +1  and  E(n.)  < « 
since  we  are  dealing  with  Compound  Poisson  input. 

To  show  that  (ii)  is  satisfied  we  must  determine  how  many  boundary 
points  there  are  in  S.  which  are  also  in  I..  Let  A = [0,c].  Then 

n J C 

Z(t)  e Ac  if  and  only  if  Z(t)  £ c.  Note  that  if  Z(t)  e Ac  then 
Z(t+s)  eAc  if  there  are  no  inputs  to  the  system  in  (t,t+s].  If  there 
is  an  input  there  are  two  possibilities 

(1)  The  process  remains  in  Ac  (hence  creates  no  boundary  points) 

(2)  The  process  jumps  above  c.  Then  we  have  a boundary  point,  and 
possibly  another  if  the  process  reenters  Ac> 

Thus  an  input  during  I.  = (x.  ,,x.]  results  in  creating  at  most  two 

J J ^ J 

boundary  points  of  SA  --  the  first  if  the  input  is  large  enough  for  Z(t) 

to  jump  above  c,  followed  by  a second  if  the  process  falls  back  below  c 

before  Xj.  Let  = ( tk_ ^ , t^]  where  t^,  t2,  ...  , t^  are  times  where 

inputs  occur  in  I , tQ  = Xj_j,  t^+1  = x^.  and  N = nj  . (Note,  if  tN  =Xj 

then  we  say  JN+^  = <J>.)  We  can  write  Ij  = J^.  CertainlyT\(I^A3:'SA)  = 

?Sa  ).  We  have 
k=l  K Mc 

?SA  ) < J 1 if  z(xj.1)  eAc 

1 2 if  Z(Tj.l)  ^AC  j 


-n(Jkn  tsa  ) < 2 


k = 2,  ...  , N 


"tt(JN+inTSA  > i J 1 if  Z^Tj)  cA 


0 if  Z(Tj)  ^Ac 


i 
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Thus 


• 

1+2(N-1)+1 

if 

z(Tj_l)  eAc  and  Z(t.)  eAc 

n(  ij0  3-sA)< 

2+2(N-l)+0 

if 

z(tj-_1  ) and  Z(t j ) $AC 

2+2(N-l)+l 

if 

z(^j_l ) <^AC  and  Z(x j ) eAc 

l+2(N-l)+0 

t> 

if 

z ( T j _ 1 ) eAc  and  Z(t..) 

= 

2nj  + 

3 J 

We  have  shown 

that  y.  = 2 n. 

+ 6. 

satisfies  the  assumptions 

Lemma  2 of  Smith  (1955),  which  implies  that  <t>^( t)  (l-F(t)}  is  of  bounded 
variation  for  sets  A of  the  type  [0,c].  This  is  sufficient  to  complete 
the  proof  of  our  theorem,  since  with  (B1)-(B3)  satisfied,  theorem  2 of 
Smith  (1955)  states  that  the  limiting  distribution  of  Z(t)  exists.  In 
fact,  the  theorem  specifies  that 

Hm  P <Z(t)  < c)  = I f 4>.  (v)  (l-F(v))  d v 

t-*»  U O « 

where  A = [0,c]. 

Sometime  after  this  theorem  was  completed  the  paper  by  Brockwell  (1977) 
appeared.  His  limit  theorem  is  similar  to  Theorem  2.5,  but  Brockwell 
assumed  from  the  outset  that  r(x)  was  nondecreasing  and  continuous,  while 
this  theorem  does  not. 


CHAPTER  3 


PERTURBATION  APPROACH  TO  THE  STORE 
WITH  NONHOMOGENEOUS  INPUT 

3.1  A Probabilistic  Method 

In  this  chapter  we  present  an  approach  which  will  allow  us  to  develop 
formulae  for  various  functions  of  a storage  system  with  time  dependent 
input.  The  method  is  direct  and  probabilistic  in  nature,  making  use  of 
the  structure  of  the  input  process. 

We  will  be  interested,  from  this  time  forward,  in  a process  with 
general  release  rule,  where  the  number  of  inputs  into  the  system  in  any 
interval  (0,t]  is  a nonhomogeneous  Poisson  random  variable,  i.e. 

P{ 1 input  in  (t,t+dt)}  = A(t)dt  + o(dt) 

P{more  than  1 input  in  (t,t+dt)}  = o(dt) 

Let  N(t)  be  the  number  of  inputs  in  (0 , t] . Then  we  have 

(1)  P{N(t)  = n}  = [A(t)]n  e'A(t)  ; A( t)  = A(u)du 

n!  0 

Note  that  if  x(t)=A,  then  A(t)=At  and  (1)  is  just  the  familiar  homogene- 
ous Poisson  probability.  If  we  let  B(x)  be  the  distribution  function  of 
input  size,  then  the  total  input  in  (0 , t] , A(t),  has  nonhomogeneous  Com- 
pound Poisson  distribution,  with 

P{A(t)  < x)  = f e“A(t)  rA(t)1n  B*"(x) 
n‘°  n! 

where  B*n  is  the  n-fold  convolution  of  B with  itself. 
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The  key  to  the  development  of  our  approach  will  be  to  assume  that 
the  intensity  X(t)  « Xo(l+e<t>(t))  where  XQ>0,e>0  and  <t>(t)  is  some  meas- 
urable nonnegative  function.  The  resulting  input  process  can  thus  be 
regarded  as  the  sum  of  two  independent  Compound  Poisson  processes:  one 
with  intensity  X1(t)=XQ;  the  other  with  intensity  X2(t)  = XQe4)(t).  We 
will  call  the  inputs  into  the  system  originating  from  these  two  sources 
Type  I and  Type  II  inputs  respectively.  Notice  that  the  Type  I input  is 
time  homogeneous.  This  will  play  an  instrumental  role  in  the  develop- 
ment of  formulae,  since,  conditional  on  knowing  when  the  Type  II  input 
occur,  the  storage  process  has  homogeneous  input. 

In  what  follows,  we  will  focus  our  attention  on  the  probability  of 
emptiness,  P{Z(t)=0}.  The  method  developed  here  can  be  applied  directly 
to  other  quantities,  such  as  the  mean  store  content  and  Laplace-Stieltjes 
Transform  of  store  content.  Let  N - ( t ) be  the  number  of  inputs  of  Type  j 

J 

which  arrive  at  the  store  in  (0,t],  j=l,2.  Then  the  probability  the 
store  is  empty  can  be  written 

P{Z(t)=0}  = kf0  P{Z(t)=0  and  N2(t)*k}. 

We  will  examine  this  series  term  by  term.  To  begin,  we  note  that 

P{Z(t)=0  and  N2(t)=0}  = P{Z(t)=0  | N2(t)=0}  P{N2(t)=0} 

If  we  know  there  have  been  no  Type  II  inputs  up  to  time  t then  the  pro- 
cess (up  to  that  time)  is  equivalent  to  the  homogeneous  store  with 
intensity  x . Let  n(t,x(*).W)  denote  the  probability  the  store  is 
empty  at  t,  for  a store  with  Compound  Poisson  input,  intensity  A(-)*  and 
initial  content  distribution  W.  Then  we  have 

p{z(t)=o  | N2(t)=o)  = n(t,x0,w0) 


L 
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and  thus 

-eX  d>  ( t) 

P{Z(t)=0  and  N2(t)=k}  = n(t,x0,W0)e  0 

t 

where  <t>(t)=  £ ,j>(u)du  and  WQ(x)  Is  the  distribution  function  of  Z(0). 

Next,  suppose  N2(t)=l.  The  Type  II  input  will  occur  at  some  point  u 
in  the  interval  (0,t].  Let  T.  be  the  time  of  the  ith  Type  II  input. 

Then 

t 

P{Z(t)=0  and  N2(t)=l>  = f P{Z(t)=0  and  N£( t)=l  and  TjedUj} 
which  is  just 

^P(Z(t)=0  | N2(t)=l  and  T^Uj)  P(N2(t)=l  and  T^du^. 

The  probability  of  a Type  II  input  in  dUj  is  e*0<MUj)dUj  + o(dUj)  and 
the  probability  that  this  is  the  only  Type  II  input  in  (0,t]  is 
U1  t t / 

-O  CV(),>dX  -i,EAo^X)dX  -0  CV  X dX 

e eXo<))(u1)du1e  1 + ofdu^eX^u^e  du^  + o(du1) 

We  need  to  determine  the  probability  the  store  is  empty  at  t,  if  we 
know  that  the  first  Type  II  input  occurs  in  dUj  and  no  more  Type  II  in- 
puts occur  before  t.  The  input  up  to  time  Uj-  is  all  homogeneous  (Type 
I)  input,  so  at  Uj-  the  distribution  of  store  content  is  W(x,Uj-,^o,WQ) , 
where  W(x,t,X,W)  denotes  the  distribution  at  time  t of  the  content  of  a 
store  with  Compound  Poisson  input,  intensity  X,  and  initial  distribution 
W(x).  At  time  u^+  there  is  an  input  to  the  store,  and  the  size  of  the 
input  is  independent  of  the  content  at  Uj.  Thus  the  distribution  of  con- 
tent at  Uj+  is  the  convolution  tf(x,u,Xo,WQ)*B(x).  Let  W(x,t,X|u)  be  the 
content  at  time  t of  a homogeneous  store  which  has  Compound  Poisson 
Input  with  intensity  X,  conditional  on  knowing  that  a Type  II  input 
occurs  in  du.  Then  we  can  write  W(x,Uj ,*0luj)  = W(x,Uj ,^0,W)*B(x) . 


We  use  the  shorthand  notation  W^(u^)  for  this  distribution  function  at 
time  Uj. 

Since  there  are  no  more  Type  II  inputs  in  (Uj,t],  the  probability 
the  store  will  be  empty  at  t is  just  the  probability  that  a store  with 
homogeneous  Compound  Poisson  input,  intensity  \ , and  initial  distribu- 
tion W(x,Uj,Aq|Uj)  will  be  empty  t-Uj  time  units  later,  or 
Il(t-u1,\0,W1(u1)).  We  are  therefore  led  to  the  equation 

-e\  <f(t)tt 

P{Z(t)=0  and  N^( t)  = l } = e t-uJ  ,A0.W1  (Uj ) )duj . 

In  a similar  fashion,  we  can  see  that  for  N^t^k,  if  the  k Type  II 
inputs  occur  in  dUj,  du^,  ...  , duk  then 


P{Z(t)=0  and  N2( t)=k } = 


U1 

/ x ) dx 


"c'o  4-,-'(x)dx 

•e  cAo«(“k^ 


t 

-C.\0  {.  *(x)dk 
k 

e n(t-uk,A0,Wk(uk))du1. 


■e\0'T(t)  ^ 

= e CV  o o l ♦(“1)--^«k)n(t-uk.A0.Wk(uk))du1...duki 

The  notation  Wk(uk)  represents  the  distribution  function 
W(x,Uk ,AQ I ui * •••  »uk)  where  W(x,t,A|Uj,  ...  un)  is  the  content  at  time 
t of  a store  which  has  homogeneous  Compound  Poisson  input  with  intensity 
\0,  conditional  on  knowing  that  Type  II  inputs  occur  in  dUj,  du0,  .... 
du^.  We  assume  that  W(x,0,\o|u^,  ...  ,u^)  = W^(x)  for  all  n. 


59 


We  have  thus  derived  the  following  formula  for  the  probability  a 


store  will  be 

n(t,A(.),wo)=e 

(2) 


empty  at  t: 

t 

Ol(t,A0,Wo)+eA0  / <l»(u1)n(t-u1,Ao>W1(u1))du 


uk  u2 

o ***o  • • •<f’(u|<)n(t-uj, ,Aq»W|i,(u|i,) )dUj . . .du^} 


1 


An  alternate  form  of  (2)  expresses  n(t,A( -) .WQ)  as  a power  series  in  e. 
To  do  this,  simply  note  that 

-eA  4>(t) 

e = jl0  (-l)J[eA0«l»(t)]J 

If 

t uk  u£ 

Then  If  we  let  Ik(t)=  ' ' ' ) . . . n( t-ufe ,XQ .«k( ufc ) )dUl . . .dufc  _ 

k>l,  I0(t)*n(t,A0,W0)  and  Ij(t)=  /t<K(u1)n(t-u1,A0,W1(u1))du1,  we  can 
rewrite  (2)  as 

(3)  n(t,A(-).w0)  - j|o  nj(t)ej 

where 

Rj(t)  = Xo  Jo  (_1)k  *k(t)  I(t) 


Similar  equations  to  (2)  and  (3)  can  be  written  for  such  quantities  as 
the  distribution  of  store  content  and  mean  content. 

_3.2  Extending  the  Approach  to  General  aft) 

All  the  work  so  far  in  this  chapter  has  required  that  4>(t)  be  a non- 
negative function,  since  e)^<j>{t)  was  taken  to  be  the  intensity  of  a 
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Compound  Poisson  process.  But  this  assumption  is  restrictive  and  unde- 
sirable in  general.  Unfortunately  we  have  been  unable  to  prove  a general 
theorem  which  applies  to  a wide  range  of  processes.  We  feel  that,  at 
the  very  least,  the  following  can  be  shown  to  hold. 

Conjecture  3.1:  For  a storage  system  with  nonhomogeneous  Compound 
Poisson  input  which  has  intensity  X ( t ) =\o ( l+ed> ( t ) ) and  a release  rule 
r(x)  which  is  a finite  degree  polynomial  in  x,  the  Laplace-Stieltjes 
Transform  of  store  content  can  be  written 

-cX  <J>(t)  t 

M(s  ,t,X(  • ) ,WQ)=e  0 {M(s,t,X0,W0)+eX0  ^ 4>(u1)M(s,t-u1,X0,(u1)  du 

(4) 

" t uk  u2 

+ k=2  e\  ooo  4>(u1)...4>(u|<)M(s,t-uk,X0,Wk(uk))du1  ...  dufc} 

where  M(s ,t,x(* ),W)  is  the  Laplace-Stieltjes  Transform  of  the  content  Z(t) 
of  the  store  with  intensity  X(-)  and  initial  distribution  W. 

Outline  of  Proof:  We  can  show  that  M(s ,t,X(* ),WQ)  satisfies  the  following 
partial  differential  equation: 

3 M(s ,t,X( • ), W)  + X(t)(l-B*(s))  M(s  ,t  ,X( • ) , W)  = 
at 

(5) 

m 

s kl0  (-Dkak  3kM(s,t,X(-),W)  - a0sn(t,x(-),w) 


ffl  I- 

where  r(x)  = k£Q  akx  . We  do  this  in  the  following  Lemma,  which  we  will 
find  useful  later  also. 


Lemma  3.1.1:  For  a storage  system  with  nonhomogeneous  Compound  Poisson 
input  which  has  intensity  X(t)  and  release  rule  r(x),  M(s,t,X( • ),VI) 
satisfies  the  equation 

(6)  3M(s ,t,A( • ),  W)  + A(t) (1-B*(s) )M(s ,t,A( • ) ,W)  = sE[r(z(t)  )e'SZ(t)]. 
at 

m k 

In  addition,  if  r(x)  = k£  a^x  then  (6)  becomes 

3M(s,t,X(-),W)  + X(t) ( 1-B*(s) )M(s ,t,A( • ),  W) 

3t 

m m l 

= S k^0  (-1)  ak3K  M(s,t,x(-),W)  - a0sn(t,X(-),W). 


Proof:  We  need  to  consider 


3M(s,t,X(-),W)  = lim  M(s ,t+x  ,X( • ) »W)  - M(s ,t,X( • ),  W ) 

x-K)  


Let  Aq  = AQ(t,x)  = {oo:N  ( t ,x  ,co) 
Aj  = A^(t,x)  = {w:N(t,x  ,u) 
A^  = A2(t,x)  = (a) : N(  t ,x  .to) 
where  N(t,x,io)  is  the  number  of 
the  parameters  X and  W from  the 
M(s,t+x)  as 


= 0} 

= 1} 

> 1} 

inputs  to  arrive  in  (t,  t+x].  Dropping 
notation  for  the  present,  we  write 


M(s,t+x)  = / e_sZ(t+T)dP 


I 
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= fl/e"sZ^t+T^dP  + ,/e"sZ^t+T^dP  + fl/e"sZ^t+T^dP 
o M1  a2 

= I0  + Ij  + I2.  say. 

We  will  examine  each  integral  separately.  Note  that  for  small  x, 
the  output  of  the  system  in  (t,t+x]  is  r(Z(t))].  If  there  is  no  input 
in  ( t , t+x ] then  the  content  at  time  t+x  can  be  written 
Z(t+x)  = Z(t)  - r(Z(t))i  if  x is  sufficiently  small.  For  such  x 

t = / -sZ(t)+s  r(Z(t))x.p 

o a' 

= A/  e"sZ(t)dP  + A /sxr(Z(t))  e"sZ(t)dP  + o(x) 
o o 

(8)  =i2XA  e"SZ(t)dP+^XA  sxr(Z(t))  e'sZ(t)dP  + o(x) 

O 0 

where  x A is  the  indicator  function  of  the  set  A.  Since  AQ(t,x)  is  inde- 
pendent of  Z(t),  (8)  can  be  written 

(9)  lo  = p^0(t’'r)}M(s,t)+sxP{A0(t,x)}E{r(Z(t))e‘sZ(t)}+o(x) 

For  a process  with  nonhomogeneous  Compound  Poisson  input,  if  x is  small, 
P{AQ(t,x)}  = (l-A(t)x)  + o(x).  Thus  we  have 

(10)  IQ  = M(s,t)-A(t)xM(s,t)+sxE{r(Z(t))e'sZ(t)}+o(x) 

Suppose  there  is  an  input  in  (t,t+i]  of  size  U.  Then  the  content  at 
time  t+x  can  be  written,  for  small  i , as  Z(t+i )=Z(t)+IJ  r(Z(t))  i.  We  use 
this  to  write  1^  as 


I 


! = A /e- sZ(t)-su+STr(Z(t))dp 

(ID 

= ./  e"sZ(t)"sudP  + sx  . / r(Z(t) ) e"sZ(t)_sudP  + o(t) 

1 A1 

Since  AjU.xJ.U  and  Z(t)  are  mutually  independent,  (11)  becomes 

Ij  - P{A1(t,r)}[B*(s)M(s,t)+B*(s)sr^r(Z(t))e'sZ(t)dP]  + o(r) 
Recall  that  PtA^t.x)}  = X(t)x  + o(t).  Thus 

(12)  Ij  = X(t)xB*(s)M(s,t)  + o(x) 

The  third  integral,  I2,  is  nonnegative  and  bounded  above  by  a func 
tion  of  order  o(x).  We  see  this  as  follows: 

(13)  I2  = A/e"sZ(t+T)dP  < /ldP  = P{A2(t,T)}  = o(x). 

We  can  combine  (10),  (12),  and  (13)  to  show  that 

3M(s,t)=  -X(t)M(s ,t)+s  E{r(Z(t))e'sZ(t)}  + x(t)B*(s)M(s,t) 

3t 

°o  k - 

Suppose  r(x)  = a^x  , and  let  r(x)  = r(x)  for  xx).  Then 

E{r(Z(t))e"sZ(t)}  = ^r(Z(t))e"sZ(t)dP 

= /r(Z(t))e‘sZ(t)dP 
(Z(t)>o) 
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= Q/r  (Z(t))e'sZ{t)dP  - r(0)ll(t) 

■ Jo  *k  f!Zk(t)e-sZ(t)dP  - aon(t) 

m 

(14)  = k^0  aR(- l)k  3Mk(s,t)  - aQn(t) 

3s 

Substituting  (14)  into  (5)  gives  the  desired  conclusion. 

Outline  of  Proof  continued:  We  also  know  that  (4)  holds  for  ^(t)No. 

This  can  be  shown  by  using  the  perturbation  argument  presented  earlier 
in  this  chapter.  Thus  M(s,t,A(* ) ,W)  as  defined  in  (4)  must  satisfy  (5). 
It  seems  clear  that  (4)  will  formally  satisfy  the  partial  differential 
equation,  regardless  of  the  range  of  values  assumed  by  <5>(t).  If  we  can 
show  that  (5)  has  a unique  solution  which  satisfies  the  boundary  condi- 
tions 


M(0,t,A( • ) ,W)  = 1 

(15) 

*x> 

M(s ,0 ,x( • ) »W)  = fQ  e"sxdW(x) 

then  we  will  have  shown  that  (4)  is  the  Laplace-Stiel tjes  Transform  we 
seek. 

For  the  store  with  r(x)  l,x'0,  Reich  (1958)  has  shown  that  there  is 
a unique  probabilistic  solution  to  (5).  For  higher  order  polynomials, 
we  have  been  unable  to  find  appropriate  criteria  for  the  existence  of  a 
unique  solution. 

Corollary  3.1.1:  For  a store  as  described  in  Conjecture  3.1  the  follow- 
ing hold: 
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-eA4>(t)  t 

(16)  n(t,X(-),W0)=e  {n(t,A0,W0)+eA0  o 4>(ui)n(t"ui»x0»wi(ui)  dui 


+ k=2ek*o  o o”k’ o — 4>(uk)n(t-uk ,X0,Wk(uk)  JdUj. . .duk> 


-eA  4>(t)  t 

(17)  y(t,x(- ),wQ)=e  {y(t,A0,Wo)+eA0  o 4>(u1)y(t-u1,X0,W1(u1))du1 


. . t u.  u0 

k,  k r r k r 2 


+ k=2eXoo  o—o  4>(u1)..^(uk)M(t-uk,xo,Wk(uk))du1...duk} 


where  y is  the  expected  store  content. 

Proof:  Recall  that  lim  M(s,t,x(*),W  )=n(t ,x( - ) ,W  ).  Taking  the  limit 

s-*»  0 0 

as  s goes  to  infinity  in  (4)  yields  (16).  Note  also  that 

lim  9M(s ,t ,X( • ) »W  )/3s=y( t ,x( • ), W ) . Hence  using  (4)  we  can  get  (17). 
s+o  0 

Some  comments  are  in  order  about  our  chief  assumption— that  the  in- 
tensity of  our  input  process  can  be  written  in  the  form  A(t)=A0(l+o}>(t)). 
At  first  glance  this  might  be  thought  to  be  rather  restrictive.  Actually, 
it  is  not.  Suppose  we  have  an  intensity  x( t ) >0.  For  any  choice  of  XQ 
and  e we  can  define  <t>(t)  as  <p( t)= [X ( t ) -X0]/eXQ. 

Of  course,  more  suitable  choice  of  XQ,e  and  <j>  can  be  made.  If 
A(t)  is  periodic,  we  would  want  XQ  to  be  the  average  intensity  per 

period.  Then  A = y“x(u)du  and  <f>(t)=A(t)~  ywA(u)du 
o o 

c£*\(  u)du 

where  £>  is  the  period.  If  we  know  that  A(t)-*-L  as  t-*»  then  a good  choice 


for  A_  would  be  L. 
o 


3.3  An  Exact  Limiting  Result 

It  is  well  known  that  in  general,  the  transient  probability 
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n(t ,X0,Wk)  is  difficult  to  calculate.  Thus  determining  IT(t,A(  • ) ,WQ) 
will  be  even  more  difficult,  as  it  involves  an  integral  of  which 
n(t>A0,w|<(uk))  is  just  a part.  For  this  reason,  we  will  concentrate  on 
results  concerning  the  limiting  value  of  II. 

Suppose  we  know  that  X ( t ) converges  to  a finite  limit.  In  particular, 
let  A(t)=Ao(l+e<!>(t))->-Ao.  We  might  expect  intuitively  that 
n(t,A(  • ),  Wo)-*-n(°°. Aq).  With  only  some  minor  assumptions  on  4>(t)  we  can 
show  that  this  is  indeed  the  case.  Theorem  3.2,  which  proves  this,  is 
essentially  an  extension  of  part  of  Theorem  3 of  Takacs  (1955)  for  the 
M/G/l  queue. 

Before  proving  the  theorem,  we  will  first  prove  a few  lemmas  which 
will  be  useful  both  here  and  elsewhere.  The  first  shows  that  if  the  in- 
tensity of  the  input  process  is  bounded,  then  the  probability  of  emptiness 
and  the  content  moments  are  bounded. 

Lemma  3.2.1:  Suppose  we  have  a storage  system  which  has  nonhomogeneous 
Compound  Poisson  input  with  intensity  A(t)  where  Mj(t)<A(t)<M2(t)  for 
some  positive  functions  M1  and  M2.  Let  (t)=Z(t,<^  ,£ ,r)  denote  the 
store  content  with  intensity  M^t),  i = 1 ,2  and  Z(t)=Z(t,u),£,r)  denote 
the  content  with  intensity  A(t).  Then 

P(Z2(t)=0}<  P{Z(t)=0}<  PtZjC t)=0} 

E{Zk( t) } < E{Zk(t)}  < E{Zk(t)} 


i 

. 


Proof:  Suppose  we  have  a process  Z3(t)=Z(t,u)3,£,r)  where  the  Z3  - process 
has  two  independent  input  sources,  both  with  Compound  Poisson  input  and 
intensities  X1(t)=A(t)-M1(t)  and  A^O-M^t).  Let  N^t)  be  the  number  of 
inputs  in  (0,t]  from  the  source  with  intensity  A.(t),  i=l,2.  If,  for  the 
Z3  - process,  we  know  that  N1(t)=0  then 

P(Z3(t)  = 0 | Nj(t)  = 0}  > P{Z3(t)  = 0} 
and  E{Z3(t)  | Nj(t)  = 0}  < E{Z^(t) } 

But  if  the  Z3  - process  only  has  input  from  its  second  source,  then 
since  X2(t)  = Mj(t) 

P{Z3(t)  = 0 | Nj(t)  = 0}  = P{Z1(t)  = 0} 
and  E{Z3(t)  | Nj(t)  = 0}  = E{Z*(t)} 

Note  that  the  overall  input  rate  of  the  Z3  - process  is  (X(t)-M^(t) )+ 

M j ( t)  = X(t),  which  implies  that 

P{Z3(t)  = 0}  = P{Z(t)  = 0} 
and  E{Z3(t) } = E{Zk( t) } 

Thus  we  have 

P{Z(t)  = 0}  < P{Zj(t)  = 0} 
and  E{Zk(t) } > E{Zk(t) > 

Similarly  we  can  obtain  the  other  inequalities  by  introducing  a pro- 
cess Z^(t)=Z(t,u)4,^,r)  where  the  Z^  - process  has  two  independent  input 


sources,  both  with  Compound  Poisson  input  and  intensities  X ^ ( t)=M^( t) -A ( t) 
and  A,,(t)=A(t).  Then 

P{Z(t)=0}=P(Z4(t)=0}>  P(Z4(t)=0  | Nj ( t ) =0 } = P{Z2(t)=0} 

E{Zk(t)}  = E{Zk(t)}  iE{Zk(t)  | Nj ( t ) =0 } = E{Zk(t)}. 

This  second  lemma  is  a basic  limiting  result  which  will  be  used  often. 
Lemma  3.2.2:  For  any  vpeL ^ 

A)  If  f ( x ) is  a function  such  that  |f|<k  and  f(x)-*-L  as  x-*»  then 

t 00 

lim  £ i^(u)f(t-u)du  = L / i^(u)du 

t-HJc  0 0 

B)  If  f ( x ) is  a periodic  function,  period  u>,  and  f(kuj+x)-*f  (t)  as 

ki+T  oo  ^ 

k-H»,  0<x<io,  then  lim  f ^(u)f(kaj+T-u)du  = ■(  ij'(u)f  (x-u)du 
k-x»  0 0 

Proof : To  prove  part  (A)  we  note  that  since  |f|<k  we  have 

(vp(u)f ( t-u)  |<  K|»|)(u)|,  and  K^»( u ) eL ^ . Also,  vp(u)f(t-u) -►Lvp ( u ) as  t-**. 

Thus  by  dominated  convergence 

t 

lim  / ijj(u)f(t-u)du  = L / (u)du. 

t-*oo 

Part  (B)  follows  in  a similar  fashion.  In  this  case  we  note  that 
kl!V ( u) f ( k£+t -u )=*( u) f ( t -u) . Thus  by  dominated  convergence 

, • ku)+T  oo 

k-»«’  o ^u)f(k£3+T"U)du  = £ *(u)f(x-u)du. 

Theorem  3.2:  Let  Z(t)  be  the  content  of  a Type  A storage  process  with 
nonhomogeneous  Compound  Poisson  input  which  has  intensity 
X(t)=A0( t ) ) ,i  >0  where  i|>(t)^0,  (|>»L,  and  «t>(t)  *0  as  t>-. 


If  the 
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corresponding  homogeneous  storage  process  with  intensity  XQ  has  a limit- 
ing distribution  of  store  content,  then 

lim  n(t,A(.),W)  = nKAj 

t-K»  0 

Proof:  By  Lemma  3.2.1,  since  X <A(t)  we  know 

o— 

(18)  n(t,x0,Ho)>  n(t,x(.),w0) 

Taking  upper  limits  of  both  sides  of  (18)  yields 

n(°°,A  ) > lim  n(t,A(- ),  (w  ) 
u t-*»  0 

where  n(°°,A  ) = lim  n(t,A  ,W  ) exists  by  assumption. 

° t-*»  0 0 

We  can  use  (16)  to  obtain  the  other  inequality.  By  Lenina  2.1.1  we 
have  P{Z(t)=0|Z(0)=0}  >_  P{Z ( t ) =0 ( Z (0 ) >0 } and  hence  n(t-u^,A0,W^)  <_ 
n(t-uk,A0,D)  where  D(x)  is  the  degenerate  distribution  function  D(0)=1. 
Then  we  have  the  inequality 


n(t,A(-),W0)<e 


-eA  *(t)  .t 

{n(t,A0,W0)+eA0  fQ  4>(u)H(t-u,Ao,D)du 


t u,  u 


k k c uk  2 

+ k=2  e Xo  o o”’o  • • - 4>(u^)n(t-u^,X0  ,D)du  ^...du^> 


-cA  *(t)  t 

= e ° (n(t,A0,W0)+  f eA04»(u)n(t-u,Art,D)du 


+ k=2ekxo  i 4>(uk)n(t-uk»X0*D)  o”K’o  ^(u1)...<J>(uk_1)du1...duk) 


uk  ,u2 


The  summation  can  be  simplified  using  the  identity 


Xk  / 1 


(20)  /•*•£  A4>(x0)  — ♦(xk_1)d>^  — dxk-l  = f i k<^(u)du5k 

R1 


This  identity  can  be  easily  proved  by  induction.  The  case  k=l  i 
Assume  (20)  is  true  for  k=n.  For  notational  convenience,  let 


xk  rxl 


♦k(,k)  * l— l »(x0)...*(xk.1)dx0...dxk_1  thm 


x„,.  xn  X, 
n+1  r n f 1 


4ntl(lVl’=  o o'"o 


■ £*ntl*(x„)['Vxn)]ndxn 

n! 


since  d<J>j(xn)  = <J>(xn)  we  have 
dx~ 


n! 


(n+l) ! 


obvious. 


as  claimed. 

Using  (20)  we  obtain 


-eX  4>(t 

n(t,x(-).w0)>e  0 (n(t,x0,w0)+GX0  4»(u)n(t-u,xo,D)du 


♦ kI2  A*  l ♦lu)^1(u)Ht-u,>0.D)du). 


Since  <|>elj  we  have  that 


t "Ik 

lim  f <t>(u)<tk-1  (u)du  = lim  [ 4>(u)du  < 
t~  0 TFTT!  t-~  Lo— n — 


i.e.  eLj.  In  addition,  as  n(t,*0,Wo)  is  bounded  and  approaches  a 
limit,  we  can  use  Lemma  3.2.2  to  show 


n(t-u,Xo,D)du  = 


n(°»,x0)  l 


i •x  1 

4>(u)<I>  (u)di 

IFTt! 


n(°°,xo)^  («). 


This  yields 


-eX  <&(<») 


ljm  n(t,x(-).w0)>e  0 {n(~,x0)+ex0n(«>,x0)<5>(«>) 


+ kf2ckXQn("°,xn)»k(<”)) 


-eX  eX  ♦(<») 

0 nKXje  0 


= n(»,A0). 

From  (19)  we  have  lim  n(t,  (’)»Wo)  < nfc,^)  and  from  (21)  we  have 
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lim  Jl(t,A( • ) >W  ) < n(°°,A_),  so  we  have 
t-x»  oo 

n(=°.x ) < lim  n(t,A(-),wj  < lTm  n(t,A(-),w  ) < n(~,A) 

° t-x»  ° t-*co  0 0 


which  gives  us  our  conclusion. 


Note  that  we  have  not  assumed  here  that  e<l.  Thus  Theorem  3.2  can 
apply  to  any  function  a( t ) which  can  be  written  in  the  form  a+b(t)  where 
a>0,b(t)>0  and  beLj,  since  we  can  write  A(t)=a(l+  ^ b(t)).  Or,  suppose 
that  d> ( t)->L  as  t-*»,  L>0.  If  <J>(t)>L  and  <j>(t)=<j>(t)-L  then  we  can  write 
A(t)=  A(l+e<j>(t))  where  A=A0(1+gL)  and  e=e/(l+eL).  If  in  addition  $eL, 
then  Theorem  3.2  applies  and  n(t,A(  • ) ,WQ)->-n(«>,A0(l+eL))  as  t-*». 


3.4  Periodic  Intensity 

The  situation  in  the  preceding  section  is  the  exception  rather  than 
the  rule.  In  most  cases  we  will  not  be  able  to  get  such  simple  answers. 
It  is  impractical  to  attempt  to  evaluate  all  of  the  terms  in  the  expan- 
sion of  IT,  or  in  fact  many  at  all.  We  can  perhaps  hope  to  calculate  a 

• L 

few  terms,  leaving  a remainder  term  of  order  0(e  ) as  r>0.  To  find 
the  probability  the  store  is  empty  at  t,  up  to  an  error  of  order  of 

Is 

magnitude  0(e  ) we  need  to  consider  at  most  k-1  inputs  of  Type  II  occur- 


i 

a 

i > 

4 • 

- 

.-  a-,.  -*- 


ring  in  (0 , t] . To  see  this,  recall  that 
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P{Z(t)=0}=  njjQ  P{Z(t)=0  and  N2(t)«n} 
k-1 

= I P{Z(t)=0  and  N?(t)=n} 
n=o  c 


+ P{Z(t)=0  | N2(t)>k}  P{N2(t)>k} 


and 


oo  -eA  <J>(t)  n . 

P{N2(t)>k>  = n^e  0 [cA0*(t)]  = 0(ek). 

fn 

We  will  only  look  at  the  first  two  terms  of  the  e-expansion,  which 
will  give  us  the  probability  the  store  is  empty  at  t,  up  to  an  error  of 

p 

0(e  );  i.e.  we  look  at 

(22)  n(t,X(.),W0)=n(t,A0,W0)(l-eX0<t(t))+eX0  ^(u )n( t-u ,A0,W,(u) )du+R2 

where  R2=0(e^). 

This  approximation,  while  appearing  simple  enough,  is  not  easy  to 
evaluate.  The  difficulty  lies  in  the  fact  that  is  a function  of  u. 
Since  we  are  interested  in  the  limiting  distribution  of  n,  this  problem 
can  be  alleviated.  As  we  have  seen  in  Chapter  2,  under  many  circumstances 
a limiting  distribution  of  store  content  will  exist,  independent  of  the 
initial  content  distribution.  Thus  an  advantageous  choice  of  initial 
distribution  could  simplify  things  considerably.  We  will  see  that  very 
often  a good  candidate  for  this  is  the  limiting  value  of  W(x,t,X0,WQ) 
which  we  will  denote  W(x,XQ).  If  this  is  taken  as  the  initial  distribu- 
tion, then  W(x,u1,XQ|u1)=W(x,Xo)*  B(x).  Note  that  this  is  not  dependent 


74 


on  Uj.  For  this  reason,  we  will  henceforth  write  Wj(uj)=Wj  in  this 
situation.  This  will  prove  helpful,  as  we  will  soon  see. 

In  order  to  be  able  to  take  advantage  of  the  limit  theorems  of 
Chapter  2,  we  assume  here  that  we  have  periodic  intensity,  with 
d> ( t+k ) =0 ( t ) and  <j>(t)  is  Riemann  integrable.  Then  $ can  be  expressed 
as  a Fourier  Series 


oo 

(23)  <l>(u)=a  + E (b  cos2i'mu+c  sin2[miu) 

' ^ o m=-«»  m m 


- E 3 Q 2irimu 
= L a e 


It  will  be  convenient  to  assume  that  £ 4>(u)du=0;  i.e.  ao=0.  If 

not,  then  we  can  define  a new  function  $(u)=<f>(u)-^<j>(x)dx.  In  this  case 

1 1 

A(u)=A(l+e$(u))  where  A=A0(l+e  £ 4>(x)dx)  and  e=e/(l+e  JQ  <|>(x)dx).  It 
will  also  simplify  matters  to  assume  <t>  can  be  expressed  as  a finite  series. 
We  can  substitute  (23)  into  (22)  to  get  the  approximation 


(24)  n(t,A(-),W)=n(t,A0,W)(l-e<J>(t)) 


Z 1 

+ m=-Me*o  o ame 


2TTimujj(  t_u  ,aq  ,Wj  )du+Rr 


Making  the  substitution  n(t,A0,W^)=n(t,A0,Wj)-n(«>,A0)  and  noting  that 
n(t,Ao,W)=n(®,Ao)  we  can  rewrite  (24)  as 


M . t 

II(t,A(-),W)=n(-,Art)+J  mcA  e27nmtam  f e' 
o m=-M  o mo 


2 iri  mu 


l(u  ,Aq  ,W j )du+R? 
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If  the  storage  process  satisfies  appropriate  conditions,  such  as 

those  of  Theorems  2. 2-2. 4 then  we  should  be  able  to  find 

n (t  ,A(  • ) )=1 1m  Il(n+r  ,A(  • ) ,W) , 0<t<1.  We  need  to  evaluate  terms  of  the 
r>*» 

form 


(25)  lim  fe’*""'  fi(u,A  ,W.)du. 

n-+on  U 0 1 

2"1mT 

We  claim  that  this  term  equals  e n (27iim,Ao>W1). 
First  we  note  that  if  fieLj  then  since 

I e-E’imu  fl(u,X0,W1)|<|S(u,X0,W1)|, 

by  Dominated  Convergence 


lim  /n+e~2TT^mu 

n-xxi  0 


n(u,A0,W^)du  = / e271imu  n(u,Ao,Wj)du 


n (2Tim,A0,W1). 


To  show  that  this  integral  is  meaningful,  we  look  at 


11m  3°(i  + le.i.K.)  - 

n-xx.  n o 1 o 

-( l/n+ifl)x 

We  Have  |e  fi(x,A0,W1)  |<|ll(x,A(),W1)  | , 

gence  again. 


n(x,Ao,W1)dx. 

so  by  Dominated  Conver- 


lim  n°(i  + le.A  ,W  ) = / lim  e 

n-*°° 


- ( l/n+i 3) x 


n(x,A0,w1)dx 
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r -iex  „ 

= o e n(x,X0,W1)dx 

= n°(i3) 

Having  determined  the  limit  in  (25),  we  can  now  write 

(26)  TI(t ,X(  • ) ) = H(»,X0)  + mLMeX0am  e2'ni”  n°(2^im,Xo,W1)  + R2 

Of  course,  the  above  limiting  formula  is  contingent  on  ITeLj.  We 
prove  this  in  the  following  theorem. 

Theorem  3.3:  Consider  a Type  A storage  process  with  Compound  Poisson  in- 
put and  release  rule  r(x)  such  that  there  exists  a level  c where 

r(x)>a>A3j  for  x>c 

r(x)<M  for  x<c 

— c — 

for  some  constants  a and  Mc>  Let  n(t,A,W)  = n(t,A,W)  - n(~,\).  If  the 
first  two  moments  of  the  input  size  distribution  (pj,  p^)  are  finite, 
then  IlcLj. 

Proof:  Let  cc  be  the  event:  ec  occurs  at  t if  Z(t)=c  and  Z(t-)>c. 

Define  the  following: 

q(t,u)  = probability  the  store  is  empty  at  t,  given  i occurred  at 

u,  and  r has  not  yet  recurred, 
c 

q(t)  = probability  the  store  is  empty  at  t and  t has  never  occurred. 


= distribution  of  the  time  to  the  first  occurrence  of  e 

c 

= distribution  of  the  time  between  successive  occurrences  of 


Note  that  since  ec  is  regenerative,  q(t,u)  is  a function  of  t-u.  We  will 


write  q(t,u)=q(t-u). 

We  have  the  following  equation: 


n(t,A,W)  = q(t)  + f q(t-u)dG(u)  + 
t t 

{ q(t-u)d[G*  F(u) ] + f q(t-u)d[G*  F*  F(u)]  + ... 


= q(t)  + £ q(t-u)d  H(u) 

where  H(u)=G(u)  + G*F(u)  + G*F*F(u)  + Notice  that  H(u)  satis- 

u 

fies  the  renewal  equation  H(u)=G(u)  + £ H(u-x)dF(x)  and  thus  is  a 
renewal  function. 

We  need  to  determine  n(«°,A  ).  We  have  shown  previously  (in  the  proof 
of  Theorem  2.5)  that  ec  is  certain,  so  q(t)<P{F-c  has  not  occurred  by 
t }-*0  as  t-*°°.  If  q ( t ) can  be  shown  to  be  of  bounded  variation  on  every 
finite  interval,  and  qeLj,  then  by  the  Key  Renewal  Theorem  we  will  have 


lim  l q(t-u)dH(u)  = - f q(u)du 
0 0 

where  = f udF(u).  We  have  q(u)<  dF(u)  = l-F(u)  and 

^ q(u)du  < ^ (1-F(u))du  = p^<™,  by  Lemma  2.5.1.  Thus  qcLj. 

To  prove  bounded  variation  we  refer  back  to  the  proof  of  Theorem  2.5. 
There  we  saw  that  for  the  class  of  inteivals  [0,a],  we  can  use  Lemma 





78 


2 of  Smith  (1955)  to  show  that  the  function  i|^(t)=4>A(t)(l-F(t))  is  of 

bounded  variation  in  every  finite  interval.  Let  A={0}.  We  have  AetV" 

Note  that  if  we  suppose  occurs  at  tQ  and  its  next  occurrence  is  at 

t +t.  then 
o 1 

<f>A(t)  = P{Z(tQ+t)  = 0 | Z(0),to,t1>t}  = q(t) 

is  the  function  we  are  interested  in. 

The  function  l-F(t)  is  monotone,  so  (l-F(t))"*  is  also  monotone. 
Then  (l-F(t))"*  is  of  bounded  variation  on  every  bounded  interval. 
This  implies  that  (t)/(l-F(t) ) is  also  of  bounded  variation,  since 
the  product  of  two  functions  of  bounded  variation  is  again  a function 
of  bounded  variation.  But  vp^(t)/(  1-F(t))=d>^(t)=q(t) . Thus  we  have 
shown  that  the  conditions  of  the  Key  Renewal  Theorem  are  met. 

With  this  we  have  that 


(27)  n(t,A  .W)  = q(t)  + l q(t-u)dH(u)  - l q(u)du 


= p(t)  + l q(t-u)d(H(x)  - -7  ) - -7  { q(u)du. 

U ^1  ^1  ^ 

OO 

— f r°°  - f 00 

Now  q(u)  < ' dG(v)  = l-G(u),  so  f q(u)du  = T (l-G(u))du  = Vj , say, 
but  Vj  £ i.e.  qcL^.  Looking  at  the  third  term  of  (27)  we  have 


f L q(u)dudt  < { (1-F(u))dudt 

0 t 0 t 


= 0 0 (1_F(u))dtdu 
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= £ u(l-F(u)du 
= \ fQ  u2dF(u) 

= \ y2  < * 

by  Lemma  2.5.1. 

For  the  second  term  in  (27)  to  be  in  L j,  we  must  have 
(2S)  I \(  q(t-x)d(H(s)  - i ) | dt  < - 

v u M j 

We  have  the  inequality 

oo  t 

f0  I f0  q(t-x)d(H(x)  - J ) | dt< 

® t 

f 1 q(t-x)  | d ( H ( x ) - £ ) I dt. 

Mj 

By  Fubinis  Theorem  this  can  be  written 

00  oo 

' { q(t-x)dt  I d(H(x)  - J )|. 

OO  oo 

The  inner  integral  is  f q(t-u)dt=  fQ  q(v)dv<  oo  since  qeLj. 

We  know  that  H(x)  is  a renewal  function  of  a delayed  renewal  process, 
where  the  distribution  function  of  the  time  until  the  initial  event  is 
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G(x).  Let  H (x)  be  the  renewal  function  of  the  ordinary  renewal  process. 


Then  H(x)=G(x)+G*  Hq(x).  Let  D(x)=HQ(x)  - — we  have 


G.0(x)  = [G.Ho(„)-i  ]♦[*  - I /G(y)dy] 


By  Theorem  8 of  Smith  (1954),  if  Uj  < 00  and  , then 


q | dHQ(x) - - | < °°  as  long  as  F(x)  has  an  absolutely  continuous  compo- 


nent. We  would  like  to  write  F(x)=L(x)+A(x)  where  A(x)  is  an  absolutely 
continuous  function.  Let  X be  the  time  between  successive  occurrences 


of  cc>  and  let  A(x)=P{X  < x and  no  input  in  (0,t*]J  where  t*  is  the  num- 


ber which  makes  the  following  hold: 


+*  - r dx 
1 o Tft) 


Thus  t*  is  the  time  the  store  will  first  empty  if  we  start  at  level  c 
and  there  are  no  inputs  in  (0,t*]. 


We  have  A(x)=0,  x < t*.  Let  Kc(x)  be  the  distribution  function  of 


the  time  from  the  first  input  to  an  empty  store  until  the  content  next 
crosses  c from  above.  Then  for  x > t* 


x-t* 


A(x)  = e~Xt*f  xe"Xu  Kc(x-t*-u)du. 


Let  Aj(x)  = l-e”Xx  and  A2(x)  = e~Xt  «c(x).  Then  Mx)  = Ai  * A;>(x-t*) 


Since  A^  is  absolutely  continuous,  so  is  A.  Thus  f | d(l(x)  | < 
Since  G(x)  <1  / ( dG  * D(x)  | < f \ dD(x)  | 


* 1 


I f wo  can  show 
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fn  I <*(£■-  h l s(y)dy)  I < 


o 1 'Uj  o 


00 

then  we  will  have  { I d(G  * H (x)  - - )L  . But 

o'  0 u j I < ® 

f° | d(-  - - /*G(y)dy)|  = f°\  I - dx 

0 1 \ 0 UVy#Ujr/  1 0 1 pj  JJj 


oo 

= r £ (1-G(x))dx  < °° 


since 


/ I <*(H(x)  - £)|  < / | dG(x)  | + / | d(G  *H  (x)  - £ )| 

00 

we  thus  have  { ( d(H(x)  - — ) | < °°.  All  this  implies  that  (28)  holds, 
o M j 

t 

so  ^ q(t-x)d(H(x)  - ^ ) el_^.  We  therefore  have  shown  that  fl(x,X,W)eL^. 


The  formula  in  (26)  is  not  as  formidable  as  it  seems.  We  look  at  a 
simple  case.  Suppose  <J>(t)=cos2iTt.  This  can  be  written 

*(t)  = 2-e"2nit  + ^e27T1t.  This  makes  (26)  look  like 

r.  -2m  T , 2m' T 

(29)  H(T,A(-))-w(»,A0)+eA0[^  e n°( -2m  .X^)-^  e n°( 2tt1  .X^Wj)  j+R, 


We  have  fi°( 2">Ti  ,Xo,W1)=n°(2iTi  ♦ X0.W1)-n(°°,^)/2iri  and  thus  (29)  becomes 


RVH'JK  KU'iJ* 
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n(x,\(.))-nKxj4€A  L’2lTiT 


o^xo  le  IW-^VWjHik* o} 


-Trrf 


+ e2lTlT  n0(2iri,X0tM1)Hl(..X0)  e2l,ix 


2th 


+ R„ 


Let  n (2TTi,Xo,W1)  = x+iy.  Then  we  can  simplify  the  above  equation  to 


(30)  n(x,X(-))  = n(~,X  ) +eX  {xcos2irr-ysin27rr-n(<»,X  )sin2irr}+  R 

0 — 2 


The  values  of  x and  y will  depend  on  the  specific  storage  process  we 
are  interested  in.  Both  the  input  size  distribution  and  the  release  rule 
are  essential  ingredients.  We  will  explore  this  more  extensively  in 
Chapters  4 and  5.  There  we  treat  the  release  rules  r(x)*c ,r(x)*a+bx  and 
r(x)-cx,  x>0.  The  expected  content  will  also  be  examined  in  these  cases. 


CHAPTER  4 


LIMITING  RESULTS  FOR  THE  STANDARD  STORE 

The  standard  release  rule  in  storage  theory  is  r(x)=>c,  x>0.  For 

such  a process,  the  content  of  the  store  is  released  at  a constant  rate 

as  long  as  the  store  is  not  empty.  In  this  chapter  we  will  investigate 

at  some  length  this  storage  process.  We  focus  on  the  particular  case  of 

nonhomogeneous  Compound  Poisson  input  with  periodic  intensity 

A(t)=X  (l+e<j>(t) ) where  /*<{>(  u)du=0.  Without  loss  of  generality  we  assume 
o o 

that  c=l.  A simple  change  of  time  scale  can  transform  a system  with 
c^l  to  one  with  c=l. 

4.1  The  Limiting  Probability  of  Emptiness 

A store  with  r(x)=l  is  a Type  A storage  process,  since  R(x,0)  = 

J*  dy  = x < «>  for  all  x < oo.  Thus  Theorem  2.3  can  be  used  to  establish 
o 

the  existence  of  a quasi-limiting  distribution  of  store  content.  We 

need  only  assume  that  0.X  = 3,  /^(ujdu  < 1. 

10  i o 

In  order  to  determine  limiting  values  of  either  the  probability  of 
an  empty  store  or  the  mean  store  content  we  will  need  to  look  at  the 
homogeneous  store.  The  approximations  derived  in  Chapter  3 require  values 
for  n(»,xo,W) ,y(°°,X0,W)  and  n°(z,X0,W).  The  first  two  quantities  are 
standard  results  in  storage  theory:  n(°°,X0,W)  = 1-AQ31  and  y(°°,X0,W)  = 
X032/2(l-X061).  (See  Appendix.) 

A value  for  the  Laplace  Transform  Jl°(z,xo,W)  can  be  determined 
through  the  aid  of  the  differential  equation  for  the  Laplace-Stieltjes 


Transform  of  store  content.  From  Lemma  3.1  this  equation  is 


(1)  3M(s,t,A0,W)  + A0(l-B*(s))M(s,t,X0,W)  = sM(s,t,X0,W)-sIT(t,X0,W) 

5T 

Taking  Laplace  Transforms,  with  respect  to  t,  of  both  sides  of  equation 

(1)  yields 

zM°(s,z,\o,W)-M(s,0,\o,W)  + X0(l-B*(s))M°(s,z,X0,W) 

= sM°(s,z,A0,W)-^n°(z,A0,W) 
or 

(2)  M°(s,z,A0,W)  = M(s,0,Ao,W)-sJ7°(z,Ao,W) 

Z-S  + X0(l-B*(s)r 

Since  M(s,t,\o,W)  is  a Laplace-Stieltjes  Transform  it  is  analytic  for 
R(s)>0.  The  Laplace  Transform  of  M is  also  analytic  in  the  right  half 
plane,  because  an  analytic  function  of  an  analytic  function  is  analytic. 
Thus  any  zero  in  R(s)>0  of  the  denominator  in  (2)  must  also  be  a zero  of 
the  numerator. 

Let  ip( s ) =z-s  + X ( 1-B*(s ) ) . We  will  use  Rouches  Theorem  to  show 
that  tjj(s)  has  a zero  in  R(s)>0.  Fix  z,  R(z)M).  Define  an  increasing 
sequence  of  real  numbers  >k,  0 < >k  < 1 such  that  >k  t 1.  For  each  inte- 
ger k,  let  rk  be  the  contour  {u:  |u-a|  = x -r  k } where  a=A  +z  and  >k  is 
chosen  such  that  1 j^X  ( • We  a^so  let  f(s)=z-s+AQ  and 
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Clearly,  both  f and  gk  are  analytic  in  and  on  r^.  On  r^, 

|f(s)|  = A0-ek  and  |gk(s)|  < AQyk  | B*(s)  | < AQYk.  By  assumption 

|f(s)|  = Vek>  Vk  - I 9k(s)  I 

Hence  by  Rouches  Theorem,  the  number  of  zeros  of  f in  rk  is  equal  to  the 
number  of  zeroes  of  f+gk  in  rk.  There  is  only  one  zero  of  f in  rk 
(at  s=XQ+z).  Thus  there  is  one  zero  of  f+gk  in  rk-  Call  the  zero  sk» 
Let  r be  the  contour  (u:  ( u-a  | = XQ}.  Certainly  sk  is  inside  r. 
Since  {sk>  is  a bounded  infinite  sequence  there  is  an  s*  and  a subse- 
quence sk  such  that  sk  -*■  s*  as  n-*°°.  The  point  s*  is  in  or  on  r.  We 
n *n 

want  to  show  that  ^(s*)=0.  Let  g(x)  = -A  B*(s)  = lim  g. (x).  We  have 

0 k-o°  K 

>Ksk  ) = (f  + gk  - gk  + g)(sR  ) 
n n n n 

* (f  + 9k  )(sk  ) + (g  - gk  )(sk  ) 
n n n n 

= (g  - gk  )(sk  ) 
n n 

Now 

(9  - 9k  )(sk  ) - -U-fk  > X B*(sk  ) 
n n n n 

and 

l(g  - gk  )(sk  )|  £ (1-Yk  ) A0  -*■  0 as  n 


-*  00. 
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Thus 

lim  4»(s.  ) = 0 = \|>(s*) ; 
n-x»  n 

i.e.  there  is  a zero  of  ij/(s)  in  or  on  r. 

If  R(s*)>0  then  M°(s*,z,A0,W)  <°°,  so  we  must  have 
M(s*,0,Aq,W)-s*  ji°(z, Aq,W)  = 0.  If  R(s*)  = 0 then  we  must  check  to  see 
if  M°(s*,Z,Ao,W)  < ».  But  I M(s*,t,A0,W)  I <1,  so  M°(s*,z,X0,W.)  < ~ < 
and  again  we  must  have  M(s*,0,Ao,W)— s*  n°(z,A0,W)  = 0.  Thus  we  have  the 
relationship: 

(3)  n°(z,A0,W)  = M(s*,0,Ao,W) 

s* 

The  actual  value  of  n°  will  depend  on  the  initial  distribution  of 
content.  Looking  at  equation  3.26  we  see  that  the  distribution  of 

interest  is  W^x)  = W(x)  * B(x)  where  W is  the  limiting  distribution  of 

content  of  the  homogeneous  store,  intensity  XQ,  and  B is  the  distribution 
of  input  size.  For  every  t,  M(s,t,A0,W)  = M(s ,»,A0,W) . Thus 
M(s,0,Xo,W1)  = M(s,»,X0,W)  * B*(s) . Smith  (1953)  showed  that 

(4)  M(s,~,A0,W)  = s(l-X0B1) 

i^yT-B*(s)) 

Recalling  that  z - s*  + Xq(1-B*(s*))  = 0,  we  can  now  write 
n°(z,X0,W1)  = (1-XpB^  B*(s*) 


z 
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The  value  of  B*(s*)  can  be  computed  in  two  ways.  The  most  obvious 
way  is  to  find  s*,  the  zero  of  vp( s ) . Another,  more  illuminating  way,  is 
to  consider  the  following.  Let  G(x,aq)  be  the  distribution  of  the  busy 
period.  It  is  well  known  that  the  relationship 

(5)  G*(z,X0)  = B*(z+A0G*(z,A0)) 

holds.  Let  u = z+a  -A  G*(z,A_).  Then  G*(z,a  ) = U ~ Z ~ Xo  . 

U O U U ■» 

" Ao 

Substituting  into  (5)  gives 


or 

(6)  -AqB*(u)  +Aq+z-u=0 
A zero  of  (6)  is  u = s*.  Thus 

(7)  B*(s*)  = G*(z,A0)  . 

Returning  to  the  original  problem,  we  can  now  write 

n°(z,A0,Wi)  = (1  - A^j)  g*(z,a0) 

z 

= (1  - X081)G°(z,A0) 

and  we  thus  have 

n(t,AQ,W1)  = (1  - A061)G(z,A0)  . 

Finally,  we  can  write  down  an  equation  for  the  limiting  value  of  the 
probability  of  emptiness  as 
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(8)  lim  n(n+rtX(.).W)s(l-X  Bj)+eX  (1-X  Bj)  S_M  ame27nmT[G*(2Trim)-l]+0(c2) 
n-*°°  2nim 


where  4>(t)  = JM  a e27Timt  . 
Y m=-M  m 


4.2  The  Mean  Store  Content 

An  approximation  to  the  mean  store  content  can  be  obtained  by  using 
the  storage  equation: 

N(t)  t 

(9)  Z(t)=Z(0)  + jZ0  U.  - f r(Z(u))du. 

Taking  expectations  of  (9)  we  have 

t 

y(t,X(-),W)  = y(0,x(-).W)  + %EN(t )-/  ( l-n(u,x(  • ) ,W)  )du 

t t 

= u(0,X(.),W)  + ^ X(u)du— t + f n(u,x(*).w)du 


*y(0,x(-),W)  + BjXjjt  + BjX^O-t  + / n(u,x(-),W)du 


We  can  substitute  for  n(u»X( - ) »W)  into  equation  (10)  to  obtain 

t 

p(t,X(-),W)  = u(0,X(-).W)  + (B1X0-l)t  + eX0B1*( t)  + f n(°°,X0)du 

+ /»  »„  /te2’1mu/V2’imx  fi(x,x  ,W,)dxdu  ♦ 0(e2) 

o m=-M  m o o oi 

Writing  n(°°,X0)  = 1-X  and  y(0,X(-)»W)  = y(°°,X0)  and  noting  that 
ll  e2"(ra(u-x)5(x>A0.W1)dxclu  = /tn(x,»0.H1)  e2l,™(t-x>l  dx  , 
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we  have 


(11)  u(t,A(*) ,W)=y(oo,Ao)+eX081$(t)+eX0  mI_M2irim/  IUx.Xq.Wj)^ e27r1m^"x^|dx+Ofe2 


Write  ^(tJ-ujI.^ci^e  . Note  ^at  following  identities  hold: 


fam/2irim 

k^O 

{ M 

|jn=-M  (V2™) 

pr 

II 

O 

M t 


r. 

L i 


w(t,X(*)  »VJ)=y(«,»X0)+eX081$(t)+eX0  mL^o  e^^^nU.X^WjJdx 

ntfO 


+eXoao  o n(x,Xo,W1)du+0(e^). 


As  in  Chapter  3,  we  have  for  m^O, 

lim  / +Te27rim(n+T-x)n(x,X0,W1)dx=e2lTlT  no^irlm.A  ,1^) 

n-+« 

=e2lTimT(i-xoe1)(G*(z,X0)-l)/2. 


When  m=0 
n+x 

lim  { 

n-x» 


MTT^  oo^ 

lim  f n(x,XQ,W1)dx=^  n(x,XQ,W1)dx 


=o  (i-Vi^s(x,V“1)dx 

= -(l-A^bj 

where  is  the  expected  length  of  the  busy  period.  Thus  we  can  finally 
write 
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limu(n+T,X(*)  ,W)  » p (°°» AQ) -«-cXoe ^<l>( t ) -eX0<  1-Xq3j ) Y 

IV® 

M „ . 9 

(12)  + eA0(l-X031)mi_MattleZlTimT(G*(2Ti1m,X0)-l)  + 0 (cZ) 

m/o  27rl"m 

We  shall  henceforth  denote  the  approximations  to  the  limits  in  (8)  and 
(12)  by  fi( x ,X(  * ) ) and  u(t  ,X(  • ) ) • 

It  does  not  follow  immediately  that  the  limit  of  u(n+x,A( . ) ,W)  is 
equal  to  the  expectation  of  the  quasi-limiting  random  variable.  It  is 
still  necessary  to  prove  a condition  such  as  uniform  integrability.  We 
do  this  by  proving  uniform  integrability  for  a sequence  of  random  vari- 
ables which  bounds  our  store  content  from  above.  According  to  Lemma 
3.2.1,  for  two  processes  with  the  same  release  rule  and  Compound  Poisson 
input 


uk(t,A(-).W)  < pk(t,X0(l+cM),W) 

A.L. 

where  4>(t)  < M and  uk(t,A( • ) .W)  denotes  the  ku  moment  of  store  content 

at  time  t of  a store  with  intensity  \(*)  and  initial  distribution  W.  If 

we  can  show  that  supu0(n+T ,A  ( 1+rM)  ,W)<  •»>  then  supuJn+r,A( • ) ,W)< , 

n t o n c 

which  implies  uniform  integrability. 

In  standard  queueing  and  storage  theory  it  has  been  shown  that  as 
long  as  XQ(  1+cM)0j<1  and  32<a,then*  If  A^X^  l+eM) , 

(13)  lim  u, (n+x.X.  ' » A.B9/2( l-A.B, ) 

n«  1 1 12  11 

where  1-AjBj)  Is  the  expectation  of  the  limiting  value  of  store 

content.  Since  store  content  is  nonnegative,  (13)  implies  uniform 
integrability  for  the  homogeneous  process.  Thus 


supu,(n+T  ,A(  • ) ,W)  < supvj9(n*T  ,\.  ,W)<  - 
n i - n c » 
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and  we  have  uniform  integrability  of  the  nonhomoqeneous  process.  This 
implies  the  anticipated  conclusion. 


4.3  Special  Case:  $(t)  ■ cos2n t 

Let  us  now  look  at  a specific  example  in  more  detail.  Suppose 
<f>(x)=cos2nx  ■ ^ e 2n’x  ♦ ^ etn’x.  Equation  (8)  becomes 


(14) 


1 e-2l,1,6*(-2*1.X0)-  ^e-2niT 

~-27T 


27T 


+ 1 e2l,1lG*(2iri.X0)-  J 

Equation  (14)  has  an  alternative  form  without  any  complex  numbers.  Let 
G*(2iri  ,AQ)=x+iy.  Then  G*(-2ir1,^))=x-1y  and  (14)  can  be  written 

(15)  II ( t ,A(.)  )>=(  1-Aq3 j )+c^0(  l-A082){xsin2irT+ycos2TTT-sin2iTT} 

—%r~ 

To  determine  p(x,A(*))  in  this  case,  first  note  that 

./  v 1 [ 2iHt  - 2tt i t ] . Since  a =0  we  do  not  need  to  know  the 

*(t)  Le  - e J o 

expected  length  of  the  busy  period  here.  The  approximation  to  the  mean 
in  this  particular  case  is: 


p(T.x(*))axoe2 


*n^vp 


*rXoP1s1fl2jt*t»o(l-Xoe1)  { jjTe~2’lHG^(-;n.) 


-sr 
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(16) 


= X 

2\ 


Xo32  “ c*0(l-*0fl1)(x“D<:os27n'  + eA0  r 2n01+(  1-X()61  Jylsi 
V_T  « ) — 7 — — 7~  L J 

1 Vr  4/  4/ 


n2irT 


4.4  A Queueing  Theory  Example 

The  actual  value  of  n and  p depend  on  the  distribution  B(x). 

Suppose  that  we  have  the  most  basic  queueing  situation-negative  exponen- 
tial service  time  with  average  service  time  1/p.  It  is  well  known  that 
the  Laplace-Stieltjes  Transform  of  the  busy  period  of  an  M/M/1  queue  is 


G*(z,Aq)=  p + Xq  + z 


2A_ 


(p+Ao+z)2-4Aop 


and  the  mean  waiting  time  in  steady  state  is  Ao/p(p-AQ). 

Figures  4. 1-4.3  present  graphs  of  the  approximation  (15)  of  the 
limiting  probability  of  emptiness  in  an  M(t)/M/1  queue  with 
A(t)  = A0(l+ecos2iTt) . We  take  p=10  for  all  three  figures,  while 
Aq  = 1.0,  5.0,  9.5  in  Figures  4. 1-4.3  respectively.  In  each  figure,  c 
takes  on  the  values  .1,  .25,  .5  and  .9. 

The  values  of  n(«>, A ( • ) ) fluctuate  symmetrically  about 
n(«,A0)  = l-A0/p,  for  each  e,  with  the  amount  of  fluctuation  dependent 
on  the  size  of  e.  This  was  to  be  expected  from  the  form  of  equation 
(15).  It  is  also  intuitively  reasonable  that  the  probability  should 
vary  about  the  emptiness  probability  of  the  M/M/1  queue  which  uses  the 
average  intensity  per  period  as  the  arrival  intensity.  In  fact,  the 
M(t)/M/1  queue  with  A(t)  = Aq( l+e<}>(t) ) i*  usually  approximated  by  the 
M/M/1  queue  with  A( t)=AQ. 


— 
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It  is  interesting  to  note  how  the  extreme  values  of  n( t ,A( • ) ) differ 
from  n(~,Ao).  This  gives  an  indication  of  how  much  the  process  is 
smoothed  by  using  a homogeneous  intensity  rather  than  a nonhomogeneous 
one.  The  extrema  of  n(x,A{.))  can  be  determined  by  differentiating  (15) 
with  respect  to  t: 

an(TtA(-))  = eA  (l-A  61)2TT{xcos2iTT-ysin2TTT-cos2TTT}. 

Setting  the  derivative  equal  to  zero,  we  find  that  the  extrema  occur  at 
points  tq  in  the  period  such  that 
tan2iTTo=(x-l)/y. 

Note  that  the  extrema  occur  at  the  same  point  in  each  period,  irrespective 
of  the  value  of  e,  which  only  effects  the  value  at  that  point.  The  ex- 
trema for  each  graph  are  shown  in  Table  4.1.  If  we  consider  the  percent- 
age deviation  of  the  extrema  from  n(°°,A0),  we  see  that  the  percentages 

increase  as  A„  increases.  These  values  are  in  Table  4.1. 
o 

It  is  also  of  interest  to  note  the  time  from  the  maximum  (minimum) 
value  of  A ( t ) to  the  corresponding  minimum  (maximum)  value  of  n(x,A(*)). 

We  shall  call  this  the  response  time,  as  it  measures  how  long  it  takes 
for  the  queue  to  fully  respond  to  extremes  in  the  arrival  pattern.  Values 
of  response  time  are  found  in  Table  4.1.  Note  that  as  AQ  increases,  the 
response  time  increases. 

Table  4.1 


xo 

Extreme 

Values 

% Deviation 
from  n(°°,A0) 

Response 

Time 

1.0 

. 90+ . 08e 

9e 

.0995 

5.0 

.50+. 23c 

46e 

.1355 

9.5 

.05+.045e 

90c 

.1655 

PROBABILITY 


PROBABILITY 
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FIGURE  4.3 

Probability  of  an  Empty  Queue 
r(x)=l 

Xq=9.5  p=10.0 
X(t)=XQ(  l+ecos2irt) 
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Figures  4. 4-4. 6 are  graphs  of  the  mean  waiting  time  of  the 
M(t)/M/1  queue,  X(t)=Xo(l+ecos2iTt).  Again,  we  let  vj*10  and  XQ*1.0,  5.0, 
9.5  respectively  for  the  three  figures,  and  e takes  on  the  values  .1, 
.25,  .5,  .9.  The  graphs  of  the  mean  waiting  time  are  similar  to  those 
of  the  emptiness  probability.  The  values  of  m(t,X(*))  fluctuate  sym- 
metrically about  p(°°,Xo).  The  extrema  occur  at  points  tq  such  that 

tan2nT()=2T\61+(  1-XQ$1  )y 
(l-X0Bj)(l-x) 

Table  4.2  gives  the  extrema,  percentage  deviation  of  the  extrema  from 
u(°°tX0),  and  the  response  time  for  each  graph. 

Table  4.2 


*0 

Extreme 

Values 

% Deviation 
from  u{°°,X0) 

Response 

Time 

1.0 

.011+.009e 

82e 

.1071 

5.0 

.10  +.06e 

60e 

.1768 

9.5 

1.90  +. 14e 

7 e 

.2434 

We  see  that  while  response  time  Increases  as  XQ  increases,  the  percentage 
deviation  of  the  extrema  from  y(°°,X0)  decreases,  in  contrast  to  what 
happened  with  n(x ,X( • ) ) . 

4.5  General  Service  Distribution-An  Approximation 

As  mentioned  earlier  in  the  chapter,  to  determine  explicitly  the 
value  of  n(x,X(* ))  and  v»(t,X(  - ))  for  a store,  we  need  to  either  know  the 
busy  period  distribution  or  be  able  to  find  the  zero  of  the  equation 
z+X(l-B*(u))-u=0,  R ( u ) >0 . This  may  not  be  an  easy  task  for  an  arbitrary 
distribution.  An  approximation  to  many  general  distributions  which  often 


MEAN 


FIGURE  4.6 


Mean  Waiting  Time 
r(x)=l 


Aq=9.5  ua10.0 
A(t)=Xo(l+c°s2'nt) 
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is  simpler  to  deal  with  is  the  Weighted-sum  Erlang  distribution  which 
has  density 

oo 

(17)  B'(t)=b(t)=k|  cke‘ytp(yt)k_1;  0<ck  < 1,  ^ ck=l 

( k"l ) ! 

With  proper  choice  of  {cm>  and  y it  can  approximate  a wide  range  of  dis- 
tributions. 

The  Weighted-sum  Erlang  distribution  can  be  seen  to  be  related  to 
the  more  basic  negative  exponential  distribution  (which  it  contains  as  a 
special  case,  when  c^=l),  deriving  from  the  relationship  its  advantage 
as  an  approximation.  This  is  easier  to  understand  in  the  context  of 
queueing  theory.  Suppose  we  let  A(t)cmdt+o(dt)  be  the  probability  that 
during  (t,t+dt)  a customer  requiring  m phases  of  service  arrives  at  the 
queue.  If  the  distribution  of  each  phase  of  service  is  negative  exponen- 
tial with  mean  1/y  then  B ( t ) as  in  (17)  is  the  resulting  (total)  service 
time  distribution  for  each  customer. 

It  is  also  interesting  to  note  that  this  same  approach  can  be  used 
to  model  a bulk  queue.  In  this  case  x(t)cmdt+o(dt)  is  to  be  interpreted 
as  the  probability  that  during  (t,t+dt)  m customers  will  arrive  at  the 
queue,  each  with  negative  exponential  service  time. 

The  Laplace-Stieltjes  Transform  of  B(t)  is  easy  to  find.  We  have 


(18) 
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where  C(x)  is  the  probability  generating  function  of  {cm}.  Since  we 
know  that  G*(z,A0)  = B*(z+A0-A0G*(z»X0) ) we  can  solve  for  G*(z,A0)  once 
we  determine  the  probability  generating  function  of  {cm} . 

4.6  A Special  Case:  cm=(l-a)am-1 

Suppose  cm  = (l-a)am~*.  m>0.  Then 

C(x)  = x(l-a) 

1-ax 


Thus 


G*(Z,A0)  = C ( z+A  - A G* ( z , A ) +p } 


i 


- (1'a)w„ ->0e«u,yn. 

1 - ay ^ 

«V»0G*(z.X0)^ 

= 

Z+A0  Aq6*(z ,AQ)+y-ay 

This  can  be  transformed  to  the  equation 

Ao[G*(z’Ao)]2-[z+Ao+y"au]G*(z’Ao)+w(1"a)=0 
Solving  this  for  G*  gives 

G*(z,A0)=z+A0+y-ay  + sJ(z+A0+y-ay  )2-4AQy(l-a) 
^ 

The  correct  solution  is  the  one  which  makes  G*(0,Ao)=l. 
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We  still  need  to  determine  n(<»,X0)  and  p(~,Xo)  to  be  able  to  find 
n(T  .X(  • ) ) and  0 (t  .X ( • ) ) . We  know  that  n(°»,X0)  ■ l-XQ0j  and 


01  “ “a?  B*(s)  I s*0 


Js  C(s+Ti)  U-0 


In  the  particular  case  c *(l-a)a  m>0. 


k i 

s+p 


(s+p-oip)^  | s=0 


p(i-a) 


"'-•V  ■ 1 ■ „tt4t 


order  to  determine  p(*»,X0)  we  need  to  know  the  value  of  0^,  where 


ds^B*(S)l, 


'!b*c(»Ms-o  • 


104 


Again,  for  our  example  we  have 


h dV  ' L‘.(l-a)-, 

i (s+u-ap)‘ 


- MtsU 

(s+p-mp)' 


rrr  t? 

p ( i ; 


Xo  2-, 
t<("v\  ) = p (l-o) 

0 ?rpypn=5T) 


u(  l-«)£p(  l-it)-\Qj 


4.7  Another  Example 

In  Figures  4.7-4.10  we  have  graphs  of  the  probability  of  emptiness 
and  mean  waiting  time  of  the  Mx(t)/M/1  queue  where  \(t)  = \o(  l+i  cos2'it) 
and  cm  = ^ We  take  p=10  for  all  figures,  and  \ * 1,  5.  We  do 

not  have  any  graphs  for  \ o * 9.5,  since  for  this  value  of  XQ, 

3 

Vi  - 9.5  x ^ > 1,  and  no  quasi-limiting  value  of  waiting  time  exists. 
As  before,  in  each  figure  t takes  on  the  values  .1.  .25,  .5,  .9. 

These  graphs  are  similar  to  the  others  we  have  seen.  We  see  that 
the  response  time  increases  as  X increases,  and  the  percent  deviation 
from  the  average  increases  for  the  probability  of  «'mpt.iness,  but  dot  reuse 
for  the  mean  waiting  time.  Tables  4.3  and  4.4  give  these  values. 
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FIGURE  4.10 

Mean  Waiting  Time 
r(x)=l 

v*(r~ 

Xq=5.0  m=10.0 
A(  t)=XQ(  l+eCOs2irt) 


CHAPTER  5 


LIMITING  RESULTS  FOR  SOME  STORES  WITH  GENERAL  RELEASE  RULES 

There  are  many  release  rules  that  are  of  interest  in  storage  theory. 
The  most  comnon,  r(x)=c,x>0  was  considered  in  the  previous  chapter.  The 
output  of  such  a store  is  not  a function  of  the  content  (except  that 
r(Z(t))=0if  Z(t)=0).  In  the  class  of  "general"  release  rules,  where 
output  is^  dependent  on  the  content  of  the  store,  the  release  rules 
r(x)=a+bx,x>0  and  r(x)=cx,x>0  are  the  most  popular  in  the  literature. 

Both  these  release  rules  will  be  discussed  in  this  chapter. 

5.1  The  Limiting  Probability  of  Emptiness;  r(x)=a+bx 

We  look  first  at  a storage  process  with  nonhomogeneous  Compound 
Poisson  input,  periodic  intensity  X(t),  and  linear  release  rate 
r(x)=a+bx  ; a>0,  b>0.  We  would  like  to  use  Theorem  2.2  to  show  that  a 
quasi -limiting  distribution  of  store  content  exists.  To  do  this,  we 
need  to  look  at  Y(u).  Recall  that 


Y(u)  = 


;u  dx 
o rTx) 


= f 


dx 


o a+bx  , in  the  present  case. 


The  assumptions  of  Theorem  2.2  are  met  as  long  as 


p 


1 
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■ 1 1 

‘ 

(1) 

£ X(u)du j 

E 

ln(l-4  U) 

L « 

Thus  we  can  conclude  that  a quasi-limiting  distribution  of  store  content 
exists  as  long  as  (1)  holds. 

1 

Suppose  X(t)*XQ(l+e<j>(t))  where  £ is  periodic  and  / 4'(u)du=0.  The 
formula  for  the  limiting  value  of  the  probability  of  emptiness  is  depend- 
ent on  the  existence  of  a limiting  distribution  of  the  content  of  the 
corresponding  homogeneous  store  with  intensity  X . From  Theorem  2.5 
this  limit  will  exist  if  there  exists  a level  c such  that  r(x)>XQBj  for 
xn:  and  r(x)vMc  for  x<c.  If  aNXQBj  then  r(x)^X0Bj  for  all  x>0  and  r(x) 

is  bounded  on  all  finite  intervals.  If  a^X  t>.  we  can  choose  c=h(\  B.-a). 

o 1 b o 1 

Then  r(x)=a+bx>X0^1  for  x>c  and  r(x)vXoBj=Mc. 

To  be  able  to  compute  n(r,\(’).W)  we  need  to  determine  n('",Xp»W) 
and  n°(2T’im.xo,W1)  for  m=+l,+2,  . . . ,+M  and  Wj (x)=W(x)*B(x)  where  W is 
the  limiting  distribution  of  store  content  for  the  homogeneous  store. 

An  examination  of  the  Laplace-Stieltjes  Transform,  M(s,t,\p,W),  will  be 
helpful.  From  Lenina  3.1  we  know  that  this  satisfies  the  following 
equation: 

(2)  3M(s,t)  4f \ (l-B*(s))-sal  M(s ,t)+bs3M(s ,t)=-asn( t) 

3t  L J 3s 

where  we  write  M(s ,t,xo.W)=M(s ,t)  and  n(t,.\  .W)»Il(t). 

Let  M(s)  be  the  Laplace-Stieltjes  Transfonn  of  the  limiting  value 
of  store  content.  M(s)  satisfies  (2)  with  aM(s)/^t*0.  Then  in  the  limit 
(tw)  we  have  the  differential  equation 


I 
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d|sI+xo[izEl|ija]S(s) 


The  solution  to  this  is 


(3)  M(s)  = e9(s) 


[C'F  "Mo 


e-0(x)dx 


where  0(s)=^  s — ^ ^ dx,  and  c is  some  constant.  When  s=0, 

M(0)=1.  This  implies  that  c=l.  As  s approaches  infinity,  M(s)  must 
approach  JI(°°),  which  is  finite.  But  exp(0(s)  )-*»  as  s-*»,  so  we  must  have 


s 

lim  1 - | n(°°)  / e'0^dx  =0 
s-*»  0 


(4)  n(«)  = 


a r 


Substituting  (4)  into  (3)  gives 


(5)  M(s)  - ee(s)  C e'e(x) 

r .-e(x) 


We  also  note  that 

lim  -y(")4nH4  + lim  X°(1"-^-(l-  N(s) 


An  application  of  l'Hopital's  rule  yields 


(6)  uH  = an(»)  - a + X^ 
b 
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It  still  remains  to  determine  n^z.A^).  We  go  back  to  equation 
(2),  where  we  take  the  Laplace  Transform  with  respect  to  t on  both  sides 
of  the  equation,  to  transform  the  partial  differential  equation  into  a 
differential  equation  in  M°(s,z).  This  gives: 


zM°(s.z)-M(s,0)+[xo(l-B*(s))-sa]  M°(s ,z)+bs3M°(s ,z)  = -san°(z) 


3s 


Reorganizing  the  equation  into  a more  familiar  form  we  have: 


(7)  3M°(s ,z)  + 
3s 


z+Aq( 1-B*(s) )-sa 


L 


bs 


M°(s,z)  = M(s,0)-saH°(z) 
bs 


Solving  (7)  we  have 
.s 


M°(s,z)=exp{-§  { ^ + 6(s)Hc -/  exp  I { ^ -0(x) 


M(x,0)  a o,  , 
bx  "b  ,l  (zl 


dx} 


=s-z/bee(s){c./s>2/be-8(x) 
0 


' ^ f "°(2)>X> 


As  s-w,  s **>  also.  Since  M°(s  ,z)-*-n°(z)<  We  must  have 


= /“  xz/b  e_0(x) 


M(x,0) 

bx 


I n°(z) 


dx. 


Thus 


(8)  M°(s,2)=s-z/b  ee(s)<fx2/'’  e'0'*’ 


. i n°(z) 
bx  b Jl  [Z> 


dx} 


Evidently,  lim  M°(s,z)=z  ^ °°.  But  lim  s"z^b  e0^=».  Then  we  must  have 

S-w  S-*™ 

that  the  terms  inside  the  braces  in  (8)  go  to  zero.  This  yields  the 
following  expression  for  n°(z): 
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(9)  n°(z)  = C xz/b~l  e"0^  M(x,o)dx 


a r xz/b  e-0(x> 

o 


The  value  of  n°(z,A0,W^)  is  dependent  on  the  initial  distribution 
Wj(x)=W(x)*B(x) . Thus  we  need  to  substitute  M(u)B*(s)  for  M(x,o)= 
M(x,o,Ao,W^)  in  equation  (9),  yielding: 


n°(i.*0.V-  C xZ/b'1  e'e(x)  e9(x)  If  e'9{u)|lu  1 f e'9(u)<lu  B*(x)dx 

a/Ve<*>  dx 


f°  x2/b-1  B*(x)  /“  e‘6*u'  du  dx 


r ^ .-aw, 


at  J x"'"  e v 'dx II ^ e “'"'dx 


5.2  A Special  Case:  Computational  Aspects  of  rr 


We  now  have  all  the  tools  to  evaluate  n(T,A(*))-  All  that  is  needed 
is  to  designate  values  for  a,  b,  B(x)  and  <f>( t ) . We  will  specialize  here 
to  the  case  <ji(t)=cos2irt,  a=l,  b=AQ  and  B(x)=g1e“x^l.  We  assume  that 
l+b$j<e  to  ensure  that  a limit  exists.  This  is  so,  since 

AqEY(U)=  E^ln(l+bU)j  <_  -g-  ln(l+bgj)<l  by  Jensen's  inequality.  We  note 

that  e(x)=£  — ^ ln(l+$^x).  Then 


n(-,A0) 


V0(x)dx 
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(ii)  = (l+Bjb)"1 


Also 


(12)  n°(z,Ao,W1)  = f u^h'^l+BjU)"1  /°  e~x/b(l+B1x)dxdu 

/V/b  e'^d+SjUjdu  f e'^a+e^du 


Computationally,  there  are  still  a few  considerations,  as  the 
evaluation  of  JI°(z)  at  z=  + 27ri  can  present  difficulty.  We  determine 
the  inner  integral  of  the  numerator  in  (12)  as 


f ( l+31x)e'x/bdx=( l+B1u)be"u/b+b2B1e"u/b 
Let  N stand  for  the  numerator  in  (12);  then 


N = b f uz/b_1  e"u/bdu  + b2e1  f e^/V^d+B^r^u 
= Nj  + N2,  say. 

is  easily  evaluated  as  proportional  to  a gamma  function,  using  the 
simple  change  of  variable  u=by;  indeed 


N,  = e'y  y2^1  dy 

(13)  > bz/b+1  r(z/b) 
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The  integral  r(+  ¥")  cannot  readily  be  evaluated,  but  r(+  ¥”  + 1) 
can,  and  we  then  can  use  the  identity  r(x+l)=xr(x) . We  have 

r(±¥+l)  =fe-V2li/b  dy 

= /°°  e-y  e±  2lTi/b  1o9  y dy 

= / e"y  cos ^ log  yjdy  + i f e"ysin[  log  y ]dy 

The  two  functions  on  the  right  are  easily  computable. 

The  evaluation  of  N2  cannot  be  done  straightforwardly.  We  write 

n2  = b2fh  f e’u/b  uZ/b_1(l+e1u)"1  du 

= b2e1  f e'u/b  uz/b"1|^(l+elu)‘1-ljdu+b2e1  /e"u/b  uz/b_1  du 
(14)  = -b20  2 /°°e"u/b  u^d+e.u)"1  du  + bg.N, 

10  1 lx 

We  note  that  the  first  integral  on  the  right  can  be  written,  for 
z=+2iri , as 


I 

i 

j 


i 

i 

I 


- b231¥e_u/b  ir^^^d+PjU)"1  du  = 

(15)  -b2312/0°e_u/bcos|^~  1 ogu J ( 1+pju) “1du+b2e12^  e"u/bsin|^^-loguj(l+e1u)"1du. 


If  we  also  write  D for  the  denominator  in  (12)  we  have 


- - t m*  ■ — — 


r 


— 1 — - 
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D = [^°°(l+01x)e'x/b  xz/bdx][/Vu/b(l+PlU)du] 

' idfr  [fe'X/b  «2/b^+».  fe-/b  «2/b+1  *>] 

■n^?[b2/btlr(l*‘>tV2/btH^>] 


z/b+2 


(16)  = nFjgr  1+ei2+lhb 


r(§  + 1} 


With  the  above  representations  of  Np  and  D we  can  compute  n°(z)  and 
thus  ff(x  ,X( * ) ) - 

5.3  Limiting  Value  of  the  Mean  Store  Content;  r(x)=a+bx 

The  mean  of  the  storage  process  is  also  of  particular  importance. 

To  derive  the  approximation  we  must  determine  the  mean  for  the  homogeneous 
store.  Clearly,  from  (7) 


lim  9M°(s ,z)  = lt°(z)  + j~  lim  M°(s,z)  + li 

S-x»  9S  s-*°° 


lm 

S-KO 


M(s,o)-(z+A0(l-B*(s)))M(s,z) 
Fs  ~ 


We  note: 


(1)  lim  M°(s,z)  - I ; 

S-KO 

(2)  By  1 'Hopital 's  rule 

lim  M(s,o)-(z+X0(l-B*(s)))M(s,z)  = -y(0)+Zp°(z)-AoB1/z 


S-+O0 


bs 


This  yields 


U°(z)(z+b)  = p(0)  + (yyaj/z  + ail°(z) 


L iJ 


which  implies 


y(t)  = pfOje^+U^-a^V^du+a  ^n(u)e‘b(t"u)du 

(17)  = fy(0)-X 0&1  + a \ e”bt+  X^-a  + a ^(uje'^'^du. 

\ ~V  b)  ~b~ 

Having  performed  this  preliminary  calculation  of  y(t),  we  substitute 
(17)  into  equation  3.17  to  obtain  the  desired  approximation.  Note  that 
u(0,Xo,W1)  = u(«>,Xo)  + Bj.  Thus 

t ? 

u(t,X(-),W)=u(”,Xo)(l-eXo<l'(t))+eXo  / <J)(t-u)u(u,Xo,W1)du+0(e  ) 

rx  B,-a  t 

=M(oo,X0)-eX0y(<»,X0)4>(t)+EX0|j— g £ <j>(u)du 

+ (y(^0)+3r^[p--)  ^4*(t-u)e'budu 

+ a £ 4>( t-u)  f n(x,XQ,W1)  e_b^u_x^dxdu] +0(e2) 


■ ^an(»,X0)-a+X0B1)-cX0  an(~,X0)4>(t) 


[1  t . v t u -b(u-x)  o 

5Tl(00,Xo)+ei j/  <^(u)e"DU"u;du+cXoa^  <fr(t-u)/  n(x,XQ,VJ1)e  dxdu+0(e  ) 


5.4  A Special  Case:  4>(u)-cos2ttu 

In  the  case  <J>(u)*cos2itu,  4>(t)=si-n^-2--t-^  and 

£ cos2iTue"b^"u^du=e"bb/ ebu[bcos2iTU-»-2Trs1n2iTu] 

bZ+47T2 
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(19) 


= bcos2TTt+2nsin2nt-be 

XT T— ~ 


-bt 


bu+4TT 


To  determine  u(t,X(*))  we  need  to  study,  from  (18), 

£ cos2ir(t-u)  £ n(x,XQ,W1)  e_t)^u”x^dxdu  = J(t),  say.  Since 

|cos2ir(t-u)e'b^u"x\[(x,Xo,W1)  |<  e'b^u"x^ , integrable  in  [o,t]x[o,u],  we 
can  interchange  integrals  to  get 
J ( t ) = f n(x,XQ,W1)  / cos2it(t-u)e'b^u"x^dudx 


= f n(x,X  ,W.)  <bcos2n(t-x)+2TTSin2TT(t-x)-be"b^b~xA 

v O 1 i n 

b^+4TT 


dx 


by  (19).  Making  the  substitution  ft(x)=n(x)-n(<»)  we  have 


l 


t } t 

1 


-b(t-x) 


J(t)=n(«,XQ)/b/  cos2TTudu+2iT^  sin2irudu>  -b^  n(x,Xo,W1)e'uv  ^ A/  dx 


2 2 

b‘+4iT 


b 1 fi(x,X  ,W,  )cos2ir( t-x)dx+  2ir  .rfi(x,X_  ,W,  )sin2ir(t-x)dx. 

T — O JT”  0 0 1 —35 ^ « 0 1 


b2+4it2 


b^  0 


To  determine  p(t,X(*))  we  need  to  determine  1 im  J(n+r).  As  we  saw 

n-*» 

in  the  previous  chapter 


n+t 


lim  lfi(x,Xo,W1)cos2Tr(t-x)dx=^e2lTiT^2TTi  .XQ,W1) 


1 -2iritir 


n-*«> 

and 


n+T 


lim  £ ft(x,Xo,Wj)sin2n(t-x)dx=^- e27lilll0(2ni  .Xjj.Wj)-^  e ^n1 ' n°(-2ni  ,XQ,W^) 


1 -2itiT, 


n-«o 
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From  Lemma  3.2.2,  since  n(t,AQ,Wj)  is  bounded  and  converges  to  a 


_ K V 

limit,  and  e eLj,  we  have  that 
n+t 


lim  f n(x,A0,W1)e'b(n+T-x)dx=n(«-,Ao)/Vbxdx=ifl(~,Ao). 


n-**> 

Thus 


p ( x , A ( • ))=g-{an(°°,A0)-a+Aog^)-»eXoan(«’,Ao)sin2TTT 

2?rb 


+ eA, 


an(<»,An)+e, 

1 bcos2TTT+2Trsin2iTT  I 

b 0 \ 

l b2+4TT2  ■* 

+ eA0an(«>,A0)  bsin2Ttx  - cos2ttt+1  - eA  an(»,A  ) 

~9 5 2tt  — 5 3-^ 

b^+4ir  b2+4ir 


+ eVe2lTlT  + If) 


b 2tt, 


b^+4? 


ey  e'27IlT  fi0(-21ri,A0,W1)  (|  - |f) 


b 2tt  \ 


,T7T 
b 4ir 


= an(“,A0)-a+A031  + j bcos2TTT+2usin2TTTj 

6 J 


(20)  + eAoa 

2(b2+4u2) 


(b-2iri )e2lTl T n°(2iri,A0,W1)+(b+21Ti)e'2lTlTfi0(-2TTi,A0,W1) 


Equation  (20)  can  alternatively  be  written  without  using  complex  num- 
bers. First  we  note  that  n°(z)=n°(z)-n(«>)/z.  If  we  let 
n0(2ni ,A0,W^)*x+iy,  then  after  some  algebraic  manipulation  we  can  write 
(20)  as 


I - ~ . 


cos 2ht 


P(t.X(-))  = all(<»,A0)-a+X0ei  + cAQ 


b +4n 


bp^+a(bx+2ity)+all('*',\o)  i coszm 


(21)  + cAq  I 2n0j-a(by-2iTx)-a  bn(“,Ao)  | sin2m 

b2+4Ti2  ] ^ 


As  before,  we  need  to  establish  that  1 ini  p(n+T,A(*)»W)  is  the  expec- 

n-*» 

tation  of  the  limiting  value  of  store  content.  We  look  at  the  homogeneous 
store  with  r(x)=a+bx  and  A1=A0( l+e)>A(t) . We  know  from  (6)  that 

p(oo»A1)=^(an('»tA1)-a+A101)  as  long  as  a+bx>A0Bj.  We  also  see  from  (17) 
that 

lim  p(t,A.,W)  = A.p.-a  + alU^.A.) 

JF~  «> 

= vi(oo*x1) 

Since  we  know  that  a limiting  distribution  of  store  content  exists, 
lim  p(t,A( • ) »W)=p(«>,A  ),  and  all  random  variables  are  nonnegative,  we 

t->co  0 

know  that  the  store  content  is  uniformly  integrable.  In  addition,  since 
by  Lemma  3.2.1  we  have  u(t,A( •) ,W)<p(t,A, ,W) , the  content  of  the  nonhomo- 
geneous  process  must  be  uniformly  integrable,  which  gives  us  what  we  seek- 
that  for  each  t,p(t ,A( • ) ) is  the  approximation  to  the  expectation  of  the 
quasi -limiting  value  of  store  content. 

5.5  Some  Graphs 


Figures  5. 1-5.6  are  examples  of  graphs  of  the  probability  of  an  empty 
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store  and  mean  store  content  of  the  nonhomogeneous  store  with  r(x)=l+bx, 
A(t)=Xo(l+ecos2irt)  and  XQ3b.  We  have  B^=.l  for  all  figures  and  Ao=l,  5, 
9.5.  In  each  figure,  e takes  on  the  four  values  .1,  .25,  .5  and  .9. 

The  probability  of  an  empty  store  when  r(x)=l+bx  has  form  similar 
to  the  earlier  examples.  The  mean  store  content,  while  similar  in  form, 
varies  considerably  more  than  other  stores.  In  fact,  at  some  values  of 
tande,  u(t,A(*))  is  less  than  0-  This  is  possible,  since  u(t,A(»))  is 
only  an  approximation,  and  many  of  the  remaining  terms  in  the  series 
expansion  of  u(t,A(-))  may  be  positive.  Of  course,  if  u(t,A(-))<0  we 
take  it  to  be  zero. 

Tables  5.1  and  5.2  give  the  extreme  values,  percentage  deviation  of 
the  extrema  from  the  average,  and  response  time  for  n(x ,X( * ) ) and 
u(t,A(-))  respectively.  While  not  so  indicated,  the  values  of  p are 
assumed  to  be  nonnegative.  It  should  be  noted  that  the  difference  in 
the  figures  reflect  not  only  the  differences  in  XQ,  but  also  the  release 
rules,  since  we  have  assumed  that  AQ=b.  Thus  we  cannot  directly  compare 
either  the  graphs  or  the  entries  in  the  table. 


Table  5. 1 


*o 

Extreme 

Values 

% Deviation 
from  n(“,AQ) 

Response 

Time 

1 

.91+. 08c 

' ' 

9c 

.0775 

5 

.67+. 27c 

41e 

.0777 

9.5 

.51+. 37c 

72c 

.0659 

Table  5.2 
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1 


X 

0 

Extreme 

Values 

% Deviation 
from  v* (<V,*XQ) 

Response 

Time 

1 

.009+. 007c 

78c 

.0399 

5 

.033+. 044c 

133c 

.4497 

9.5 

.049+. 038c 

78c 

.4249 

5.6  A Type  B Process 

Up  to  now  we  have  spent  a good  deal  of  time  examining  the  emptiness 
of  a store.  But  there  are  some  stores  which  cannot  empty  in  a finite 
amount  of  time.  These  are  Type  B processes.  The  content  of  a Type  B 
process  can  get  arbitrarily  close  to  zero,  but  Tl( t ,X ( • ) .W)  =0.  What  is  of 
interest  here  is  the  mean  store  content.  A typical  release  rule  for  a 
Type  B process  is  r(x)=cx.  For  this  particular  store  we  will  see  that 
we  can  get  more  than  just  an  approximation  to  the  mean  content. 

We  begin  with  the  basic  storage  equation 

t 

Z(t)«Z(0)  + A(t)  - c / Z(u)du 

where  A(t)  is  the  total  input  in  (0,t].  Assuming  the  input  is  Compound 
Poisson,  X(t)=Xo(l+e4'(t)),  we  take  expectations  of  both  sides  of  (22)  to 


p(t,X(*),W)=n(0,X(.),W)+B1N(t)-c  £ u(x,X( • ) ,W)dx 

t 

(23)  m(0,X(-),W)+B1(Xo+fXo4'(t))-C  / }i(x,X(  • ) ,W)dx 

We  take  Laplace-Transforms  with  respect  to  t,  of  both  sides  of  (23)  to 
get 

u°(s  ,x<  • ) ,w)m(  0,x(- ) ,w)+B1x0+cB1x0<f0(s)-cu°(s  ,x(  • ) ,w) 


s 


s 


ILITY 
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FIGURE  5.2 

Probability  of  an  Empty  Store 

r(x)=l+5x 

A =5.0  3, 10 
o l 

A(t)=Xo(l+ecos2n) 


MEAN 
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FIGURE  5.4 

Mean  Store  Content 
r(x)=l+x 
XQ=1.0  6^.10 

X(t)=AQ(  l+ecos27rt) 


TftU 


MEAN 
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u(t,X(.),W)  - v<(o,X(-),W)e'ct+P1Xo  / e'^du+eX^j  ^(u)e'c(t"u)du 
(24)  = [n(o,X(-).W)-l]e"ct+01Xo  +cXo  <j>(u)ecudu 


We  take  note  here  of  a few  Important  points.  First,  equation  (24) 

t , , 

only  requires  that  £ <Mu)e’c'  ~u'du  exist  for  all  t.  Second,  we  have 
shown  that  the  entire  expansion  of  u In  c is  contained  in  the  first  two 
terms.  We  have  no  need  to  rely  on  an  approximation.  Third,  we  see  that 
the  mean  store  content  depends  on  the  input  size  distribution  only 
through  the  mean  of  the  distribution. 

We  now  restrict  ourselves  to  the  case  of  periodic  d>( t) . The  last 
term  in  (24)  becomes 


■V!  e‘cto  ♦(u)eC',du  ■ \ i ‘‘"’‘"e'"* 


-Ct  V . r*"  2niku  cu 


rX  ft  £ a 2n1kt  £ 

EVl  k=-M  ake  e L-  1 


e k*-M  ak 


OT+C 


If  we  go  to  the  limit  as  n goes  to  infinity  we  get 


(25)  u(x*X(*))“11m  p(n+T,X(*)»W)“B,X  ♦ rX  p.  a.e2l1lkl 

n-x-  -±-2.  0 1 K~"M  k (2itiV+cT 


We  would  like  to  be  able  to  say  that  w(t,X('))  does  Indeed  equal  the 
mean  of  the  quasi-limiting  distribution  of  store  content.  First,  we  show 


3 


that  a quasi-limiting  distribution  exists,  as  long  as  *03j<c,  by  virtue 
of  Theorem  2.4.  For 

Y(U’c)  = ^ ^ = ^ tlnU“1nc]* 

and 

ur  ■ X(u)<lu JeY ( U ,e) 

* x0  i [ElnU-lnc]. 

The  function  f(x)=lnx  is  concave,  so  by  Jensen's  inequality 
ElnU  < 1 n [eu]  < E U * 8j 

Thus,  for  any  e>0 

and  the  assumptions  of  Theorem  2.4  are  satisfied. 

Let  us  now  consider  the  homogeneous  store  with  X^=\0(l+rM),  where 
<j>(t)<M.  The  Laplace-Stieltjes  Transform  of  store  content  can  be  deter- 
mined from  the  partial  differential  equation 


(26)  3M(s ,t ,Xj  ,W)+Xj( 1-B*(s) )M(s  ,t ,A j ,W)=  - sc3M(s ,t  ,Xj ,W) 

ait  3s 

Let  M (s,X^)  denote  the  Laplace  Stieltjes  Transform  of  the  limiting  value 
of  store  content.  We  know  this  exists,  since  r(x)NAjflj  for  x^A^,  so 
Theorem  2.5  applies.  Then  M*(s,A.)  satisfies 


(27)  d M*(s .Aj)  = -A1(l-B*(s))M*(s,A1) 


To  find  u(<»,A1)  we  take  the  limit  of  as  s*0,  in  (27).  We  get 

u(°°»A  j)  = XBj/c. 

But  we  already  know  that 


y(t,Aj,W)  = ufO.Aj.W)  e"ct+  BjAj  (l-e'ct) 


11m  u(t,A. ,W) 
t-*o  1 


^(“.Aj)  = AjBj 


Since  the  store  content  and  limiting  value  of  store  content  In  the 
homogeneous  store  are  both  nonnegative  random  variables,  we  have  that  the 
content  Is  uniformly  integrable.  By  Lemma  3.2.2,  since  A(t)<A1  we  have 
that  u(t,A( • ) ,W)  < u(t,Aj ,W) , from  which  we  can  infer  that  the  content 
in  the  nonhomogeneous  store  is  uniformly  integrable.  With  this  we  can 
conclude  that  for  the  nonhomogeneous  store 


Urn  E 


rz(n+x)|  = E lim  Z(n+x) 

L J In-*30 


5.7  An  Example:  <fr(t)=cos2nt 

When  <j>(t)*cos2nt  equation  (25)  becomes 

u(x.X(-))  * 6,xn  + tU,  fi  e"2"'  + i e2’*T 

JjS  [2  EXST  2 c5®T  J 
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= B1A.*e*oe1  [c  cos2irx+2Tisin2TfT] 
c c^+4t? 

We  now  look  at  a specific  example,  the  store  with  r(x)  = | x, 

A(t)  = Xo(l+ecos2Ttt) . Figures  5. 7-5. 9 are  graphs  of  the  mean  content  for 

this  store.  We  let  8j=.l  and  AQ=1,  5,  9.5  for  the  three  figures,  while  e 

takes  on  the  values  .1,  .25,  .5  and  .9  in  each  figure. 

The  graphs  of  the  mean  store  content  are  symmetric  about  u(®,A0). 

The  extrema  occur  at  points  t such  that 

o 

tan2irr  =2tt  . 

0 T 

We  note  that  tq  does  not  depend  on  either  the  input  rate  or  input  size 
distribution.  Thus  the  response  time  is  only  a function  of  c.  The 
actual  value  of  the  extrema  points,  though,  do  depend  on  the  mean  input 
size  and  the  average  intensity  \Q. 

Table  5.3  gives  the  extrema,  percent  deviation  of  extrema  from 
y(“tX0)  and  the  response  time  for  each  graph.  The  response  time,  of 
course,  is  the  same  for  all  A0,  while  the  percent  deviations  of  the 
extrema  from  u(°°,Ao)  remain  virtually  constant.  Thus,  except  for  a scale 
factor,  for  a given  e the  graphs  are  essentially  identical. 


Table  5.3 


*0 

Extreme 

Values 

% Deviation 
from  u(un,AQ) 

Response 

Time 

1 

.3  +.02e 

6.6c 

.2416 

5 

1.5  +. le 

6.6  e 

.2416 

9.5 

2.85+. 2c 

7c 

.2416 

MHRN 
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FIGURE  5.9 

Mean  Store  Content 
r(x)=jx 

X =9.5  6,*. 10 
o 1 

X(t)=X0(l+ecos2"t) 


CHAPTER  6 


A MATRIX  APPROACH  TO  THE  NONHOMOGENEOUS  QUEUE 


6.1  A Step  Function  Approximation  to  x(t) 

In  this  chapter  we  depart  from  our  previous  perturbation  approach 
to  the  nonhomogeneous  store  and  look  at  an  alternative  method  for  dealing 
with  the  M(t)/M/1  queue  which  relies  heavily  on  the  computer.  The  method 
was  discussed  in  a paper  by  Leese  (1967).  Basically,  the  method  involves 
using  a step  function  approximation  to  the  arrival  intensity  to  generate 
transition  matrices  for  each  time  period  during  which  the  arrival  inten- 
sity is  constant.  More  specifically,  suppose  we  assume  that  the  nonnega- 
tive real  line  [o,°°)  is  divided  into  intervals  In=[tn_j,tn)  where  tQ=0, 

CO 

lk=[o,o°),  and  on  each  interval  the  arrival  intensity  is  approximated 

by  a constant,  i.e.  x(t)~Ak  for  tel^.  If  this  is  done,  then  on  each 

interval  1^  we  essentially  have  an  M/M/1  queue  with  arrival  intensity  X^. 

1/ 

We  can  compute  the  probability  transition  matrix  P (t),  which  has 
entries 

PN,n(t)  = p{^k^)  = n I = N}*  teIk 

k/ 

where  Q (t)  is  the  number  of  customers  in  the  system  at  time  tel^,  using 
the  well  known  formula  for  the  transient  probability: 


(t)  = 


-(A.+u)t 
e * 


(N-n)/2 

(2jy~  t) 


n+N+1 


{'-m 


? 

l=n+N+2 


(2JV l) 

{\i'\  (uv  *) 


Since  the  system  is  Markovian,  by  repeated  multiplication  of  the  Pk 

we  can  find  the  system  size  probabilities.  For  example,  suppose  te^ 

and  we  want  to  know  the  probability  there  are  n in  the  system  at  time  t 

if  we  started  at  time  0 with  N.  Looking  at  tj  we  see  that  there  might 

have  been  j customers  in  the  system  at  that  time  with  probability 

PN  j^l^*  The  Pr°babi  1 i ty  of  going  from  j to  n in  the  system  in  the 

? 

remaining  t-tj  time  is  ^t-tj).  This  gives 


1 2 

This  is  just  the  (N,n)  entry  in  the  matrix  product  P (t^P  (t-t^).  It 
follows  in  general  that  if  tel^  then  P{Q(t)=n  jQ(0)=N>  is  just  the 
(N,n)  entry  in  the  matrix  product. 


(2)  P^tj^ytj).  . .Pk'1(tk_1-tk_2)pk(t-tk_1) 


In  his  paper,  Leese  complained  that  even  with  a judicious  choice  of 
step  function  and  form  of  the  formula  to  compute  (1),  the  amount  of 
computation  required  would  be  formidable  and  perhaps  excessive.  Thus 
he  did  not  recommend  this  technique  as  a practical  approach  to  the  non- 
homogeneous  queue.  Since  Leese  wrote  his  paper  more  than  10  years  ago. 
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much  progress  has  been  made  in  computer  technology.  It  is  likely  that 
in  many  instances  the  amount  of  computation  required  would  now  be  con- 
sidered routine. 

In  the  particular  case  of  periodic  arrival  intensity,  this  matrix 

method  seems  well  suited,  even  perhaps  in  Leese's  time.  For  periodic 

X ( t ) the  step  function  need  only  be  fitted  to  the  first  period, 
n 

I=[o,Sj).  If  I = kU j I ^ and  X ( t) =X^  for  tel^  then  we  only  need  N matrices 
to  find  P(u>).  To  find  the  limiting  value  of  the  size  of  the  system,  we 
continue  squaring  P=P(w)  until  the  limit  is  reached.  This  should  not 
require  many  multiplications,  as  we  have: 


P(2u>)  = P • P = P 2 

(3)  P(4w)  = P(2S>)  • P(2G)  = P4 


P(2ku>)  = P(2k_1u))  • P(2k'1w)  = P?k 

L 

Thus  to  get  P(2  S)  required  k multiplications  once  P(u>)  is  obtained. 

6.2  Construction  of  the  Transition  Matrices 

Let  us  look  at  this  method  in  greater  detail.  We  shall  first  con- 

L 

sider  the  construction  of  the  basic  transition  matrix  P (t).  The  for- 
mula  for  the  entries  in  this  matrix  was  given  by  (1).  This  form  is  not 
well  suited  for  numerical  work,  as  it  is  not  clear  how  many  terms  of  the 


infinite  series  of  Bessel  functions  are  needed  for  any  particular  level 
of  accuracy.  Leese  suggests  rewriting  the  formula  as  a power  series  in 
X,  but  leaving  the  factor  This  gives  the  much  simpler  formula: 
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(4)  PKi  n(t)  = Z 
N,n'  ' m=o 


(U]V^  Fj  n(t) 


m! 


N,n 


where 


FN,n(t)  = Sm+N-n(t)  + Vn(t)'Tm-n-l(t)  * M(m+N+1>t) 


Tj(t)  " 


(ut)j 

j ! j >0 

0 


j>0 


(5) 


Sj(t)  . 


3-1 

z 


= uE.  Su(t)  = 1-  ^ Su(t) 


It  is  interesting  and  helpful  to  note  that  the  terms  in  (4)  can  be 


interpreted  probabilistically.  The  factor  (xt)me"Xt  is  the  probability 

m! 


of  exactly  m arrivals  in  [0,t],  while  the  factor  F™  (t)  is  the  probability 
of  there  being  n in  the  system  at  time  t,  given  that  there  were  N in  the 
system  at  time  0,  and  there  were  exactly  m arrivals  in  [0,t].  Since  F is 
a probability  we  can  make  use  of  the  fact  that  it  is  bounded  above  by  1 to 
determine  the  number  of  terms  we  need  to  use  in  (4).  We  have  that 


P ftl  - T e t rm  /.»  y (At)me 

PN,N(t)  — m=o  FN,n(t)  + m=k+l  ml 


m -At 
e 


The  error  term 


■ am 


. 


j- 
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■ 


z Ut)VXt  , , -At  z (At)m 
m=k+l  ml  m=o  m! 


can  easily  be  calculated,  and  an  appropriate  value  of  k chosen  for  any 
given  level  of  accuracy.  Another  nice  aspect  of  (4)  is  that,  if  we  can 
choose  a step  function  with  equal-sized  time  intervals  then  only  one  set 
of  F values  is  needed. 

The  transition  matrix  P is,  theoretically,  doubly  infinite.  Of 

course,  computationally  that  is  not  feasible.  Leese  did  not  deal  with 

this  particular  problem  in  the  method.  Naturally,  we  must  use  a finite 

matrix.  What  we  have  done  here  is  to  truncate  the  matrix,  making  it 

nxn.  To  retain  the  stochastic  nature  of  the  matrix,  we  replace  the  last 

entry  in  each  row  by  n-2  Thus  the  last  entry  is  in  a 

1 - j£o 


sense  the  probability  that,  starting  with  i customers,  the  queue  will 
have  n-1  or  more  customers  at  time  t.  Since  we  also  limit  the  number  of 
rows,  we  are  actually  restricting  the  system  to  be  no  bigger  than  n-1, 
and  hence  the  value  of  IT  will  be  somewhat  overestimated,  and  y under- 
estimated. 

At  what  size  should  P be  truncated?  This  is  a difficult  question 
to  answer.  Here  we  only  begin  to  explore  the  answer.  Since  we  are 
essentially  limiting  the  system  size  by  truncating,  we  should  look  at 
the  loss  involved.  For  the  homogeneous  queue  in  steady  state 
P(Q>n}=  pn  where  p=A8j  is  the  traffic  intensity.  If,  say,  p=.95  then 
p{Q>45>=. 099 ; i.e.,  if  we  use  a 45x45  transition  matrix,  allowing  at 
most  44  customers  in  the  system,  we  "turn  away"  or  ignore  1 in  10  custo- 
mers. Even  if  n=100,  P{Q>100}=.006.  It  seems  that  in  some  cases  the 


* 
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size  of  the  matrix  required  to  give  reasonable  results  will  surpass  the 
limits  of  many  computers.  Consider  further  that  to  model  a nonhomo- 
geneous  intensity  will  in  general  require  a number  of  these  matrices. 
The  following  table  shows  the  "loss"  involved  in  truncating  a queue  at 
size  n-1,  for  a few  values  of  p and  n.  It  can  be  used  as  an  aid  in 
determining  the  size  of  matrix  needed  to  be  within  a particular  level 
of  accuracy. 

Table  6.) 


\ n 

p\ 

10 

20 

40 

50 

75 

100 

140 

.95 

.599 

.358 

.129 

.077 

.021 

.066 

<.001 

.90 

.349 

.122 

.015 

.005 

<.001 

.75 

.056 

.003 

<.001 

.50 

.001 

<.001 

.30 

<.001 

6.3 

An  Example: 

A ( t ) a 

Step  Function 

It  is  interesting  to  note  an  additional  use  of  the  matrix  approach. 
For  an  M(t)/M/1  queue  where  A(t)  is  a step  function,  the  matrix  approach 
should  yield  values  of  n and  y which  are  very  nearly  exact.  If,  in 
addition,  A(t)  is  a periodic  function,  we  can  compare  the  values  of 
n(t,A(- ))  and  y(x,A(*))  with  the  approximations  n(x,A(*))  and  y(x,A(»)) 
presented  earlier. 

We  now  turn  our  attention  to  a specific  example.  Suppose  A(t)  is 
a step  function  on  each  cycle.  The  simplest  example  is  where  A(t)  takes 
on  only  two  values:  A ( t ) =A ^ for  tcfk.k+w^  and  A(t)=A2  for*  tr[k+Wj ,k+l). 
Here  we  will  set  A^=9,  A2=3,  y=10  and  Wj=.l.  This  could  be  a crude  model 
for  a situation  where  for  a short  amount  of  time  we  have  a very  high 

rate  (busy  or  peak  period)  while  the  rest  of  the  day  the  arrival 


arrival 


L 


rate  is  moderate.  We  find  from  Table  6.1  that  a 40x40  matrix  should  be 
sufficient  for  our  purpose.  Only  5 multiplications  of  P(l)  were  necessary 
to  reach  convergence.  To  find  the  mean  waiting  time,  we  can  use  the 
relationship 

(6)  p(t,X(-))  = E[Q(t) 1 

p 

where  Q(r)  is  the  quasi-limiting  value  of  system  size  at  time  t which  was 
shown  by  Luchak  (1956).  Note  that  n(x,X( • ) )=P[Q(t)=0]. 

We  will  compare  these  two  quantities,  H(t,A( - ) ) and  p(r,X(*))  with 
the  approximations  yielded  by  the  perturbation  approach.  We  need  to 
write  A(t)=AQ( l+ed>(t) ) , where  <f>(t)  is  periodic.  If  we  have  that 

A(t)  = JXj  k<t<K+wx  k=0,l  ,2,. . . 

|x2  k+Wj<t<l 

we  can  choose  XQ,  e and  <j>. , i = 1 , 2 such  that 

4>(  t ) = f<t>1  k<t<k+Wj  k=0,l,2,... 

|4>2  k+Wj£t<l 

and  A ( t ) = Aq(  l+ed>( t ) ) . 

We  would  like  to  be  able  to  write  <t>(t)  in  the  form 

<7>  ♦<*>  ■ j-M 

or,  in  the  more  traditional  form  of  the  Fourier  series 

00 

(8)  4>(x)  = cQ  + kT1(ckcos2Ttkx+dksin27rkx) 
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Term  by  term  integration  of  (8)  gives 


f |(x)dx  * cQ 


1 « 1 
^ |(x)cos2nnxdx  = f cos2nkxcos2unxdx 


w /* 

+ ^ d^  £ sin2nkxcos2nnxdx 


^ cos2irkxcos2Trnxdx  = f 0 k/n 


f cos2nkxsin2nnxdx  = 0. 


Thus  we  have  that  for  each  n 


2 / |(x)cos2nnx  = c . 
o n 


Similarly,  we  can  show 


2 ^ <J>(x)sin2nnx  = d^. 


For  our  particular  4>(t), 


W1  1 

cn  - 2 f cos2nnxdx  + 2 ^ 4tpCos2nnxdx 


li  Ip 

a sin2nnw.  - sin2nnw. 
nn  1 nn  1 


(li-lp) 

— — — sin2nnw, 
nn  1 
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W1  r1 

= 2 f <f>,sin2Trnxdx  + 2 ■>  4ssin2irnxdx 
o 1 w,  2 


^ i 4*1  4*2  4*2 

cos2irnw1  cos2nnw1 

irn  1 irn  irn  irn  -l 


4>2  _ 4*i 

— ^ — - (cos2Tinw1-l) 


W1  1 

c = f 4>,dx  + f <b0dx 
o o T1  Wj  ^2 

(4>^  ~ 4*2^1  + 4*2 • 

We  can  now  write  <j> ( t ) as  in  (7)  by  noting  that 


, < x _ _ . y „ i 2Trikx  , -2irikXx  , . / 2ui kx  -2nikx» 

4>(x)  = c0  + ck(e  + e ) + dt(e ;_e ) 


where 


= ? a e2irikx 
k=-»  k e 


ck  . dk 

5~  + 7F 


ck  - dk 

T IT 


But  our  method  requires  that  the  series  be  finite.  We  can  choose  an  N 
such  that 


? 

k=N+l 


2-nikx 

e 


+ a 


-k 


e-2"ikxj 


is  small. 


1 


To  use  the  perturbation  method,  we  also  need  £ 4>(x)dx=0.  This 
requires  a transformation  of  A(t)  in  the  following  way.  We  have  now 
that  /*<f>(u)du=a  =c  . Let  ^ ( t ) =<j>( t ) -a  and 

o ° o o 


A(t)  = X0(l+cd»(t) ) 

= A0(l+e(*(t)+a0)) 
= A^l+c^t)) 


where 


e = e/(l+caQ) 

With  the  revised  values  of  Ao,  c and  ^ we  can  now  use  the  formulae  4.8 
and  4.12  to  determine  n(t,A( * ) ) and  p(t,A(*)).  We  can  compute  these  two 
quantities  graphically  with  ll(i,A(*))  and  u ( t , A ( * ) ) . 

The  question  arises,  though,  whether  the  results  will  depend  on  the 
particular  choice  of  c,  Aq  and  , i=l,2.  Let  us  look  at  the  value  of 
Aq.  By  definition 

A = A (1+e.aJ 
o o o 


= A0(1+c(4>1-4>2)w1+c4>2) 
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Recall  that  X , e and  a.  were  chosen  such  that  X.=X  ( l+e4>  ),  1=1,2, 
0 1 10  0 

which  means  also  that  X,-X0aeX  (a.-*,.).  Thus 

1 £ o i Z 

X^  - Xg+f^ j ~X2 )wi  - q X(u)du. 


So  whatever  the  choice  of  parameters,  XQ  will  be  the  same. 

The  other  terms  in  n(x ,X( * ) ) and  p(t ,X ( • ) ) which  involve  our  parame 
ters  are  of  the  form  eXQan,  cXQbn  and  eX  <j>(t).  We  have  that 

k£t<Wj+k 

(^2“^1 ) W1  k+w^<t<k+l 


and 


( 4*  j~4,2  ) ( 1 j ) 

(<^2“^  j )Wj 


k<t<w^+k 

k+Wj<t<k+l 


eX  <j»(t)  = e X (1+ea  )$(t) 
1+ea.  0 0 


= gXq4>(  t) 


(X1-X2)(l-w1) 


k<t<w^+k 

k+w^t<k+l 


Similarly  we  can  show 


■ - ■ — 


c\Q  = (X^-X^)  sin2TrnWj 
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and 


itn 


eX0bn  = (^2-^1)(cos27rnw1-l) 

-Tm 

Thus  we  see  that  our  results  will  be  independent  of  the  choice  of  c, 

XQ,  and  , i=l,2. 

In  Figure  6.1  we  present  the  graphs  of  n(t,X(*))  and  n(t.X(* ))  for 
the  parameters  A^=9,  A2=3,  P=10.  The  solid  line  is  the  graph  produced 
by  the  matrix  approach,  while  the  triangles  form  the  graph  of  the  approx- 
imation. Figure  6.2  displays  the  two  graphs  of  the  mean  waiting  time. 
Again,  the  solid  line  is  the  graph  produced  by  the  matrix  approach,  while 
the  triangles  form  the  graph  of  the  approximation.  To  quantify  the  com- 
parison we  have  calculated  the  average,  maximum  and  minimum  values  of  the 
vertical  distance  between  the  two  graphs  in  each  figure.  These  numbers 
are  based  on  the  values  at  r=k/50,  k=l,2,...50  and  are  found  in  Table 
6. 1 , below. 


Table  6. 1 


ln(T,A(.))-n(T,x(.))l 

|m(t,X(*))-u(t,X(*))| 


Average  Maximum  Minimum 


.00317 

.02622 

.00005 

.00106 

| .00219 

.00060 

It  is  clear  that  in  this  example,  at  least,  the  approximations  ll(i,X(-)) 
and  p(x ,X( • ) ) are  very  good. 
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FIGURE  6.1 


Probability  of  an  Empty  Queue 
*1=9  X2=3  p=10 


A e -Approximation 
- Matrix  Approach 


TAU 


V .... 


6.4  Another  Example:  Periodic  X(t) 

We  would  like  to  be  able  to  judge  the  adequacy  of  our  approximations 
for  situations  where  X(t)  is  not  a step  function.  In  such  cases  the 
matrix  approach  to  the  M(t)/M/1  queue  can  only  give  an  approximation 
itself.  Suppose  we  were  to  approximate  X(t)  above  and  below  by  step 
functions  Xl'(t)  and  Xu(t),  where  X^(t)<X(tH  Xu(t).  We  know  from  Lemma 
3.2.1  that 

n(T,xu(.))in(x,x(*))i  (x ,xL( • ) ) 

and 

u(x,xL(-))si(x,x('))i}i(x,xu(* ))• 

The  matrix  approach , at  the  least,  can  provide  upper  and  lower  bounds  to 
n and  p.  We  also  know  that  the  closer  we  approximate  the  intensity  with 
a step  function,  the  closer  that  approximation  will  be  to  the  "true" 
value.  So  as  we  let  Xu(t)-XL(t)  become  small,  n(x,X(*))  and  p(t,X(*)) 
will  be  squeezed  in  between  their  upper  and  lower  bounds,  and  we  can 
evaluate  the  approximations  fi(x,X(* ))  and  u(x,X(*))* 

We  shall  look  here  at  the  example  in  Chapter  4 where  4>(x)=cos2nx. 

As  mentioned  earlier,  it  is  computationally  efficient  to  choose  equal 
sized  time  intervals  for  the  step  function.  We  let  I(c=^l("^/2N,k/2N) 
where  2N  is  the  number  of  intervals  into  which  we  divide  the  interval 
[0,1).  Thus  we  choose  for  our  step  functions 

XL(t)=Xo(l+c<0^)  , tclj 
XU(t)=X0(l+e4>“)  , tc^ 

A 

! i 

IS 
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where 

= | cos ( 2tt( k-l)/2n) 
l^N+l-k 
and 

4>J^  = Jcos(2Trk/2N) 

1^2N+l-k 

Figures  6.3  and  6.4  display  the  graphs  of  the  probability  of  empti- 
ness and  mean  waiting  time,  respectively,  for  Aq=5,  u=10,  and  e=.5.  The 
solid  lines  are  the  graphs  of  the  upper  and  lower  bounds,  while  the  tri- 
angles form  the  graphs  of  n(t,A(-))  and  y(x, *(•))•  We  see  from  the 
figures  that,  at  least  for  the  approximations,  n and  il  are  quite  good. 
Here  we  have  used  80  intervals  for  the  step  function.  Since  the  average 
difference  between  the  upper  and  lower  bound  is  .0125  for  IT  and  .005  for 
p,  this  seems  to  be  a sufficiently  fine  step  function  approximation. 

We  can  make  a few  comparisons  of  interest.  First  we  look  at  the 
deviation  between  the  perturbation  approximations  and  their  bounds.  It 
is  also  interesting  to  consider  the  mean  waiting  time  and  probability  of 
emptiness  of  a queue  which  has  the  homogeneous  intensity  A - the  average 
intensity  over  a period.  This  is  the  queue  that  is  generally  used  as  a 
substitute  for  the  periodic  queue.  This  can  give  us  a basis  of  compari- 
son to  judge  the  error  our  approximation  gives  us.  In  the  following  two 
tables  we  present  the  average,  maximum  and  minimum  absolute  differences 
between  the  bounds  of  I1(t,A(*))  and  p(t,A(*))  and  the  values  obtained 
from  the  two  queueing  models.  We  let 


k=l,  ...  , N 
k=N+l 2N 

k=l ,2,  ...  , N 
k=N+l , ...  , 2N 
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Dj  = |4>(t,AL(*))-^(t,X0)  1 

D2  = J^(t,XL(*))-^(t,X(*)) | 
D3  = | iP(t  ,XU(  • ) ,XQ)  J 
°4  = |^(t,AU(*))-^t,a(*))|. 

where  4;=n  or  u,  as  appropriate 


Table  6.2  Probability  of  Emptiness 


°1 

°2 

°3 

°4 

Average 

Minimum 

Maximum 

.07334 

.00076 

.12728 

.00687 

.00025 

.01415 

.07219 
.00286 
.11394  1 

.00625 

.00064 

.01243 

Table  6.3 

Mean  Waiting  Time 

D1 

°2 

°3 

°4 

Average 

Minimum 

Maximum 

.02013 

.00052 

.03722 

.00093 

.00033 

.00153 

.01966 

.00012 

.03229 

.00596 

.00509 

.00690 

Since  0^^  and  it  seems  apparent  that  the  approximations  for 

the  periodic  queue  are  better  to  use  than  the  results  for  the  homogeneous 
queue.  It  also  seems  that  n(x,X(*))  and  y(x,X(*))  are  sufficiently  good 
approximations  that  they  should  be  the  choice  over  the  matrix  results, 
particularly  when  ease  of  computation  is  considered. 
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CHAPTER  7 


STORAGE  SYSTEMS  WITH  STOCHASTIC  INPUT  INTENSITY 

7.1  The  Method 

The  input  process  of  a storage  system  will  not  always  be  a deter- 
ministic function  of  time.  It  might  seem  more  realistic  to  assume  that 
the  rate  of  input  varies  randomly  over  time.  This  would  be  the  case,  for 
example,  if  inputs  to  the  system  originate  from  a variety  of  different 
sources  and  the  number  of  sources  is  a random  variable.  We  would  like 
to  be  able  to  examine  processes  such  as  these  in  which  the  input  rate  is 
generated  by  some  stochastic  mechanism.  In  particular,  it  would  be  good 
if  the  method  employed  previously  on  storage  processes  with  nonhomegen- 
eous  Compound  Poisson  input  could  be  adapted  to  allow  for  stochastic 
A ( t ) . Fortunately,  it  is  a relatively  simple  matter  to  do  this. 

As  before,  we  model  our  arrival  intensity  as  A(t)=Ao(l+e<f>(t)).  We 
assume  now,  though,  that  <j>  (t)=<j>(t,w)  is  a measurable  stochastic  process 
such  that  E{<b(t ,oo) } exists  and  / E | cp ( x ,oj)  |dx  <°°,  t>0.  If  we  have  a 

realization  of  the  process,  <f>  (t),  then  A ( t)=  A (l+e<J>u)(t) ) is  just  as 

w o 

before  in  the  deterministic  case.  Thus  the  probability  of  an  empty  store, 
given  that  realization,  is  just  the  conditional  probability 


yt,A(*),W0)  = P{z(t)=0|<ys)  0<s<t} 


(1)  = n(t,A0,W0)+eA0  / <|)(j)(u)fi(t-u,Ao,W1)du+0(e2), 
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t 

As  long  as  E<j>(i(u)  exists  and  f E |<J>( ^(u)  |du  < «>,  the  integral  in  (1)  will 
exist.  We  can  make  this  claim,  since  | fi<t> | < | <J> | implies  that 
Et^(i(u)fi(t-u,Xo,W^) } exists  and  /tE|4)(ii(u)fi(t-u,A^W1)  |du  < <“.  We  note 
that  II.  is  itself  now  a random  variable. 

To  obtain  the  unconditional  probability  we  can  take  the  expectation 
of  11^  with  respect  to  the  distribution  of  Then 


n(t,A(-),Wn)=EP{Z(t)=0|4»(il(s),  0<s<t} 


' t 

= n(t,A0,wo)+cXoE[/  * (u)n(t-uix0,w1)du 


+0{e2) 


By  Fubinis  Theorem  we  can  move  the  expectation  inside  the  integral  to  get 


(2)  n(t,A(-),Wo)=IT(t,Ao,Wo)+  cAq  / m(u)fl(t-u,Xo,W1)du+0(c2) 
where  m(-)  is  the  mean  function. 

We  can  also  find  the  mean  content  of  the  storage  process  by  using 
the  relationship  EZ(t)=E^E(Z(t)  |^w)  j.  If  E{<|>(0(u)p(t-u,A0,Wj) } exists 
and  ^tE^|<|>w(u)p(t-u,X0,W1)  | Jdu  < «,  then 


(3) 


p(t,A(  • ) » Wq ) ~vi ( ^ * X 0 


t ? 

,W0)+eX0  £ m(u)p(t-u,Xo,W1)du+0(c  ) 


It  appears  that  the  complex  behavior  of  the  stochastic  input  process 
is  captured  in  the  mean  of  the  arrival  rate.  If  m(t)  0 then 
ll(t.X(-),Wo)  = Il(t,Xo.Wo)+0(,2)  and  p(  t ,x(  • ) ,Wq)=,i(  t *XQ,Wo)+0(i  2).  Unfor- 
tunately, in  most  situations  of  interest  this  will  not  be  the  case.  We 
note  that  the  formulae  in  (2)  and  (3)  are  now  the  same  as  those  in  Chapter 
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3,  with  m(.)  in  place  of  $(.).  Thus  to  look  at  the  limiting  distribu- 
tion of  n and  u we  are  now  in  the  same  situation  as  that  Chapter. 

It  is  unlikely  that  the  mean  function  will  be  periodic,  so  much  of 
the  results  we  have  looked  at  will  not  apply.  A more  likely  event  is 
that  m(t)=a+b(t)  where  b(t)cL^.  Then 

t 

n(t,A(.),w0)*n(t,x0,w0)+cA0a/  n(t-u,xo,w1)du 

t ? 

+ e\0  f b(u)n(t-u,Ao,Wj)du+0(e  ). 

Since  i is  a bounded  function  which  converges  to  a limit  (which  is  0)  we 
can  use  Lenina  3.2.2  to  show  that 

t A 

lim  £ b(u)n(t-u,A0,W1)du=n(o",A0)^  b(u)du=0 


This  leads  to  the  limit 

<x>  p 

(4)  lim  n(t,\(.),W0)=n (“’.A0)+eVo  fi(u»VVdu+°(c  ) 

t-*on 

7.2  An  Example:  A Queue  with  Two  Phases 

The  most  commonly  studied  process  with  stochastic  arrival  rates  is 
the  queue  in  which  the  arrival  rate  depends  on  the  "phase"  of  the  queue 
at  the  time  of  arrival,  and  the  amount  of  time  spent  in  a phase  is  a 
negative  exponential  random  variable.  This  essentially  yields  an  inten- 
sity A(.)  which  is  a step  function  which  changes  levels  at  random  times. 
This  queueing  process  has  been  discussed  by  Yechiali  and  Naor,  Eisen  and 
Tainiter  and  Neuts,  to  name  a few.  While  the  simple  queueing  problem 
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with  only  two  phases  has  been  studied  extensively  with  traditional 
methods,  it  is  informative  to  look  at  this  example  in  the  light  of  the 
alternative  method  presented  here. 

Suppose  we  have  a situation  where  the  arrival  rate  fluctuates  be- 
tween two  levels  A^A^l+e^)  and  A2=Ao(  1+c4>2 ) » where  the  time  spent  at 
level  i has  a negative  exponential  distribution  with  mean  1/n-.  To  find 
IT  and  u we  need  to  determine  the  mean  of  <p.  We  can  write  m(t)  as 
m(t)=4>jP1(t)+<t>2P2(t)  where 

P.j (t)=P{system  is  in  phase  i at  time  t} 


The  phase  process  forms  a simple  two  state  Markov  Chain.  We  can  use 
the  Kolmogorov  forward  equations  to  determine  P ^.  ( t ) , 1 = 1,2,  as  follows. 

P1(t+dt)=P1(t)(l-n1dt)+P2(t)n2dt+o(dt) 

(5)  P|(t)  = -n1P1(t)+n2P2(t) 
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We  can  now  write  m(t)  as 


m(t)  = if2+(4'1-'l>2)pi(t) 


3 ¥°>) 


nl+n2 


If  we  let  a^+Uj-^^/fnj+r^)  and 


b(t)=-(^1-l|>2) 


n2^(nl+n2^“Pl^°) 


- (n  j-^-rip)  t 

e then  we  see  from  (4)  that 


(7)  limn(t,A( . ) ,W  )=II(<«,A  )+eA 

t-K»  0 00 


<p2+(^l_<p2^n2 


r)l+n2 


OO  ry 

f ll(u,Ao,W1)du+0(t-  ) 


We  recall  from  Chapter  4 that  for  an  M/M/1  queue,  if  WQ=W  then 

£ n(u,xo,w1)du  = -Y1n(~,x0) 

= -d-A0/u) 


= - I 

u 

Substituting  into  (7)  we  have 


(8)  11mn(t,A(.),W) 

t"X» 


4’2+^i-tf’2^l12 

V^“ 


+ 0(e2) 


We  can  use  the  relationship  between  the  parameters  A^ , <t>^  and  r to 


rewrite  this  as 
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(9)  - x - 0(c2) 

We  notice  that,  as  we  would  expect,  the  limit  does  not  depend  on  the 

choice  of  e,  A , <(>.. 

o 1 

It  is  interesting  to  compare  the  approximation  (9)  with  the  actual 
value  as  found  in  the  literature.  Yechiali  and  Naor  (1971)  present  the 
relationship 

(10)  PjqPj  + P2oM2  = 

where  p.Q=P{Z(t)=0  and  system  is  in  phase  i } , i = 1 ,2  in  steau,  e, 
X=P1.X1+P2.X2  » p=pl.pl+p2.p2  and  Pi.=lim  p^(t).  If  U1=Pr= 

^-KJO 

(10)  simplifies  to 

<“>  plo  * »2o  ■ l4/“- 

But  plo+p2O=P(Z(t)=0}.  Thus  we  have  that 

Urn  n(t,An,W)  = 1 - A/p. 
t— 

Notice  that  from  equation  (6)  we  have  that 

(12)  pj  = n2/(n1+n2)  and  p2  =n1/(n1+n2) 

If  we  substitute  (12)  into  (9)  and  rearrange  some  terms,  we  find  that 

0(e2)=0. 

Yechiali  and  Naor  (1971)  state  that  Pi0+P2o=1-^/p  if  and  only  if 
X1^P1*X2^P2*  wh^e  if  1s  true  that  A 1A*i=A2/u2  implies  that 
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W1  -A/m,  we  have  already  shown  in  (11)  that  this  is  not  the  only 
circumstance  under  which  this  relationship  holds.  Their  theorem  does 
indeed  hold  true  when  But  when  the  proof  fails  because  at 

one  point  they  inadvertently  divide  by  zero. 


APPENDIX 


A REVIEW  OF  STANDARD  RESULTS 

The  standard  storage  theory  model  is  the  following:  Inputs  to  the 
system  occur  at  points  tQ,  tj,  . . . where  the  time  between  inputs 
Xi=trti_i  are  i.i.d  random  variables.  The  size  of  each  input,  U. , is 
governed  by  a distribution  B(x)  with  B(0+)=0.  Thus  input  to  the  store 
can  be  described  by  the  sequence  (Xp  Up  Xp  Up  . . In  many  in- 
stances the  X..  are  assumed  to  be  from  a negative  exponential  distribution. 
In  such  cases  the  total  input  into  the  system  is  Compound  Poisson. 

The  content  of  the  store,  in  standard  theory,  is  released  at  a con- 
stant rate  c (r(x)=c)  as  long  as  the  store  is  not  empty.  Thus  in  an 
interval  of  length  y,  the  output  of  the  store  will  be  at  most  c y,  with 
the  maximum  being  achieved  if  the  store  has  positive  content  throughout 
the  interval.  When  the  rate  of  release  is  1,  the  storage  process  is 
equivalent  to  a single  server  queue.  In  this  case  the  Xi  are  interpreted 
as  interarrival  times  and  B(x)  is  the  service  time  distribution.  A 
simple  change  of  time  scale  will  convert  a storage  process  with  r(x)=c 
into  one  with  r(x)=l.  Thus  without  loss  of  generality  we  can  focus  our 
attention  on  the  single  server  queue,  where  much  of  the  literature  is 
found. 

We  cite  here  results  of  the  M/G/l  queue  which  concern  us  in  this 
research.  We  look  to  the  dissertation  of  Smith  (1953)  to  summarize  the 
results  for  this  appendix. 


The  following  notation  will  be  needed: 

G(x)  is  the  distribution  function  of  the  busy  period 
W(x,t)  is  the  distribution  function  of  the  waiting  time  at  t 
ft.  is  the  jth  moment  of  service  time 

J 

u(t)  is  the  mean  of  the  waiting  time  at  t 
X is  the  arrival  intensity. 

This  notation  agrees  with  the  notation  used  in  later  chapters  and  is 
different  from  that  used  originally  in  Smith  (1953). 

The  first  theorem  and  its  corollary  provide  information  about  the 
event  e,  which  occurs  at  t if  Z(t)=0  and  Z(t-)>0,  through  an  examination 
of  the  busy  period. 

Theorem  28:  If  ,\B.<1,  G(x)  possesses  as  many  moments  as  B(x),  and  if  the 
jth  moment  of  B(x)  is  6.  the  moments  v.  of  G(x)  are  given  by 

J J 

Yj  = Bj/d-XBj) 

y2  = Bg/d-XBj)3 

y3  ■ (Sg+axe^/d-xe^Vd-xGj)4 

Corollary  28.1:  Assuming  B(+«')=l  we  have  that 

(i)  If  G(+«’)=l 

(ii)  If  X81>1  • GM=B*(o0) 

where  o0  is  the  root  in  R{s)>0  of  the  equation  s-X(l-B*(s))=0. 

We  see  from  these  theorems  that  the  condition  XB^<1  is  sufficient 
to  ensure  that  the  event  e is  certain.  The  next  theorem  demonstrates 
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conditions  under  which  a limiting  distribution  of  waiting  time  exists. 


Theorem  30:  (i)  If  X3^<1  and  W(°°,0)  =1  then  there  is  proper  distribution 

function  W(x)  such  that 

lim  W(x,t)  = W(x) 
t-*» 

where 

W*(s)  = f e‘sx  dW(x)  = sd-XBj) 

s -X ( 1 -B* ( s ) ) 

(ii)  If  then  for  all  x>0 

lim  W(x,t)  = 0 
t-*» 


We  can  gain  some  additional  information  about  the  queue  which  we 
will  need  from  W*(s).  First,  if  we  take  the  limit  of  W*(s)  as  s goes 
to  infinity  we  find  that  for  the  queue  in  steady  state,  the  probability 
of  an  empty  queue  is  (l-XBj).  Second,  if  we  look  at 

_dw* ( s ) 

lim  — we  ^1’nd  that  the  mean  waiting  time  in  steady  state  is 

S-KO 

X62/2(1-X3j) . These  are  two  well  known  and  often  used  facts. 

The  last  theorem  we  present  is  an  estimate  of  I(t),  the  expected 
idle  time  in  (0,t)  and  p(t). 


Theorem  32:  If  Xp^l,  32<00»  p(0)<  00  then  as  t-*» 

I(t)  = (1-X3x)t  + xe2  - p(0)  + o(l) 
211^X3^ 


p(t)  = X32 

2(1-X3j) 


and 


+ o(l). 
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